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We present a review of recent research on quantum entanglement, with special emphasis on entanglement be- 
tween single atoms, processing of an encoded entanglement and its temporary evolution. Analysis based on the 
density matrix formalism are described. We give a simple description of the entangling procedure and explore 
the role of the environment in creation of entanglement and in disentanglement of atomic systems. A particular 
process we will focus on is spontaneous emission, usually recognized as an irreversible loss of information and 
entanglement encoded in the internal states of the system. We illustrate some certain circumstances where this 
irreversible process can in fact induce entanglement between separated systems. We also show how spontaneous 
emission reveals a competition between the Bell states of a two qubit system that leads to the recently discov- 
ered "sudden" features in the temporal evolution of entanglement. An another problem illustrated in details is a 
deterministic preparation of atoms and atomic ensembles in long-lived stationary squeezed states and entangled 
cluster states. We then determine how to trigger the evolution of the stable entanglement and also address the 
issue of a steered evolution of entanglement between desired pairs of qubits that can be achieved simply by 
varying the parameters of a given system. 
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I. INTRODUCTION 

The studies of quantum entanglement engineering between 
trapped atoms and controlled transmission of information 
from one group of atoms to another are of fundamental im- 
portance and vital to the development of quantum informa- 
tion technologies |1, 2]. In the simplest way, entanglement 
or non-separability can be described as a property of a sys- 
tem of objects that cannot be identified without reference to 
each other - regardless of their spatial separations. In terms of 
quantum-mechanical states of a system of objects, entangled 
states are linear superpositions of the internal states of the sys- 
tem that cannot be separated into product states of the individ- 
ual objects. Entanglement is recognized as entirely quantum- 
mechanical effect which is not only at the heart of the distinc- 
tion between quantum and classical mechanics but also have 
played a crucial role in many discussions on fundamental is- 
sues of quantum mechanics |3- 5]. Apart from the importance 
to understand the fundamentals of quantum physics, entan- 
glement provides a fascinating avenue of research in different 
areas of physics, mathematics and biology. It is now well es- 
tablished that entanglement is an essential ingredient for im- 
plementation of quantum communication and quantum com- 
putation, emerging branches of science and technology based 
on the laws of quantum mechanics. 

The universality of entanglement is that it can be created be- 
tween different objects such as individual photons, atoms, nu- 
clear spins and even between biological living cells. Entangle- 
ment can be used as a resource for transmission of quantum 
information over long distances. However, a transfer of the 
information requires ingredients for the control and transmis- 
sion in the form of long-lived entangled states that are immune 
to environmental noise and decoherence. Theoretical and ex- 
perimental studies have demonstrated that atoms are promis- 
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ing candidates to achieve of all these specific requirements. 
The transfer of information has been accomplished between 
light, which appears as a carrier of the information, and atoms 
whose stable ground states serve as a storage system I^^J^. 
The ability to store quantum information in long-lived atomic 
ground states opens possibilities for a long time controlled 
processing of information and communication of the informa- 
tion on demand. 

A key model for study of entanglement creation, stor- 
age and processing is a system composed of two-level 
atoms— quantum bits, called qubits. A single qubit, often 
called a fuel for quantum information and quantum compu- 
tation, is modeled as a system composed of two energy lev- 
els that are directly or indirectly coupled to each other via 
eg. electric dipole transitions, or transitions through auxilary 
levels. If the qubit interacts with an another qubit, an initial 
excitation can be periodically exchanged between the qubits 
leading to a coherence or quantum correlations being created, 
and then the qubits become entangled. 

An important issue in the creation of entanglement between 
the qubits is the presence of unavoidable decoherence which 
is a source of an irreversible loss of the coherence. It is of- 
ten regarded as the main obstacle in practical implementation 
of coherent effects and entanglement. A typical source of de- 
coherence is spontaneous emission resulting from the interac- 
tion of a system with an external environment. For a com- 
posed system, each part of the system can interact with own 
independent environments or all the parts can be coupled to a 
common environment. In both cases, an initial entanglement 
encoded into the system degradates during the evolution. Nev- 
ertheless, the degradation process can be much slower when 
the parts of the system are coupled to a common rather than 
separate environments and, contrary to our intuition, might 
even entangle initially unentangled qubits. This effect, called 
the environment induced entanglement has been studied for 
discrete (atomic) [9-15], as well as for continuous variable 
systems limflTll . A crucial parameter in the case of the atoms 
coupled to the same external environment is the collective 
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damping | fT8l - l20ll . which results from an incoherent sponta- 
neous exchange of photons between the qubits. The effect of 
the collective damping is to slow down the spontaneous emis- 
sion from the effective collective states of the system ll2lll22ll . 
Under some circumstances, the spontaneous emission can be 
completely eliminated from some of the collective states. The 
states form a sub-space called decoherence-free subspace that 
contains long-lived entangled states that are immune to dis- 
sipation and decoherence ll24l - [26ll . A challenge then is to 
prepare a system of qubits in entangled states which belong 
only to the decoherence-free subspace. 
The destructive effect of spontaneous emission on entangle- 
ment encoded into qubits can take different time scales. The 
decoherence time depends on the damping rate of the state 
in which the entanglement was initially encoded and usually 
the decay process induced by spontaneous emission occurs 
exponentially in time. However, some entangled states can 
have interesting "strange" decoherence properties that an ini- 
tial entanglement can vanish abruptly in a finite time that is 
much shorter than the exponential decoherence time of spon- 
taneous emission 11271 - 129 1. This drastic non-asymptotic fea- 
ture of entanglement has been termed as the "entanglement 
sudden death", and is characteristic of the dynamics of a spe- 
cial class of initial entangled states. The qubits may remain 
separable for the rest of the evolution or entanglement may 
revival after a finite time |30]. The phenomenon of entangle- 
ment sudden death presents a fascinating example of a dynam- 
ical process in which spontaneous emission affects entangle- 
ment and coherence in very different ways. A recent experi- 
ment of Almeida et al. [i3 1,1 with correlated H and ^/-polarized 
photons has shown evidence of the sudden death of entangle- 
ment under the influence of independent environments. The 
requked initial correlations were created by the parametric 
down-conversion process. 

Although the sudden death feature is concerned with the dis- 
entangled properties of spontaneous emission there is an in- 
teresting "sudden" feature in the temporal creation of entan- 
glement from initially independent qubits fsH [lH . The phe- 
nomenon termed as sudden birth of entanglement, as it is 
opposite to the sudden death of entanglement arises dynami- 
cally during the spontaneous evolution of an initially separate 
qubits. The sudden birth of entanglement is now intensively 
studied as it would provide a resource for a controled creation 
of entanglement on demand in the presence of a dissipative 
environment ll34ll . 

Entanglement and coherence effects can be preserved from 
the destructive spontaneous emission by placing a system of 
atoms inside a cavity or by sending atoms across the cavity 
field. Cavity Quantum Electrodynamics (CQED), the study 
of the interaction of atoms with the cavity modified electro- 
magnetic field, provides a different context to explore entan- 
glement properties and its dynamics. Different types of cav- 
ities have been proposed and implemented raging from mi- 
crowave (Hit], and high-Q optical cavities ll36[ IstIi . to cavities 
engineered inside photonic band gaps materials pSi]. When 
atoms are placed inside or are slowly passing through a cav- 
ity, they interact with a modified electromagnetic field that is 
limited in space to a finite solid angle centered around the cav- 



ity axis. As a result, their radiative properties change, so that 
spontaneous emission can be modified or even suppressed. 
The suppression arises from a reduction of the number of field 
modes to which atoms can be coupled. With the help of a cav- 
ity, entanglement of atomic systems can be tested by probing 
the atomic states of the atoms lea ving the cavity, or by a mea- 
surement of the cavity photons f24', '39']. For example, the 
absence of fluorescence can be a signature for maximally en- 
tangled atom pairs which, in principle, can be used to create 
cluster states for one-way quantum computing fl40l] . 
An another scheme based on interactions in CQED has been 
proposed that uses a combination of the cavity field with 
atomic Raman transitions between a pair of long-lived atomic 
ground states ]41-45]. The Raman transitions are mediated 
by two laser fields that are highly detuned from the transitions 
between the atomic ground and the excited states. The tran- 
sitions occur through a virtual level without population of the 
atomic excited states, thereby avoiding spontaneous emission. 
In the CQED systems there is an another loss mechanism, the 
dissipation of the cavity mode though imperfections of the 
cavity mirrors. Although the cavity loss cannot be avoided 
in a similar fashion to spontaneous emission, it is possible 
to achieve conditions in which the dissipation of the cavity 
mode could be constructive rather than destructive for entan- 
glement creation and protection from decoherence. In par- 
ticular, the cavity damping can cause an atomic system to 
decay to decoherence-free entangled states. For example, it 
is well known that cavity QED measurements on photons 
can be an effective tool for manipulating and observing the 
states of atom-cavity systems and have applications in quan- 
tum feedback ll46ll . Cavity QED has an arsenal of detection 
methods such as homodyne, intensity second-order correla- 
tion, intensity-field third-order correlation techniques which 
available to use and has been successfully applied for quan- 
tum dynamics of atoms interacting with a cavity field. 
Bearing in mind all of those difficulties with the elimination of 
spontaneous emission, we may conclude that the best for en- 
tanglement prevention from decoherence is to encode it into 
states which belong to a group of states contained inside a 
decoherence-free subspace. However, this creates a problem 
of how to access the internal states of the subspace to trigger 
an evolution of the encoded information. It has been proposed 
to implement a dynamical technique of addressing individual 
atoms by well stabilized laser fields to achieve controlled ma- 
nipulation of the internal states of the atoms and evolution 
of entanglement ll47ll48ll . However, under realistic conditions 
the problem of addressing individual atoms poses a substantial 
experimental challenge. Thus, more hopeful systems from the 
experimental point of view have been proposed where the evo- 
lution of an initially stable entanglement could be triggered 
by a change of internal parameter of the system ll49ll . The pa- 
rameters that could be changed are, for example, the coupling 
constants between the atoms and the field modes. Alterna- 
tively, it could be done by changing frequencies of the cavity 
modes to detune them from the atomic resonance frequencies. 
These ideas could then be applied to a controlled steering of 
the evolution of an initial entanglement to a desired pair of 
qubits including a situation where qubits are completely iso- 
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lated, such as completely separated cavities each containing a 
single atom. 

Apart from two qubit systems, various attempts have been 
made to produce entanglement in more complex systems. An 
extensive research have been done on generation of multipar- 
tite continuous variable entangled states of atomic ensembles 
form multi-mode squeezed light Il50l - l53ll . The reason for us- 
ing atomic ensembles is twofold. On the one hand, atomic 
ensembles are macroscopic systems that are easily created in 
the laboratory. On the other hand, the collective behavior 
of the atoms enables to achieve almost a perfect coupling of 
the ensembles to external squeezed fields without the need to 
achieve a strong field-single atom coupling. Particularly inter- 
esting are schemes involving ensembles of hot or cold atoms 
trapped inside a high-Q ring cavity lisil - lsill . The schemes are 
based on a suitable driving of four-level atoms with two ex- 
ternal laser fields and coupling to a damped cavity mode that 
prepares the atoms in a pure entangled state. Similar schemes 
have been proposed to realize an effective Dicke model oper- 
ating in the phase transition regime ll43ll . to create a station- 
ary subradiant state in an ultracold atomic gas (sf], and to 
prepare trapped and cooled ions in pure entangled vibrational 
states ifssll . Simple procedures have been proposed to pre- 
pare ensembles of cold atoms in pure multi-mode entangled 
states and ensembles of hot atoms in pure entangled cluster 
states m. 

In this paper we review the major processes for entanglement 
creation and processing of entangled states in atomic systems. 
The plan of this review is as follows. We begin in Sec. 2 with 
an introduction of the concept of quantum entanglement and 
some of its features. In Sec. 3, we give a brief description 
of the formalism we adopt to overview the research on en- 
tanglement and disentanglement. The formalism is based en- 
tirely on the density matrix representation of systems in mixed 
states. We then in Sec. 4 introduce measures of entanglement. 
We give a short review of measures of entanglement between 
two qubits, but we mostly focus on concurrence, the neces- 
sary and sufficient conditions for a two-qubit entanglement. 
In Sec. 5, we introduce the elementary models for study en- 
tanglement. We specialize to two-atom systems and pay par- 
ticular attention to the dynamics of atoms interacting with a 
common environment. In a series of examples, we examine 
specific systems and processes for entanglement creation be- 
tween two atoms. In Sec. 6, we employ the dynamics of the 
atoms to illustrate some "strange" temporary behaviors of the 
disentanglement process, the sudden death and sudden birth of 
entanglement. We demonstrate under what kind of conditions 
the entanglement could revival during the evolution. Then, we 
establish a connection between sudden birth and the presence 
of long-lived states of the two-atom system. In Sees. 7 and 
8, we specialize to a cavity QED situation, where two atoms 
could be located inside a single standing-wave cavity or in 
two completely isolated cavities. The cases of a good and bad 
cavity are investigated and the role of an imperfect match- 
ing of the atoms to the cavity modes is discussed. Next, in 
Sec. 9, we tackle the problem of triggered evolution of entan- 
glement. Our focus is on how one could trigger an evolution 
of a stabile or sometimes referred to as "locked" or "frozen" 



entangled state without the presence of any external fields. We 
consider three practical schemes where the atoms are coupled 
to a single mode cavity field or are subject of the interaction 
with an environment whose the internal modes are in their 
vacuum states. In Sec. 10, we establish how one could steer 
the evolution of entanglement to a desired state. The role of 
the rotating-wave (RWA) approximation in the evolution of an 
initial single excitation entangled state is discussed in Sec. 1 1 . 
In the RWA, the so-called counter rotation terms in the atom- 
field interaction Hamiltonian are ignored, and in this section 
we inquire how and under which circumstances, the time evo- 
lution of the initial entanglement could be modified by the 
counter-rotation terms. The remaining two sections. Sees. 12 
and 13, are devoted to discuss the recently proposed proce- 
dures for deterministic creation of pure continuous variable 
entangled states and cluster states in atomic ensembles located 
inside a high-Q ring cavity. 

II. QUANTUM ENTANGLEMENT: A SIMPLE PICTURE 

We devote this introductory section to an elementary discus- 
sion of the concept of quantum entanglement or simply non- 
separability of systems composed of qubits— quantum bits. A 
qubit may be simply a single atom for example, or a quan- 
tum dot, or even a large, complicated molecule with a single 
transition within its nondegenerate energy levels. For the pur- 
pose of this review, we restrict ourselves to two-qubit systems 
only, which is the simplest system for study entanglement. We 
shall characterize entanglement in terms of pure and mixed 
states of a given system, and explore various properties of en- 
tanglement, which have a particular significance in describing 
unusual behaviors of entanglement and lead to further under- 
standing of the notion of entanglement. 
Our focus is on a system composed of two quantum subsys- 
tems, two qubits labeled by the suffices A and B, and prepared 
in a pure quantum state described by the state vector \'^). We 
adopt the following definition of entanglement. The qubits are 
said to be entangled or non-separable, if the state vector 
cannot be factorized into the tensor product of the state vec- 
tors I^'a) and \'^b) of the individual qubits 

I*) ^ 1^-^)0 1*5). (1) 

Otherwise, when the factorization is possible, the qubits are 
separable. An example of entangled state of two qubits is a 
superposition state 

1^-) = a|ei) ® I52) +6|5i) ® |e2), (2) 

where |ei) (X) \g2) and \gi) (g) |e2) are the basis states of the 
Hilbert space of the system and \ei) and \gi) represent respec- 
tively the excited and ground state of the ith qubit. The pa- 
rameters a and b which, in general are complex numbers and 
|ap + |6p = 1, determine the degree of superposition of the 
two product states. When a 7^ 5, we say that the state (|2| is 
a non-maximally entangled state, and reduces to a maximally 
entangled state when a — b. 

In order to develop our discussion to the entangled properties 
of a system consisting of two parts and prepared in a mixed 
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state, we consider a pair of qubits prepared in a quantum state 
described by the density matrix p. The reduced density matrix 
describing the properties of the A{B) subsystem is obtained 
by tracing the total density matrix p over the B{A) subsystem 

PA = Tr_Bp, pB = TrAp. (3) 

If the two-qubit density matrix p can be written as the tensor 
product of the density matrices of the individual qubits 

P = PA®PB, (4) 

then the qubits are separable. If such a factorization is not 
possible, then the qubits are entangled. 
For further understanding of entanglement, consider some 
characteristic properties and fundamental features of entan- 
glement. It is often asserted that entanglement is a quan- 
tum phenomenon. By this it is meant that entanglement is a 
property resulting from the superposition principle, the basic 
principle of quantum mechanics. The characteristic property 
of entanglement is the impossibility to distinguish between 
two objects prepared in a maximally entangled state. This is 
manifested in a simultaneous measurement of the state of the 
qubits. If the qubits are prepared in the state (|2|i with a = b, 
then the probability of finding the qubit A in the excited state 
and simultaneously the qubit B in the ground state is equal 
to 1/2, and the same results is found for the probability of 
finding the qubit A in the ground state and simultaneously the 
qubit B in the excited state. 

A fundamental feature of entanglement is that it can easily 
achieve the optimal value if a system of qubits is prepared in 
a pure state, and rapidly degrades when it is shared between 
different states or when the qubits are coupled to an another 
system of qubits or to a surrounding environment. It may re- 
sult in a mixed state of the system. We can easy understand 
that property just by analyzing a simple example of entan- 
glement sharing between two entangled states. Consider two 
Bell states, the prototypes of two-qubit maximally entangled 
states fislin 

l*.)-^(|ei)®l52) + l5i)®|e2)), 

|vl/,> = i=(|ei)0|52>-|gi)®|e2)). (5) 

The states (|5]i are symmetric and antisymmetric linear su- 
perpositions of the product states of individual atomic states 
which cannot be separated into product states of the individual 
atoms. 

If the system of the two atoms is prepared in one of these 
states, say then the atoms are maximally entangled. If, 
on other hand, the system is prepared in an equal superposi- 
tion of these two states then the resulting state, that is a super- 
position of these two states is 

|T> = ^(|vI/.) + |*a» = |ei)§5|<?2), (6) 

which is a product state of the individual qubits. Thus, this 
simple example clearly shows that the process of sharing of an 



entanglement between different entangled states can diminish 
entanglement or even can result in the complete separability 
of the qubits. 

Entanglement sharing may take place not only among states 
of the same system but also among combined states of one 
system and another, e.g. an another qubit or a set of qubits. 
The entanglement may result from the direct or indirect inter- 
action between them. We will deal with this issue latter in the 
Sec. 5. 



III. DENSITY OPERATOR REPRESENTATION 

Traditionally, the concept of entanglement in an atomic sys- 
tem is introduced and understood in terms of the atomic den- 
sity operator. Here, we provide a general discussion of the use 
of the density operator formalism in the description of entan- 
gled systems. We discuss general properties of the density op- 
erator and the basic conditions that the density operator must 
satisfy. Then, we introduce different orthogonal basis that are 
employed in the description of the density operator, using a 
two-atom system as an illustration. 

The possible states of an atomic system are specified by the 
atomic density (matrix) operator p. The operator p of any 
quantum system must satisfy the three basic conditions: 

1. p is Hermitian, p = p\ 

2. p has unit trace, Trp = 1, 

3. p is positive. Thus for any state \^): {'^\p\^>) > 0. 

It is useful to consider the eigenvalues and normalized 
eigenvectors \Xi) of p: 

p\X,)=K\X,), p = ^A,|A,)(A,|. (7) 

i 

Then, the two basic conditions hold as follows 

1. < A, < 1, 

Other properties also follow 

• The determinant of p and all principle minors are posi- 
tive. 

•Trp2 = E.Af <1. 

• For any state |*): (*|p|*) < 1. 

A possible situation is to have all the eigenvalues of p equal to 
zero, except one which is unity. In this case, the system is in 
a pure state characterized by 

p=|A)(A|, p = p2 and Trp^ = 1. (8) 

The subjects of entanglement generation and its dynamics can 
be studied in any complete set of basis states of a given sys- 
tem. For example, in the representation of an orthogonal ba- 
sis, the density matrix for any quantum system can be 
defined by its matrix elements 

P^J = mP\j), (9) 
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where the states \i) and \j) satisfy the orthonormahty and 
completeness relations, — Sij and = 1' 

spectively. 

Usually, we choose the basis that provide equations of mo- 
tions for the density matrix elements in a simple mathemat- 
ical structure that allows for a particularly transparent phys- 
ical interpretation. Moreover, in order to get a better inside 
into physics involved in a specific process, we often change 
the analysis to another orthogonal basis, which is done via a 
transformation 



a 



(10) 



where is transformation coefficients satisfy the normal- 
ization condition |CmP — 1- 

In terms of the new basis states, the density matrix elements 
are given as 



(11) 



which shows that the new density matrix elements can be 
easily related to the old density matrix elements through the 
transformation coefficients. 

It is natural to choose a set of orthogonal states as basis for 
the representation of the density operator. If we deal with a 
finite-dimensional system consisting of two two-level atoms, 
we often determine the density matrix of the system in the 
basis spanned by four product (separable) states 



= ® 1.92), |^'2) = |<7i)«)|e2), 
l*3> = |ei) ® 132), 1*4) = |ei) «) lea). 



(12) 



Of course, the form of the density matrix written in the ba- 
sis ( fTSl i will depend on a specific process and on the detailed 
dynamics involved. The density matrix can have a diagonal 
or non-diagonal form and the presence of any non- diagonal 
terms simply indicates the existence of coherence effects in 
the system. The coherence effects are essentially a manifesta- 
tion of the correlation which may exist between the atoms. 
Various type of coherences can be distinguished depending on 
which of the density matrix elements are nonzero. In general, 
the density matrix of a system of two atoms written in the 
basis (fTSl i is of the form 



(13) 



In all cases we will refer to the atomic density matrix el- 
ements as follows. The one-photon atomic coherences are 

Pl2,Pl3,P21,P31,P23,/032,P42,/024,P34,P43, the tWO-photon 

atomic coherences are pi4, pn, and the populations of the 
states are pn, P22, P33, Pa- Among the populations, we dis- 
tinguish the population pn of the ground state with zero pho- 
tons, excited one-photon states with populations P22, P33, and 
a two-photon excited state with the population P44. The ap- 
propriateness of this terminology describing atomic density 
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matrix elements follows from the examination of the number 
of photons involved in each state and transitions between the 
states. 

For majority of situations considered in this review, the den- 
sity matrix of a two-atom system will occur in a simplified 
form, the so-called X-state form 115911 
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where non-zero matrix elements occur only along the main 
diagonal and anti-diagonal. Physically, the X-state form cor- 
responds to a situation where all one-photon coherences be- 
tween the ground and one-photon states and between one- and 
two-photon states are zero. The X-state density matrix can be 
easily created by an appropriate initial preparation of a two- 
atom system, e.g. by the pumping of the atoms by correlated 
(squeezed) light beams as, for example, can be obtained from 
the output of a parametric amplifier Also, one can find pro- 
cesses that not only retain the initial X-state form of the den- 
sity matrix, but even could lead to the X-state form under the 
evolution. 

A variety of situations considered here will include the cou- 
pling of atoms to a common environment (reservoir). In this 
case, the product states (fTSl i do not correspond in general to 
the eigenstates of the system. A different choice of basis 
states is found particularly useful to work with, the basis of 
collective states of the system, or the Dicke states, defined 

as ifilfil 



I*i) = |gi)®l52), 

I*s) = ^(|ei)®|ff2) + |5i) 

I*a)-^(|ei)®l52)-|gi) 
1*4) = |ei) «) |e2). 



\e2)), 

162)), 



(15) 



We see that the collective basis contains two states, j^I'^) and 
l^'a) that are linear symmetric and antisymmetric superposi- 
tions of the product states, respectively. The most important is 
that the states are in the form of maximally entangled states. 
As we shall see in Sec. 5, the entangled states are created by 
the interaction of the atoms through the dipole-dipole poten- 
tial induced by the coupling of the atoms to the same environ- 
ment. The density matrix, if diagonal in the collective basis 
indicates that the dynamics of the states are then independent 
of each other Nevertheless, we will see that under some cir- 
cumstances an entanglement, initially encoded in one of the 
entangled states may not evolve independently of the evolu- 
tion of the remaining states. 

The Dicke states are found as a useful basis in situations where 
only one-photon coherences are present. When two-photon 
coherences are involved, it is often found convenient to work 
in the Bell state basis that is composed of four maximally en- 
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tangled states 



is the von Neumann entropy, and 
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(l.9i) ® 



152/ + Iffi) 

I52) - |.gi) 

I52) + |ei) 
I52) - |ei) 



|e2)), 
|e2)). 



(19) 
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We shall refer to the states |5's) and \'i'a) as one-photon Bell 
states, since they are characterized by states involving a single 
photon, and to the states |$s) and as two-photon Bell 
states, since they are characterized by states involving two 
photons. The states |5's) and jvPa) are often called Bell states 
with anti-correlated spins, whereas and |$a) are called 
Bell states with correlated spins. The advantage of working in 
the Bell state basis is the most evident in the case when two- 
photon coherences pi4 and p4i are different from zero. In this 
case, the density matrix, when written in the basis formed by 
the Bell states may be represented in the diagonal form. It is of 
interest to mention that the transformation is also useful even 
if the resulting transformed matrix is still in a non-diagonal 
form. In this situation, one could find the matrix useful, for 
example, for analyzes of the departure of the states of a given 
system from the Bell states. 

In closing this section, let us conclude that a typical reason of 
making a transformation from one basis to another is in the 
simplification of the density matrix involved. A convenient 
change of the basis may transfer the density matrix from a 
non-diagonal to diagonal form. 



IV. MEASURES OF ENTANGLEMENT 

In order to determine the amount of entanglement for a pure 
or a general mixed state of quantum systems of an arbitrary 
dimension, it is essential to have an appropriate measure of 
entanglement. A "good" entanglement measure should vary 
from zero, for separable states, to unity for maximally entan- 
gled states 1 60]. It is not difficult to give a good entanglement 
measure for a system in a pure state |61-66]. The classic ex- 
ample of pure states for study entanglement are the Bell states 
that represent a class of maximally entangled states, and thus 
a good measure of entanglement should give unity for such 
states. 

A measure, which satisfies the requirements for being a good 
measure for entanglement is the pure state entropy of entan- 
glement 



(17) 



where | ^P) is a pure state of the system composed of two sub- 
systems A and B, 



are reduced density operators of the subsystems. 
The quantity which determines the amount of entan- 

glement in a pure state |^') ranges from zero for a product 
state to maximum E{^) = 1 for a Bell state. 
As is well known, the entropy is a measure of the disorder in 
a system and obviously depends on the correlation properties. 
It is probably for this reason that the measure of entropy has 
played an important role in the development of the theory of 
entanglement measures. 

The problem of defining a good measure of entanglement is 
more complex for a system in a mixed state. There are some 
difficulties with ordering the states according to various en- 
tanglement measures that different measures can give differ- 
ent orderings of pairs of mixed states. Therefore, there is a 
problem of the definition of the maximally entangled mixed 
state. 

A number of measures of entanglement of mixed states have 
been proposed among which concurrence and negativity sat- 
isfy the necessary conditions for being good measures of en- 
tanglement I67-69J. Both measures are unity for Bell states 
and zero for product states. For intermediate values of concur- 
rence or negativity, there is some entanglement in the system, 
and the value of either of the two measures can be treated as a 
degree of entanglement. Concurrence and negativity are now 
widely accepted measures of entanglement. 
For the purpose of this review article, we adopt concurrence 
as the measure to quantify entanglement between qubits. The 
advantage of using concurrence is that it is relatively simple to 
calculate. It involves a standard procedure of multiplication of 
matrices. We begin with a brief outline of the approach for the 
calculation of concurrence. Next, we employ the concurrence 
to derive general expressions for entanglement of a class of 
physical systems determined by a common density matrix. 
For any system composed of two qubits and determined by the 
density operator p, the entanglement of formation is defined 
by 



£ (p) = inf 



E 



(20) 



where < p„ < 1 is a probability distribution and the infi- 
mum is taken over all pure state decompositions of pf^. The 
entanglement of formation for a mixed quantum system can 
be written in terms of the Shannon entropy and concurrence 

as 



e{p) = 



^Q + ^Vl-C^ (P)), (21) 



where H{a) is the Shannon entropy and C (p) is called con- 
currence. 

The concurrence can be calculated explicitly from p as fol- 
lows iiSi 



-Tr (p, log2 p^) , 



(18) 



C{p) ~ max ^0, \/Xi 



/A4 



(22) 
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where the quantities A; are the the eigenvalues in decreasing 
order of the matrix 



R = P{<jy® (Ty) p* {ay ^ ay) 



(23) 



where p* denotes the complex conjugate of p given in the ba- 
sis of the product states (fTSl i and ay is the Pauli matrix ex- 
pressed in the same basis as 



-i 

1 



(24) 



Concurrence varies from C = for separable qubits to C = 1 
for maximally entangled qubits, and the intermediate cases 
< C < 1 characterize a partly entangled qubits. 
The concurrence is specified by the density matrix of a given 
system and thus can be determined from the knowledge of the 
density matrix elements. In a general case such as described 
by the density matrix ( fTST l. the concurrence appears to have 
a very complicated form, not convenient in analytical discus- 
sion. The difficulty consists in obtaining the roots of the fourth 
order polynomial that determine the four eigenvalues of the 
density matrix. In the following, we consider some specific 
situations giving simple analytical expressions for the concur- 
rence suitable for discussion and interpretation. A particularly 
simple expression for the concurrence is obtained when the 
density matrix of a given system, written in the vector space 
spanned by the product states (fT2l l. has a block diagonal form 



Pit) 



Piiit) 

P22{t) P23{t) 

P32{t) P33{t) 

pM{t) 



(25) 



in which all the coherences except P23{t) and P32{t) are equal 
to zero. In practice, it could correspond to a situation of two 
qubits coupled to a noisy reservoir with no initial and external 
coherences. In this simplified case, the concurrence can be 
easily computed as 



C{t) = 2 max {O, \p23m-y/pii{t) Piiit)] 



(26) 



It is seen that a non-zero coherence between the | ^'2 ) and | ^'3) 
states is the necessary condition for entanglement, but not in 
general sufficient one since there is also a threshold term in 
the concurrence, as seen from Eq. ( |23] |. involving the popu- 
lations /Oii(t) and pii{t). For some situations, the quantity 
\/ Piiit)p4:iit) will be different from zero, and \p23\ may be 
positive but not large enough to enhance the concurrence C 
above the threshold for entanglement. Thus, the necessary 
and sufficient condition for entanglement between two qubits 
whose the dynamics are determined by the density matrix (|25] | 
is 



|P23WI^0 and \p23{t)\ > ^ piiit)pii{t). (27) 

We can see immediately that the threshold behavior depends 
on the population of the two-photon state |5'4). Thus, no 
threshold features can be observed if the entanglement cre- 
ation and the evolution of the system involves a single photon 



only, which rules out the possibility to populate the state \ ^i). 
Nevertheless, the threshold behavior still could be possible to 
observe. It will be seen by explicit calculations, which we 
will do very soon, that the non-Rotating Wave Approximation 
effects will give some interesting results on the threshold be- 
havior of entanglement even if only a single photon is present 
initially in the system. 

The threshold for entanglement is a sort of a "boundry" be- 
tween classical and quantum behavior of a two-qubit system. 
Above the threshold, the behavior of the system is determined 
in terms of superposition states, the basic principe of quan- 
tum mechanics. Below the threshold, the qubits are separable 
and no superposition principle is required to determined their 
properties. 

A particularly interesting behavior of concurrence is found for 
a system of qubits determined by the density matrix of the X- 
state form, Eq. ( fT4] i. Note that the X-state form arises natu- 
rally in a variety of practical situations that we will discuss in 
details in the following sections. First, we find that the square 
roots of the eigenvalues of the matrix R are 




|Pl4 



\P23 



{t)\ ± v/pilWP44(i), 



(28) 



from which it is easily verified that for a particular value of 
the matrix elements there are two possibilities for the largest 
eigenvalue, either \/Xi or VAa- The two possibilities result in 
the concurrence of the form 



C{t) =2max{0, Ci(t), C2(t)} 



where 



and 



Ci{t)^\puit)\~ ^P22it)p33it), 



C2W = \p23it)\ - Vpllit) Piiit). 



(29) 



(30) 



(31) 



From this it is clear that the concurrence C{t) can always be 
regarded as being made up of the sum of nonnegative contri- 
butions, Ci{t) and C2(t) associated with two different coher- 
ences that can be generated in a two qubit system. From the 
forms of Ci{t) and C2(i), it is obvious that Ci{t) provides a 
measure of an entanglement produced by the two-photon co- 
herence pii, whereas C2(i) provides a measure of an entan- 
glement produced by the one-photon coherence p23- The con- 
tributions Ci{t) and C2(t) are traditionally understood as the 
criteria for one and two-photon entanglement, respectively. 
They are often identified as weight functions for the distribu- 
tion of entanglement between the two coherences. 
Inspection of Eq. (fl4l) shows that the two separate contri- 
butions to the concurrence are associated with two separate 
blocks the density matrix elements are grouped. We can dis- 
tinguish an inner block composed of the one-photon popu- 
lations, p22,P33 together with the coherences p23,P32, and 
an outer block composed of the populations pn , P44 together 
with the two-photon coherences pi4,pii. Further analysis 
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show that the inner block could be related to the one-photon 
Bell states, whereas the outer to the two-photon Bell states. 
Thus, the threshold behaviors of the concurrence, seen in 
Eqs. (l30t and (ISTT l. can be interpreted as a competition be- 
tween one and two-photon processes or equivalently between 
one and two-photon Bell states in the creation of entangle- 
ment of the X-state form in the two-qubit system. If an entan- 
glement is created among the one-photon states, the amount 
of entanglement created is limited by a population of the 
two-photon states, and vice versa, if entanglement is created 
among the two-photon states, the entanglement is limited by 
a population of the one-photon states. 

The previous discussion on the concurrence focused on the 
detailed analysis of the relations between the two different 
groups the density matrix elements could be divided, with- 
out specifying the state of the system. There are specific cases 
of interest, where two qubits are found in a pure entangled 
state, or in a mixed state involving only a single or two en- 
tangled states. It is worth discussing these specific cases be- 
cause properties of the concurrence can be easily determined 
from the nature of the entangled states and their populations. 
For example, consider a system of two qubits prepared in a 
pure non-maximally entangled state of the form of the wave 
function (|2]l. This state is associated with the one-photon co- 
herence and could be created between the states of the inner 
block of the matrix ( fT4] i. Note that in the special case of 
\a\ = \b\ = 1/^/2, the state coincides with the Dicke state 
Ivts). We can write the state Q in an equivalent form 



1 



(|ei)® |g2) + a|5i)® |e2)), (32) 



where a — b/a. A pure-state density matrix formed from the 
state ( l32l i is of the form 



1 



P = 



(33) 



^1 + |q|2 T+H^ 



and leads to a particularly simple expression for the concur- 
rence 



C = 



2\a\ 



1+ a 



(34) 



The concurrence involves only a cross term, the product of the 
off-diagonal elements of the density matrix. Note an interest- 
ing feature of the concurrence ( l34l l. namely, the concurrence 
of the pure state is determined only by a single parameter \a\, 
the ratio of the probability amplitudes of the product states. 
If, in addition to the entangled state, the Hilbert space of the 
system is expanded to include few separable states, then the 
expanded system could be found in a mixed state. Neverthe- 
less, the concurrence still could be determined by the param- 
eters characteristic of the entangled state. A consequence of 
the presence of the separable states could be in a possibility 
of the threshold behavior for entanglement. A familiar ex- 
ample of this situation is the Dicke model, which refers to 
two atoms coupled to a common reservoir and confined to 



a region much smaller than the resonance wavelength of the 
atomic transitions ifTsl [T9I1 . The Hilbert space of the Dicke 
model is spanned by three collective state vectors, l^fi), |^'s) 
and 1 5*4). In the absence of the coherences, the density matrix 
of the system, written in the basis of the three state vectors, 
takes a diagonal form 



Pii 

P = \ pss 
P44 

and then the concurrence is given by 



C = 



:{0, Pss - 2^piip44} . 



(35) 



(36) 



We see that positive values of the concurrence are brought by 
a non-zero population of the entangled state and unlike the 
pure state situation, the concurrence is not determined by a 
single parameter The involvement of the other density matrix 
elements can lead to a threshold behavior for entanglement. 
It happens whenever the term 2^^11/344 overweights the pop- 
ulation Pss of the state It is interesting to note that the 
threshold could occur only if the population is shared between 
the entangled state and both of the separable states. If only one 
of the separable states is populated, no threshold for entangle- 
ment is observed. 

A further very instructive example is a situation when the 
Hilbert space of the system is spanned by state vectors among 
which two states are entangled and the other are separable. An 
example of a physical system corresponding to this situation is 
the case of two qubits coupled to the same environment. The 
Hilbert space of this system is composed of the four collective 
states ( fTSl ), and the resulting density has the following form 



P - 



/Pii 

Pss Psa 

Pas Paa 

V P44 



(37) 



With the density matrix dJTb . the concurrence is of the form 

C==max{0, C2}, (38) 

with 



Co 



\J [Pss - Paaf' ~ [psa - Pasf - 2^piiP44. (39) 



It is seen that the criterion for entanglement becomes very dif- 
ferent when two or more entangled states are present. One 
would expect that a larger number of entangled states involved 
in the description of a system should result in a better entan- 
glement. This is not the case. According to Eq. ( [39l ), the 
concurrence depends crucially on the population difference 
rather than the sum of the populations of the two entangled 
states involved. A further analysis of Eq. ( |39] | results in two 
important conclusions. First of all, any two qubis represented 
by two maximally entangled states, which are equally popu- 
lated, are completely separable. This is true independent of 
whether the entangled states are correlated or not. Hence, two 
equally populated Dicke entangled states will not give rise to 
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entanglement even though they may be strongly correlated. 
Secondly, a close look at Eq. ( l39b reveals that the correlations 
between the entangled states have a further destructive rather 
than constructive effect on entanglement. This is in contrast 
to the correlations between product states that always appear 
as constructive correlations for entanglement. 
Finally it may be mentioned, that threshold behavior of entan- 
glement is also possible in this system, as indicated by the 
presence of the term 2^piip44 in Eq. ( [39] l. Thus, similar 
to the above example, a threshold of entanglement could be 
observed only if the redistribution of the population involves 
both of the separable states. 

The examples presented in this section show that we can find 
a variety of physical systems which exhibit the same entan- 
gled properties. We have established the relation between the 
concurrence and the form of the density matrix of a physical 
system. We hope that the general situations considered here 
will lead to a better understanding of the detailed dynamics 
of entanglement discussed in the succeeding sections, where 
we specialize to specific situations of two atoms and atomic 
ensembles evolving under the influence of noisy environment 
and/or external fields. 



V. ELEMENTARY MODELS FOR ENTANGLEMENT 

The theoretical study of multi-atom systems is generally com- 
plicated and the physics involved is often difficult to grasp in 
the necessarily complicated analysis. Although less impor- 
tant from the point of view of practical applications, studying 
simpler systems such as single two-level atoms is significant 
in terms of enabling analytic treatments to be undertaken and 
thus facilitates a better insight into the fundamentals of the 
atomic dynamics. It is our purpose in this section to special- 
ize the preceding general discussion on entanglement creation 
and processing to the simplest examples of potential atomic 
systems for entanglement. Four cases are studied: (1) Two 
atoms in free space and interacting with a common environ- 
ment, (2) two atoms interacting with a single-mode cavity 
field, (3) two atoms located in separate cavities and interacting 
with local single-mode fields, (4) atomic ensembles interact- 
ing with radiation modes of a high-Q ring cavity. 
The systems considered here involve a set of identical atoms 
located at fixed positions r^, and separated by a distance 
'''ij — l^j ~ '^il l^ge compared to the atomic diameter, so 
that overlap between the atoms can be ignored. The atoms 
are modeled as two-level systems with ground states \gi) and 
excited states \ei), separated by transition frequencies cjq and 
connected by a transition dipole moment jl. Both atoms are 
assumed to be damped with the same rates 7— the sponta- 
neous emission rates arising from the coupling of the atoms to 
a common environment. The atoms can be prepared initially 
in an arbitrary state that could be a separable (product) or a 
superposition (entangled) state. The initial state may evolve 
in time under the action of the Hamiltonian of the system, or 
its evolution may be "frozen" for all times. 
The dynamics of the atoms can be studied in a complete set 
of basis states of a given system. For the first scheme of two 



atoms interacting with a common environment, we will work 
in the basis of the Dicke states ( fTSb . For the second scheme 
of two atoms located inside a single-mode cavity field, it will 
prove convenient to study the dynamics of the atoms in the ba- 
sis of four product states ( fT3] l. The third scheme involves two 
atoms located in two separate cavities. We shall consider two 
cases, where there is only a single or two excitations present in 
the system. Thus, in the first case, we will choose the Hilbert 
space of the system spanned by four product state vectors, de- 
fined as follows 

la) = |ei)<^|52)® |0)i® 10)2, 
16) = ® |e2) ® |0)i <E> 10)2, 
16) = I5i)<^l52)'^|l)i'^|0)2, 
1^4) = |5i)^ 152)^ |0)i® 11)2. (40) 

Here, for example, |ei)(8) |(?2) |1)]^ <8) |0)2 represents the state 
in which the atom 1 is in the excited state, the atom 2 is in the 
ground state, one photon is present in the mode of the cavity 1, 
and no photons are present in the mode of the cavity 2. 
In the second case of this scheme, we will assume that there 
are two excitations present, one in each sub-system, and 
choose the Hilbert space of the system spanned by the vec- 
tors 

Ixi) = |ei)® lea)® |0)i® 10)2, 
1X2) = |ei)®|ff2)® |0)i® 11)2, 
1X3) - |5i)<® |e2)(8)|l)i(8)|0)2, 
1X4) = |5i)'»l52)<» 11)1® 11)2- (41) 

In this case we also include the ground state 

Ixo) = \9i) ® I52) ® |0)i ® 10)2, (42) 

for which there is no excitation present in the system. In prac- 
tice, this state could result from a loss of the excitation due 
to spontaneous emission from the atoms or through the damp- 
ing of the cavity modes. The possibility of the presence of 
this state will allow us to study the case where the initial state 
is in the form of a non-maximally entangled Bell state with 
correlated spins 

IXsd) = cosalxi) +e*^sina|xo), (43) 
where a and f3 are real numbers. 

In the fourth scheme, which involves atomic ensembles cou- 
pled to radiation modes of a ring cavity, we discuss entangle- 
ment procedures in terms of the field (bosonic) representation 
of the atomic operators. In this representation, we first express 
each atomic ensemble in terms of collective atomic operators 

N N 
Jtn-T.^tne''--'^", J.n=Y.S^n, (44) 

where N is the number of atoms in the nth ensemble, km is the 
propagation vector of the mth field mode, Vjn is the position 
vector of the 7th atom in the nth ensemble, St„ , S7„ , and S^„ 
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are the raising, lowering and the energy difference operators, 
respectively, of the jth atom in the nth ensemble. 
Next, we reformulate the operators in terms of boson variables 
by applying the Holstein-Primakoff representation of the col- 
lective atomic operators |70J. The representation transforms 
the collective operators into harmonic oscillator annihilation 
and creation operators, c„ and cjj of a single bosonic mode 



fet,„ -4\/^-clc„, J^„ = 4c„ - 7V/2. (45) 



J, 



Because the ensembles are independent of each other, the 
bosonic operators satisfy the commutation relation 



(46) 



where the subscripts n and m enumerate the ensembles. The 
commutation property of the modes will allow us to prepare 
each mode separately in a desired state and independent of 
the other modes. In other words, an arbitrary transformation 
performed on the operators of a given mode will not affect the 
remaining modes. 



A. Two atoms in free space 

The simplest model for the study of entanglement is a sys- 
tem composed of two two-level atoms interacting in free space 
with a quantized multimode electromagnetic field. The field 
appears as a vacuum reservoir (environment) to the atoms 
that causes initial atomic excitation and coherence to decay 
in time. The atoms may decay independently or their radi- 
ation field could exert a strong dynamical influence on one 
another through the vacuum field modes. This would result in 
cooperative behavior of the system. In this section we concen- 
trate on the problem of entanglement creation in the process 
of spontaneous emission from an excited initially unentangled 
state. We shall demonstrate under what kind of conditions the 
initially separated atoms become entangled during the sponta- 
neous emission process. We explore the role of the collective 
states of the system in the entanglement creation via sponta- 
neous emission. Our analysis reveals the necessity of includ- 
ing the antisymmetric state into the dynamics of the system if 
the spontaneous decay of an initial excitation is to generate an 
entanglement during the evolution of the system. 
The first step is to define the approach we take to study the 
spontaneous dynamics of the system. We shall work in the 
density matrix formaUsm and adopt the master equation tech- 
nique to study the time evolution of the density matrix of the 
system. For two two-level atoms interacting in free space with 
a vacuum field at zero temperature, the evolution of the den- 
sity matrix of the system is given by the Lehmberg— Agarwal 
master equation 1I21I I22I1 



dp 

dt 



1 2 

7^T.^'^^([P^t^^i] + [St,S-p]), (47) 



where where S^, S~ , and Sf are the dipole raising, lower- 
ing, and population difference operators, respectively, of the 
ith atom, and = 7 are the spontaneous decay rates of 
the atoms, equal to the Einstein A coefficient for spontaneous 
emission. The terms in the master equation that depend on jij 
and rtij {i 7^ j) are the so-called collective terms because they 
result form the mutual exchange of photons between the atoms 
and thus determine the atomic interaction. The parameter -fij 
represents the collective damping which results from an inco- 
herent exchange of photons between the atom. The collective 
damping leads in general to a change in the lifetime of the col- 
lective states from the single-atom radiative lifetime. The pa- 
rameter riij represents the collective shift of the atomic levels 
and results from a coherent exchange of photons, the dipole- 
dipole interaction between the atoms. The effect of ili2 on the 
atomic system is the shift of the energy of the single excitation 
collective states from the single-atom energy. The collective 
parameters are given by the expressions lll9l - l23l ItiII 



7ij 



>7 



i-ifi-f-ijY 



l-3(/i • f^jY 



cos{krij) sin(fcry) 
(krij) {krij) 



(48) 



and 



cosikrij) 
kvj/, 



l-3(/i • fijY 



sin(fcrij) ^ cos(fcrjj) 
[kn^f {krijf 



(49) 



where ft is the unit vector along the dipole moments of the 
atoms, which we have assumed to be parallel {jj, = pi = pj), 
fij is the unit vector in the direction of Vij , k = ujq/c, and 7-y 
is the distance between the atoms. 

The parameters ( |48] ) and ( l49b are oscillatory functions of krij 
multiplied by inverse powers of fcr^ ranging from {krij ) ^'^ to 
(krij)^^. In the two limiting cases of fcr^ ^ 1 and krij ^ 1, 
the parameters simplify to the following forms 



7ij 



and 



where 



-0{kr,j] 
0{kr,,Y 



0{kr,jY 



(kr.j < 1) 
(krij > 1), 



{kr,j < 1) 
(fcry > 1), 



37 



4{kr,jY 



1 - 3 (/i • njY 



(50) 



(51) 



(52) 



For small krij, the collective damping 7^ becomes equal to 
7, and reduces to the quasistatic dipole-dipole interaction 
potential. For large krij, the collective effects become negli- 
gible and the system reduces to two independent atoms. 
We choose the collective states ( fTSl ). as a convenient basis for 
the representation of the density operator and the analysis of 
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the spontaneous dynamics of the system. Consider first the 
evolution of the diagonal matrix elements, which correspond 
to the populations of the collective states. The equations of 
motion are found from the master equation ( l47l l. and are given 
by 

PiA = -27P44, 

Pss = - (7 + 712) {Pss - P44) , 

Paa = - (7 - 712) {Paa " P44) • (53) 



two-photon state |^'4). The density matrix elements evolve as 

Pas = - (7 + 2ir2l2) Pas, 

Pi4 = - (7 - 2iwo) Pi4- (54) 

The other set contains four coupled equations for the one- 
photon coherence of the cascade transitions from the upper 
two-photon state l^i^) through the intermediate states \'^s) 
and l^'a) to the ground state l^*!). The density matrix ele- 
ments evolve as 



Equations ( l53T l are in the form of rate equations for the popula- 
tions of the collective states. The equations give us an impor- 
tant information about the transitions rates between the col- 
lective states 1 721 47411 . One can see that the transitions rates to 
and from the symmetric and antisymmetric states are changed 
by the collective damping 712. The transitions to and from the 
symmetric state occur with an enhanced rate 7 + 712, whereas 
the transitions to and from the antisymmetric state occur with 
a reduced rate 7 — 712. The collective damping depends on 
the distance between the atoms. For small kri2, the state j^I's) 
becomes supermdiant with a decay rate double that of the sin- 
gle atom 7, and the state \ '^a) becomes submdiant, with a de- 
cay rate of order (fcri2)7 which vanishes in the limit of small 
distances fcri2 <C 1. For large fcri2, the effects of the collec- 
tive damping vanish as (fcri2)^^ and become negligible for 
fcri2 ^ 1. In this limit, the states decay with a single rate 
identical to that of a single atom. 



|W4> 



|Ws> 




|Wa> 



|Wl> 



FIG. 1 : The collective states of a two-atom system with possible one 
(solid arrows) and two-photon (dashed arrows) transitions. Transi- 
tions — >■ \^ s) — > I'I'i) through the symmetric state occur with 
an enhanced rate 7+712, while the transitions 1^4) \^a) ^ Y^i) 
through the antisymmetric state occur with a reduced rate 7 — 712. 
The intermediate states are shifted from their resonant positions by 
Q,\2, the energy of the dipole-dipole interaction. 



Using Eq. ( |47] |. we find that the off-diagonal matrix elements 
are decoupled from the diagonal elements and obey equations 
of motion that can be grouped into two independent sets: One 
containing two independent equations for the coherence pas 
between the symmetric and antisymmetric states, and the two- 
photon coherence pi4 between the ground state Y^i) and the 



Ps4 = — riPs4, Pai = —^2Pa4, 

Pis = -T3PIS + (7 + 712) Ps4, 

Pla = -r4Pla - (7 - 712) Pa4, (55) 

where, for simplicity, we have introduced the notation 

ri^2 = ^ (37 ± 712) -i{ojoT ^12) , 

r3,4 = i (7± 712) - i (wo ± ^^12) • (56) 

The one-photon coherences ( fSST i oscillate with frequencies 
shifted from the unperturbed frequency ljq by the dipole- 
dipole interaction energy r2i2. Similar to the decay of the 
populations, the coherence involving the symmetric state de- 
cay with an enhanced rate whereas the coherence involving 
the antisymmetric state decays with a reduced rate. 
The collective states of the system are shown in Fig.[T] where 
it is seen that the intermediate states are shifted from their 
unperturbed energies by the dipole-dipole interaction energy, 
and there are two transition channels [^'4) — > j'l's) — ^ 
and |4'4) — |^'a) — > l^*!), each with two cascade nondegen- 
erate transitions. The transition channels are uncorrected, but 
the transitions inside these channels are damped with different 
rates. 

Our objective is to calculate the time evolution of the con- 
currence under the spontaneous decay of an initial excitation. 
To this end, we need the time evolution of the density matrix 
elements, which is obtained by solving the equations of mo- 
tion (l53])-(l55b. The general solution valid for arbitrary initial 
conditions is given by 

P44(t) = P44(0)e-2^*, 



7 + 712 



Paa{t) 



P44(0) 

7 - 712 



^(7-712)* 



+ /O44(0) 



7 - 712 



7 + 712 



3(7+712)* 



1 



-27* 



,-27* 



(57) 



andpii(t) = 1 - P44(t) - pss{t) ~ Paa{t)- 

We see from Eq. ( ISTl i that the decay of the populations de- 
pends strongly on the initial state of the system. When the 
system is initially prepared in the state s), the population 
of the initial state decays exponentially with an enhanced rate 
7 + 712, while the initial population of the antisymmetric state 
decays with a reduced rate 7 — 712. This occurs because the 
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photons emitted from the excited atom can be absorbed by the 
atom in the ground state, so that the photons do not escape 
immediately from the system. For a general initial state that 
includes in the state |^'4), the populations of the symmetric 
and the antisymmetric states do not decay with a single expo- 
nential. 

Similarly, Eqs. ( l54l i and i55[ are straightforward to solve, giv- 
ing 

Pi4{t) = pi4(0)e-(^-2*'-»)*, 
pUt) ^ pu{0)e-^'' 

Plait) = Pla{0)e'^'' 

7 - 712 



-Pa4(0) 



-r,* 



7 - 2ifli2 



(58) 



The density matrix elements depend on the initial state |/) 
of the system. Since we are interested in creation of entan- 
glement from an initial separable state, we take for the initial 
state of our system the single-excitation state |/) = l'^^) = 
|ei) (8) I52), which corresponds to atom 1 in the excited state 
and atom 2 in the ground state. The initial state is, of course, 
a separable state, in other words there is no entanglement in 
the system at t = 0. For |/) — | 'I' 3), it is easily verified that 
the only non-vanishing matrix elements are 



PssiO) = Paa(O) = Psa(O) 



,(0)4. 



(59) 



and then the initial density matrix has the following form 



P(0) 



/O 

PssiO) PsaiO) 

PasiO) PaaiO) 

\ 



(60) 



According to the solutions for the density matrix elements, 
Eqs. $5% and ( l58b , the matrix elements which are zero at the 
initial time t = will remain zero for all time, except the pop- 
ulation of the ground state which will buildup during the 
evolution. Hence, the diagonal form of the density matrix will 
be preserved during the evolution. Moreover, the dynamics of 
the systems can be confined to the subspace spanned by three 
state vectors only, l^^i), jvPs) and \^a)- 
Our next and, in fact, the major problem is to determine if 
an entanglement can be generated in the system ifTillTSl - lTSll . 
To examine the occurrence of entanglement, we must con- 
sider the concurrence, Eq. ( |29] l which, on the other side, is 
described by the density matrix elements. Looking at the den- 
sity matrix (|60] |, two conclusions can be made. First of all, 
since the coherence pi4(t) is equal to zero for all times t, we 
see from Eq. (|30] | that the criterion Ci (t) is always negative. 
Consequently, Ci{t) will not contribute to the concurrence. 
Secondly, the two-photon state {'^4) is not involved in the dy- 
namics of the system which, according to Eq. dSTT i rules out 



a possibility for the threshold behavior of the criterion C2{t). 
Under this situation, the concurrence is determined only by 
the criterion C2{t) which would always be positive. From this, 
it follows that 



C{t) EE C2{t) - \/[pss{t) - Paa{t)f + ^[Impsaf 

= e-'''Jsmh^{ji2t)+sm\2ni2t). (61) 



There are two contributions to the time evolution of the con- 
currence. First, there is an imaginary part of the coherence 
between the states \'^s) and |^'a). Since this is off-resonance 
coupling, it leads to oscillation in the concurrence with fre- 
quency 2J7i2, the frequency difference between the two states. 
The second contribution is the difference between the popu- 
lations of the states j'l's) and |^'a). This contribution is non- 
oscillatory. We see that the effect of the non-oscillatory term 
is apparently to lengthen the lifetime of the concurrence. 




FIG. 2: Transient behavior of the concurrence (solid line) and the 
populations of the antisymmetric (dashed line) and the symmetric 
(dashed-dotted line) states for the atoms prepared initially in the sep- 
arable state I'I's) — |ei) ® \g2), with p. 1. fi2, and fcri2 — 7r/5. 

The features described above are easily seen in Fig. |2] where 
we plot the time evolution of the concurrence C(i), calculated 
from Eq. dMT l. together with the populations pss (t) and pau (t) 
of the symmetric and the antisymmetric states, respectively. 
The time evolution of the concurrence reflects the time evo- 
lution of the entanglement between the atoms. It is seen that 
the concurrence builds up immediately after t — and re- 
mains positive for all time, which indicates that spontaneous 
emission can indeed create entanglement between the initially 
unentangled atoms. The buildup of the concurrence in time 
generally consists of oscillatory and non-oscillatory compo- 
nents, so that two time scales of completely different behavior 
of the concurrence can be distinguished. At early times, the 
concurrence builds up in an oscillatory manner, and the oscil- 
latory structure is smoothed out on a time scale (7 + 712)"^, 
the lifetime of the symmetric state. The oscillations vanish 
at time close to the point where the symmetric state becomes 
depopulated. At later times, the concuiTence evolves in a non- 
oscillatory manner and overlaps with the population of the an- 
tisymmetric state. As a result, the concurrence decays slowly 
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in time with the reduced rate 7 — 712. Although the concur- 
rence dMT l involves both symmetric and antisymmetric states, 
it is clear that crucial for the entanglement is the presence of 
the antisymmetric state. 

The decay time of the population of the antisymmetric state, 
so that the transient entanglement seen in Fig. |2] varies with 
the distance between the atoms. The time goes to infinity 
when fcri2 — ^ 0. In this limit the transition rates to and from 
the antisymmetric state vanish and the state decouples from 
the remaining states. Hence, any initial population encoded 
into the state will remain there for all times. For example, 
with the initial state j'l's) and t — 7> 00, half of the popula- 
tion, that initially encoded into the antisymmetric state, still 
remains in the atomic system, and half of the population, that 
initially encoded into the symmetric state, is emitted into the 
field. As a result, the concurrence evolves to its stationary 
value of C{t) = 0.5 indicating that in the limit of t — ^ 00, 
the spontaneous emission could produce steady state entan- 
glement to the degree of 50% with the corresponding pure 
state of the system. 

The creation of the transient entanglement can be understood 
by considering the properties of the density matrix of the sys- 
tem. For the initial state of only one atom excited the den- 
sity matrix is not diagonal due to the presence of coherences 
/Osa(O) and pas(O). Since the form of the matrix is preserved 
during the evolution, it remains non-diagonal for all times. 
Consequently, the collective states are no longer the eigen- 
states of the system. The density matrix can be diagonalized 
to give new diagonal states. It is easy to verify that the states 
l^'i) and [^'4) remain unchanged, whereas the states \'^s) and 
1 4* a) recombine into new diagonal symmetric 



[p++{t)-pssm'fa)+Pas{t)\'i>s 



(62) 



[[p++{t)-pUt)f + \Pas{t)\^}" 

and antisymmetric 

|^_^ ^ Pasima) + [P—it)-Paaitms) ^^3^ 
{[p^^it)-paa{t)f + \pas{t)\'}^ ' 

States, with the eigenvalues p++{t) and p (t), the popula- 
tions of the diagonal states, given by 

= ^ [Paa{t) + Pss{t)] 



^ 1 

-{[Paait)-pssit)f+MPas{t)\''y , 



P—it) = ^[Paa{t)+Pssit)] 

I 1 

"2 {[Paa{t)-pss{t)f+MPa.{t)\'y. (64) 

It follows from Eq. ( l63l l that the coherences psa {t) and pas {t) 
cause the system to evolve between the ground state and 
two "new" one-photon states and l^*-), which are lin- 
ear combinations of the collective states \'$s) and |5'a)- It is 
easy verified that for the initial condition ( |59] l, the population 



p (t) = for all times, whereas the population p++ (t) is 

different from zero and equals to the sum of the populations 
Pss{t) and Paa{t). The lack of population in the states l^*-) 
together with no population in the state reduces the four- 
level system to an effective two-level system with the excited 
nonmaximally entangled state |^'+) and the separable ground 
state 1^1). In this case, the density matrix has a simple diag- 
onal form 



= p++(<)|*+)(*+|+pii(t)|*i)(*i 



(65) 



Since p [t) ~ for all times t, we can easily find 

from Eq. (|64] | that in this case \pas{t)\^ — Paa{t)Pss{t), and 
then the state |^+) can be written as 



1*4 



^Pss{t)\^s) + ^Paa{t)\^a) 
y/Paa{t) + Pss{t) 



(66) 



When the collective states \"i's) and |*a) are equally popu- 
lated, the state reduces to a separable state 



= lei 



1.92) 



(67) 



On the other hand, the state |^'+) reduces to a maximally en- 
tangled state, I^Ps) or |^'a), when either pss{t) or Paa{t) is 
equal to zero. Since the population of the symmetric state de- 
cays faster than the antisymmetric state, see Eq. (|53] |. at time 
when the state \'^s) becomes depopulated, the state |^'+) re- 
duces to the maximally entangled antisymmetric state |^'a)- 
This explains why at later times the evolution of the concur- 
rence follows the evolution of the population of the state | *a) ■ 
The above analysis give clear evidence that the creation of 
transient entanglement from the separable state by sponta- 
neous emission depends crucially on the presence of the anti- 
symmetric state. 



B. The Dicke model 

The most familiar model to study collective effects in sponta- 
neous emission by a system of two or more identical atoms is 
the Dicke model 1 18]. It was originally introduced by Dicke 
in his famous article published in 1954, and currently a very 
large literature exists on a wide variety of problems involving 
the model, in particular on the concepts of super-radiance and 
the directional propagation of light in atomic ensembles ll79l - 
IsTfl . Recently, the model has been employed in the studies 
of entanglement, in particular for creation of entanglement in 
atomic ensembles cornposed of a large number of trapped and 
cooled atoms 1*411 14114511 . 

The two-atom Dicke model is a simplified form of the two- 
atom system considered in the preceding section. It is often 
called, a small sample model and corresponds to a system of 
two atoms confined to a region much smaller than the radia- 
tion wavelength of the atomic transitions. In other words, the 
model assumes that the atoms are close enough that we can 
ignore any effects resulting from different spatial positions of 
the atoms. Mathematically, it is equivalent to set the phase 
factors exp(ifc • 7\) associated with atomic positions, equal to 
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one. In terms of the collective states of a two-atom system, 
the Dicke model involves only the triplet states |'I'4), \'^s), 
and l^*!), with the antisymmetric state I'J'q) totally decoupled 
from the triplet states and thus not participating in the evolu- 
tion of the two-atom system. The antisymmetric state appears 
as a metastable or trapping state. In the terminology of mod- 
ern quantum optics, the state is an example of a decoherence- 
free state that any initial population will remain in the state for 
all times. 

The trapping property of the antisymmetric state could have a 
significant effect on entanglement in the small sample model 
of two atoms. According to Eq. ( |39] |, the concurrence depends 
on the population of the antisymmetric state. Thus, an initially 
entangled or separable system would evolve to a steady state 
whose the entangled properties are determined by the initial 
population of the antisymmetric state. In this way, one could 
create a stable entanglement between the atoms simply by a 
suitable preparation of the atoms in the non-radiating antisym- 
metric state. For example, if the system is initially prepared 
in the antisymmetric state, €2(0) ~ Paa{Q) = 1, and then 
due to the trapping property of the antisymmetric state that 
the population in the state does not change in time, the system 
could stay maximally entangled for all < > and would never 
disentangle. 

If the system is initially prepared in a state with no population 
in the antisymmetric state, Paa{0) = Paa{t) = 0, and then 
the concurrence is determined only by the populations of the 
triplet states 



condition, it follows from Eqs. (I68I I and (|69] l that 



C{t) = max{0, (r - 2Ve^ - 1 - r) e"^} , 



(70) 



Cit) = max {0, pss{t) - 2^pnit)pMit)} , 



(68) 



Hence, one would expect that the absence of the population 
in the antisymmetric state could result in a large transient en- 
tanglement buildup through the population of the symmetric 
state. 

We demonstrate a somehow surprising result that in the Dicke 
model, entanglement cannot be created by spontaneous emis- 
sion if one disregards the antisymmetric state. In the Dicke 
model, the time evolution of the diagonal density matrix el- 
ements under the spontaneous emission is determined by the 
following equations 1 14] 

P44(i) =/344(0)e-2^*, 

P. At) = PssiO) e-^-" + 27tp44(0) 

Paait) = PaaiO). (69) 

The population of the symmetric state decays with a rate 
double that of the single atom, and for an initial state with 
/544(0) 7^ 0, the population does not decay exponentially. The 
time evolution of the population is a convolution of a linear 
function of time and an exponential decay. 
Let us assume that the atoms are initially prepare in the sep- 
arable state 1^4), which implies that ^44(0) = 1 and all 
other density matrix elements equal to zero. Since in the 
Dicke model, the population of the upper state can only decay 
through the symmetric state, the population will accumulate 
in this state during the evolution. As a result, an entanglement 
may emerge in the system. However, for the ^44(0) — 1 initial 



where r = 2'jt. 

The exponent under the square root in Eq. ( iTOl i may be devel- 
oped into series and one then obtains the following expression 
for the concurrence 



C{t) ^ max<^ 0, re""" 



(71) 



Since the term in the square brackets is always negative, this 
shows that no entanglement can be created during the sponta- 
neous decay of the initial excitation. This is a surprising result, 
because the symmetric state can be significantly populated 
during the spontaneous decay and still no entanglement can 
be created in the system. It is easy to find from Eq. ( l69b that 
at early times the population pss {t) increases linearly from its 
initial value of Pss{0) ~ 0, attains the maximum of 0.42 at a 
time t = 1/ (27), and then for t > 1/ (27) decreases exponen- 
tially with the rate 27. 

The failure of the spontaneous creation of entanglement in the 
Dicke model is linked to the threshol d behavior of the concur- 
rence. Simply, the threshold term 2y/pijjt)p44jt), which de- 
pends on the redistribution of the population between the two 
separable states is dominant in the concurrence for all time 
despite the fact that the symmetric state can be significantly 
populated during the evolution. In the coming sections, we 
shall explore in more details the role of the threshold in the 
evolution of initial entangled states and a delayed creation of 
entanglement from an initial separable state. In the terminol- 
ogy of sudden death of entanglement that will be introduced 
in the next section, we may conclude that the Dicke model is 
"dead" for creation of entanglement by spontaneous emission. 



VI. STRANGE BEHAVIORS OF ENTANGLEMENT 

In the preceding section we have developed a simple model 
for creation of entanglement in a system composed of two 
two-level atoms interacting with a common environment. We 
have showed how the spontaneous decay of an initial excita- 
tion encoded into a separable state of the system can create a 
transient entanglement between the atoms. We will now ex- 
amine the opposite situation, a spontaneous disentanglement 
of initially entangled atoms. Typically, an initial entangle- 
ment is expected to decays exponentially in time. However, 
we point out some "strange" behavior of the disentanglement 
process of a two-atom system undergoing spontaneous evolu- 
tion that an initial entanglement encoded into the system can 
be lost in a very different way compared to an exponential de- 
cay. We provide a discussion of this unusual phenomenon in 
terms of the density matrix elements and show the connection 
of the phenomenon with the threshold behavior of the concur- 
rence. We also discuss other unusual features of entanglement 
and disentanglement. Namely, we explore an interesting phe- 
nomenon of entanglement revival that spontaneous emission 
can lead to a revival of the entanglement that has already been 
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destroyed. In the connection to the threshold behavior of the 
concurrence, we will also return to the problem of the dynam- 
ical creation of entanglement from an initial separable state. 
We discuss a phenomenon of delayed sudden birth of entan- 
glement that the spontaneous creation of entanglement may 
be postponed to later times even if the correlation between the 
atoms exists for all time. 



A. Sudden death of entanglement 

The most familiar of the "sudden" features of entanglement 
and disentanglement is the phenomenon of entanglement sud- 
den death, i.e., abrupt disappearance of the entanglement at 
a finite time even if the correlation between the atoms ex- 
ists for all time. The subject received its initial stimulus in 
an article by Yu and Eberly iIztI l28ll . in which they for the 
first time introduced the concept of sudden death of entan- 
glement. Many authors since have dealt with the entangle- 
ment sudden death in systems composed of two atoms or two 
harmonic oscillators. Most of the subsequent literature can 
be divided into two categories: (i) studies which deal with 
independent atoms interacting with local environments lf82ll . 
and (ii) studies which deal with interacting atoms coupled to 
a common environment [83]. We postpone the studies of the 
latter group to the following section, where we will mostly 
focus on the phenomenon of entanglement revival. Here, we 
focus on the original concept of Yu and Eberly, and discuss 
in details the phenomenon of sudden death of entanglement 
in a system of two independent atoms interacting with local 
environments. This model is also applicable to a situation of 
two distant atoms interacting with a common environment. At 
large distances, the collective parameters 712 and fii2 are very 
small, so that the interaction between the atoms can be ignored 
and the atoms can be treated as independent sub-systems. 
Suppose that at t = the system of two independent atoms is 
prepared in a non-maximally entangled state of the form 



|To> = V^lei) ® |e2) + Vl^l.9i) I52), 



(72) 



where g is a positive real number such that < q < I- The 
state corresponds to an excitation of the system into a coherent 
superposition of its product states in which both or neither 
of the atoms is excited. In the special case of q = 1/2, the 
state ( |72] ) reduces to the maximally entangled Bell state |^'s). 
Let us next consider the time evolution of the concurrence 
when the system is initially prepared in the state (|72] ). It is 
not difficult to verify that the initial values for the density ma- 
trix elements are 



P44(0) = q, pi4(0) - ^9(1- Pll(O) = 1 



(73) 



and the other matrix elements, the populations of the symmet- 
ric and antisymmetric states, and all one-photon coherences 
are zero, i.e. pss{0) = Paa(O) = and pes{0) = Pea{0) = 
Psg(O) = Pag(O) = PasiO) 0. According to Eq. 
the coherences will remain zero for all time, that they can- 
not be produced by spontaneous decay. However, the pop- 
ulations pss{t) and Paa{t) can buildup during the evolution. 



This implies that for all times, the density matrix of the sys- 
tem spanned in the basis of the collective states ( fTSl l. is in the 
X-state form 



Pit) 



( Pii{t) 

Pss{t) 



\pii{t) 




(74) 



with the density matrix elements evolving as 



P44(i) 

Pli{t) 
Pss{t) 



qe 



'2'yt 



Paa{t) =q{l 



(75) 



subject to conservation of the trace of p{t): pii{t) = 1 — 
Pss (t) — paa (t) — P44 (t). It is to be noticcd that the symmetric 
and antisymmetric states are equally populated for all time. 
This holds for any initial state and results from the fact that the 
atoms radiate independently from each other. In this case, the 
populations of the collective states decay with the same rate, 
equal to the single-atom damping rate. Therefore, the initial 
relation between the populations cannot be changed during the 
spontaneous emission. 

The density matrix (|74] | leads to a particularly simple expres- 
sion for the concurence. In general, the concurrence is given 
in terms of two entanglement criteria Ci (t) and C2 {t), as seen 
from Eq. ( |29] |. However, with the initial state (l72b . the two- 
photon coherence puit) is different from zero and the sym- 
metric and antisymmetric states are equally populated for all 
time. As a consequence, the criterion C2 (t) is always negative, 
irrespective of q and times t. Therefore, entangled properties 
of the system are solely determined by the criterion Ci (t). On 
substituting from Eq. ( fTSl ) into Eq. ( l30b , we obtain the follow- 
ing expression for the concurrence 



at) 



{0, Ait) e"'^*} 



where 



Ait) = 2^9(1-9) 



1 - 



1-q 



-ft 



(76) 



(77) 



This shows that the major features of the entanglement are de- 
termined by the properties of A(<) which, on the other hand, is 
dependent on the parameter q. We see that there is a threshold 
for values ofq;q—l /2, below which A(t) is always positive. 
However, above the threshold, A(t) can take negative values 
indicating that the initial entanglement can vanish at a finite 
time. Consequently, the sudden death of the entanglement is 
possible for initial states with q > 1/2. Since (044(0) — q, we 
can conclude that the entanglement sudden death is ruled out 
for the initially not inverted system. 

Figure |3] shows the concurrence C(t), calculated from 
Eq. (TTSl l. as a function of time for two different values of the 
parameter q. It is evident from the figure that for q < 1/2 
the initial entanglement decays exponentially in time without 
any discontinuity. The entanglement sudden death appears for 
q > 1/2 that the concurrence decays in a non-exponential 
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way and vanishes at a finite time. In addition, we plot the 
two-photon coherence |pi4(i)| for g — 2/3. It is apparent that 
the coherence decays exponentially in time, which clearly il- 
lustrates that the entanglement disappear at finite time despite 
the fact that the two-photon coherence is different from zero 
for all time. 

As we have already stated, time at which the entanglement 
disappear is a sensitive function of the initial conditions deter- 
mined by the parameter q. It is easily verified from Eq. ( f76] l 
that the time id at which the entanglement disappears is given 
by 



1 \ 2g-l 



9) 



(78) 



The time td gives the collapse time of the entanglement be- 
yond which the entanglement disappears. The dead zone of 
the entanglement continues till infinity that the entanglement 
never revive. It may continue for a finite rather than infinite 
time that under some circumstances the already dead entan- 
glement may revive after some finite time. A revival of the 
entanglement may occur when the atoms directly interact with 
each other. We leave the discussion of this problem to the fol- 
lowing section. 



conserving, phase matched modes. Since the pairs of po- 
larized photons are spatially indistinguishable, they are de- 
scribed by a pure state 



\^) = \o.\\lili) 



(79) 



where the coefficients |a| and and the phase b were ad- 
justed by applying half- and quarter-wave plates in the laser 
beam pumping the crystals to control the creation of pairs of 
a desired polarization. 



A for p = 3 a, experiment 
C for p = a/3 experiment 

A for p = 3 a, tlieory 

C for p = a/3 theory 

C for p = 3 a, theory 





FIG. 3: Time evolution of the concurrence C(i) for two values of the 
parameter q: g = 1/3 (solid line) and g = 2/3 (dashed line). We 
also plot the time evolution of the two-photon correlation |pi4(t)j 
(dashed-dotted line), with g = 2/3. 

Here, we would like to point out that first demonstration 
of entanglement sudden death has recently been reported by 
Almeida et al. IsTll . The apparatus used in the experiment 
involved a tomographic reconstruction of the density matrix 
and from it the concurrence by measuring polarization en- 
tangled photon pairs produced in the process of spontaneous 
parametric down-conversion by a system composed of two ad- 
jacent nonlinear crystals. One of the crystals produced pho- 
ton pairs with V^-polarization and the other produced pairs 
with i7-polarization. Parametric down conversion is a non- 
linear process used to produce polarization entangled photon 
pairs, which are manifested by the simultaneous or nearly 
simultaneous production of pairs of photons in momentum- 



FIG. 4: Results of the Almeida et al. experiment demonstrating the 
sudden death of entanglement. The theoretical prediction for the con- 
currence C is plotted as a function of p = 1 — exp(— 7t) for two val- 
ues of the ratio |/?|/q|: |/?|/q| = 3 (solid line) and |/3|/|q| = 0.333 
(dashed line). The squares and triangles are experimentally mea- 
sured values for the concurrence. Reproduced from Ref. Bill with 
permission. 

In the experiment, they measured the decay of a single polar- 
ized beam, serving as a qubit, which was monitored by gener- 
ating pairs of photons of the same polarization and registering 
coincidence counts with one photon propagating through the 
interferometer and the other serving as a trigger. 
Figure |4] shows the results of the measured concurrence for 
two different values of the ratio |/3|/|a|. The solid and dashed 
lines represent the theoretically predicted concurrence. The 
measured values for the concurrence are found to be in good 
agreement with the theoretical predictions. Thus, it was con- 
firmed that entanglement may display the sudden death fea- 
ture and that the spontaneous decay of the initial entanglement 
depends on the relation between the coefficients \a and \(3\. 
For 1/3 1 < |a| the entanglement decays exponentially in time, 
while for > \a\, the entanglement vanishes at finite time. 
In closing this section, we summarize the research on the en- 
tanglement sudden death. The phenomenon has been exten- 
sively studied in recent years and it has been demonstrated that 
the entanglement sudden death can be obtained in a variety 
of systems including atoms coupled to single-mode cavities, 
atoms coupled to local environments appearing as multi-mode 
reservoirs to the atoms, atoms interacting with a common en- 
vironment. Other studies have been carried out for indepen- 
dent and also for interacting harmonic oscillators. Studies 
have also been carried out for non-Markovian and non-RWA 
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situations where the sudden death of entanglement can be ob- 
served in completely different regime of the parameters. 
Finally, we briefly comment on the evolution of an entan- 
glement initially encoded into spin-anticorrelated states. We 
have just shown that the phenomenon of entanglement sud- 
den death is characteristic of two-photon or spin correlated 
entangled states. It is easy to conclude from Eq. (|39] l that an 
initial entanglement encoded in a spin anti-correlared state, 
the symmetric or antisymmetric states, would decay asymp- 
totically in time without any discontinuity. Expression (3% 
clearly shows that necessary condition for entanglement sud- 
den death is a simultaneous population of the ground j'l'i) and 
the two-photon [^'4) states. 



B. Revival of entanglement 

We have already shown that under appropriate conditions, two 
initially entangled and afterwards not interacting qubits can 
become completely disentangled in a finite time. The term 
"entanglement sudden death" has been introduced to describe 
this unusual feature. In this section, we extend these analysis 
to include the direct interactions between the atoms and study 
in details the time evolution of an initial entanglement. As 
we shall see, the inter-atomic interactions can lead to transient 
features of a spontaneously evolving entanglement which are 
quite distinct to those seen in their absence. The interactions 
impose strong qualitative changes to the time evolution of en- 
tanglement that the irreversible spontaneous decay can lead 
to a revival of the entanglement that has already been de- 
stroyed 1 30]. We pay particular attention to the role of the 
collective damping 712 and the dipole-dipole interaction ili2 
in entanglement revival and find that the feature depends crit- 
ically on whether or not the collective damping is present. We 
study in detail the dependence of the revival time on the initial 
state of the system and on the separation between the atoms. 
Consider again the initial state of the system given by Eq. ( l72b . 
but now assume that the atoms can interact with each other. In 
this case, the initial values of the density matrix of the system 
as the same as for independent atoms, but the time evolution 
of the density matrix elements is different. It is now given by 



Pss{t) = q 

Paa{t) = q 



7 + 712 
7 - 712 
7 - 712 
7 + 712 



3(7-712)* 



3(7+712)* 



-27* 



-27* 



(80) 



and pii(t) = 1 - pn{t) - pss{t) - Paa{t). As before for 
independent atoms, the remaining density matrix elements 
are equal to zero. By comparing Eq. dSOl l with Eq. ( iTST i for 
independent atoms, we see that a major difference between 
the time evolution of the density matrix elements for inde- 
pendent and interacting atoms is that in the later the sym- 
metric and antisymmetric states are populated with different 
rates. As a result, the transient buildup of the populations 



from Pss(O) = Pqq(O) = at t = will lead to unequal 
populations of the states for all later times t > Q. 
A direct consequence of unequal populations of the symmet- 
ric and antisymmetric states is seen in the transient evolu- 
tion of an initial entanglement which may now depend on 
both criteria Ci{t) and C2(i), that are needed to construct the 
concurrence C{t), rather than on one of them. According to 
Eqs. ( |29] l and (|73] |. a two-atom system initially prepared in the 
state ( f72] i can be entangled according to the criterion Ci{t), 
and the degree to which the system is initally entangled is 
Ci(0) = 2^J q{l — q). The other criterion, C2(t), is negative at 
t = 0. Nevertheless, during the spontaneous evolution, a pop- 
ulation builds up in the symmetric and antisymmetric states 
and according to Eq. dSOl l. pss (t) 7^ Paa (t) for all t > 0. As 
discussed in Sec. 4, an unequal population of the states gives 
rise to positive values of the criterion C2 (t) which then may 
result in an entanglement of the atoms. Of course, in order to 
create the entanglement, the the population difference has to 
overweight the threshold term in the criterion C2{t). If this is 
the case, we could see an entanglement having its origin in the 
correlation determined by the criterion C2{t). 
According to the above discussion, the concurrence, in gen- 
eral, is determined by both, Ci (t) and C2 (t) criteria 



where 



C{t) = max{0, Ci{t), C2{t)} 



Ci{t) ^ 2\pi4t)\ - [pss{t) + paa{t)] , 



and 



C2{t) = \pss{t) - Paa{t)\ - 2 ^Pll (i)p44 W • 



(81) 



(82) 



(83) 



We point out again, that initially at t = 0, C2 (0) is negative 
and then for t > it may rise to positive values. Thus, an 
entanglement can be generated during the spontaneous evolu- 
tion of the system. This kind of entanglement is an example 
of spontaneously generated entanglement. 
Let us illustrate the time evolution of the concurrence for a 
collective system with 712 7^ 0. Figure|5]shows the deviation 
of the time evolution of the concurrence for two interacting 
atoms (712 7^ 0) from that of independent atoms (712 = 0). 
In both cases, the initial entanglement falls sharply in time. 
For independent atoms we observe the collapse of the entan- 
glement without any revivals. However, for interacting atoms, 
the system collapses over a short time and remains disentan- 
gled until a time tr w 1. 5/7 at which, somewhat counter- 
intuitively, the entanglement revives. This revival then decays 
asymptotically to zero, but after some period of time it revivals 
again. Thus, important features of the transient evolution of 
the entanglement is that there are two time intervals at which 
the entanglement vanishes and two time intervals at which the 
entanglement revives. 

The origin of the sudden death and revivals of the entan- 
glement can be explained in terms of the populations of the 
collective states and the rates with which the populations 
and the two-photon coherence decay. It is easily verified 
from Eq. dSOb that at early times of 7t <C 1, the population 
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Paa{t) « 0, but pss{t) is large. Moreover, the two-photon 
coherence puit) is also large at that short time. Thus, the 
entanglement behavior can be understood entirely in terms of 
the population pss{t) and the coherence pi4{t). In this con- 
nection, we compare in Fig. |6] the transient behavior of the 
concurrence C{t), with the transient properties of the popula- 
tion Pss{t), and the coherence puit) for the same choice of 
parameters as in Fig. |5] As can be seen from the graphs, the 
entanglement vanishes at the time where pss {t) is maximal, 
and remains zero until the time U at which pss (t) becomes 
smaller than pi4{t). Clearly, the first dead zone of the en- 
tanglement arises due to the significant accumulation of the 
population in the symmetric state. 

The reason for the first revival of the entanglement is in the 
unequal decay rates of the population pss{t) and the coher- 
ence pi4{t). According to the criterion Ci{t), entanglement 
would be generated if the sum of the populations of the single- 
photon collective states is small while at the same time the 
two-photon coherence is large. Since the coherence puit) de- 
cays more slowly than the population of the symmetric state, 
and keeping in mind that the population of the antisymmetric 
state is negligibly small at that time, the two-photon coherence 
can become significant and entanglement generated over some 
period of time during the decay. Note, that this is the same co- 
herence that produced the initial entanglement. Therefore, we 
may call the first revival as an "echo" of the initial entangle- 
ment that has been unmasked by destroying the population of 
the symmetric state. 

The revival of the entanglements depends on the initial state 
of the system, which is determined by the parameter q. There 
is a range of values of q for which the concurrence is positive. 
This range defines the values of q for which the first revival oc- 
curs. The range depends on the values of 712 and the revival 
is most pronounced for q over the range 0.85 < q < 0.95. For 
values outside this range, no revival is observed. Note that the 
values of q at which a pronounced revival exists are close to 
one. This is not surprised because for q > 1/2 the system is 



initially inverted that increases the rate of spontaneous emis- 
sion. 

We may estimate approximate values for the first death and 
revival times of the entanglement. Using Eqs. (l30l ) and ( l80l ) 
we find that in the case of small inter-atomic separations, at 
which 712 « 7, the entanglement criterion Ci (<) attains zeroth 
at times satisfying the relation 

-ft exp(-7i) = y^^, (84) 

which for q > 0.85 has two nondegenerate solutions, td and 
tr > td. The solution td gives the death time of the entan- 
glement beyond which it disappears. At this time, the atoms 
becomes disentangled and remain separable until the time tr, 
at which the entanglement revives. The analytic expression 
( |84| ) makes clear that for the parameters of Fig. |5] the sys- 
tem disentangles at td — O.6/7 and entangles again at the 
time tr = 1.7/7. 

Figure |5] shows that entanglement revivals not once but twice, 
and the second revival occurs at long times. The long time en- 
tanglement has completely different origin than that at short 
times. At long times, all the density matrix elements are al- 
most zero except Paa{t), which remains large even for long 
times due to the small decay rate of the antisymmetric state. 
These considerations imply that the long time entanglement 
comes from the population paa {t) which determines positive 
values of the criterion C2 {t). The entanglement decays asymp- 
totically with the reduced rate 7 — 712, and vanishes at 

td2 ~ -arcsinh ( a - — - — — — ) . (85) 

7 VV 9 7-712/ 

This explicitly shows that lifetime of the entanglement de- 
pends on the collective damping parameter 712 and ap- 
proaches infinity when 712 — > 7. 

We may summarize this section that the most interesting con- 
sequence of the collective damping is the possibility of entan- 
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FIG. 5: Transient evolution of the concurrence C{t) for the initial 
state ( 172b with q = 0.9. The solid line represents C{t) for the collec- 
tive system with the interatomic separation ri2 = A/20. The dashed 
line shows C{t) for independent atoms, 712 = 0. 



FIG. 6: The transient behavior of the entanglement criterion Ci{t) 
(solid line), the two-photon coherence 2|pi4(f)| (dashed line) and 
the population p3s(t) (dashed-dotted line) plotted for the same pa- 
rameters as in Fig.|5] 
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glement revival. What is even more remarkable, the revival of 
entanglement occurs under the spontaneous evolution of the 
system without the presence of any external coherent fields. 
We may say that entanglement is reversible even if it evolves 
under the irreversible process. 

Finally, we would like to point out that entanglement revival 
is associated not only with the collective damping, but may 
occur in a number of other situations. For example, it has 
been found that entanglement revival may occur under non- 
Markovian dynamics of two independent atoms coupled to 
local reservoirs |84]. Here, the memory effects of the non- 
Markovian dynamics may result in the "return" of the cor- 
relation from the reservoirs to the atoms. Entanglement re- 
vival has also been found in reversible systems whose dy- 
namics are determined by the Jaynes-Cummings Hamilto- 
nian ll49l [ssl [86ll . In these models, an initial entanglement 
encoded into atoms undergoes a coherent transfer forth and 
back to the field modes, so it returns periodically to the atoms. 
Some examples of such coherent oscillations of entanglement 
will be discussed in details in Sections 8 and 9. 



C. Sudden birth of entanglement 

We have already seen that entanglement can exhibit strange 
features, sudden death and revival, which are characteristic of 
the dynamics of a special class of initial two-photon entan- 
gled states. In this section, we discuss an another noteworthy 
strange feature of transient entanglement, a delayed creation 
{sudden birth) of entanglement | S^, ■ The phenomenon 
arises for a different set of initial states of the system, and is 
characteristic of the transient dynamics of one-photon entan- 
gled states. As in the Sec. 5.1, we study the time evolution 
of the concurrence of a two-atom system initially prepared in 
a separable state. The model differs from that discussed in 
Sec. 5.1 in the initial state, where we considered only the case 
of the initial one-photon state l^s). The system, of course, 
can be prepared in other separable states. Here, we choose 
the separable two-photon state [^'4) as the initial state of the 
system and consider the transient spontaneous dynamics of 
the concurrence. The model is interesting in that it provides 
an example of the creation of entanglement on demand in the 
presence of a dissipative environment. 

Suppose that ait — the system is prepared in the separa- 
ble two-photon state |'I'4). In this case, ^44(0) = 1 and all 
the remaining density matrix elements are zero. According 
to Eq. ( l58T l, the off-diagonal terms (coherences) will remain 
zero for all time, but the populations Pss{t), Paa{t) and pii(t) 
will buildup during the evolution. This implies that for < > 0, 
the density matrix of the system spanned in the basis of the 
collective states (fTSl ). is of a diagonal form 



with the density matrix elements given by 



Pit) 



Pii{t) 

pss{t) 

paa{t) 

P44(<) 



(86) 
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andpii(t) = 1 - Pii{t) - pss{t) - Paa{t). 

After these preliminaries, we now turn to determine which of 
the two criteria, Ci (t) and/or C2 (<) could be positive. Looking 
at the density matrix ( |86] |, one can easily find that in this case 
Ci (t) is always negative. Thus, positive values of the concur- 
rence are determined only by the criterion C2{t): 



C = 



:{0,C2(i)} 



with 



C2{t) = \pssit) - Paait)\ - 2 y/ pn{t) Piiit) ■ 



(88) 



(89) 



We now use Eq. (188b to discuss the ability of the system to 
generate entanglement in the process of spontaneous emis- 
sion. It is clear from Eq. (|89] l that entanglement, if any, may 
only result from an unequal population of the symmetric and 
antisymmetric states. When the system is prepared in the state 
1 5*4), the resulting spontaneous transitions are cascades: The 
system decays first to the intermediate states and |^'a), 
from which then decays to the ground state Since the 

transition rates to and from the states |5's) and l^tf^) are dif- 
ferent when 712 7^ 0, there appears unbalanced population 
distribution between these states. According to Eq. ( [89] l, this 
may result in a transient entanglement between the atoms. 
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FIG. 7: The time evolution of the concurrence and its dependence on 
the distance between two, initially inverted qubits and fi _L f 12. 

To visualize the behavior of the concurrence, we plot C{t) 
as a function of time and the distance between the atoms in 
Fig. I2I We see that there is no entanglemnet at earlier times, 
but suddenly at some finite time an entanglement emerges. 
However, it happens only for a limited range of the distances 
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ri2. It is easily verified that the "islands" of entanglement 
seen in Fig. [7] appear at distances for which 712 is different 
from zero. Again, it reflects the role of the collective damping 
which is responsible for creation of entanglement by spon- 
taneous emission from an initial separable state. This result 
also demonstrates how the collective damping leading to en- 
tanglement between the atoms becomes less important as the 
interatomic distance is increased. Atoms separated by more 
than just a few wave-lengths become separable. 
We have seen that 712 modifies the damping rates of the tran- 
sitions between the collective states and, on the other hand, 
introduces two time scales for the decay of the population of 
the system. One can easily find that the entanglement seen 
in Fig. I2] decays out on a time scale (7 — 712)^^ that is the 
time scale of the population decay from the antisymmetric 
state. Again, this shows that crucial for entanglement creation 
by spontaneous emission is the presence of the antisymmet- 
ric state. This is perhaps not surprising, as the population of 
the antisymmetric state builds up on a much longer time scale 
than the population of the symmetric state. At long times, the 
antisymmetric state will posse a large population with with 
no population left in the symmetric state. What is surprising 
and, in fact, is very similar to the situation found in the Dicke 
model that at early times 7i <C 1 and 712 7^ 0, the symmetric 
state is significantly populated with almost no population in 
the antisymmetric state, and no entanglement is created. 




FIG. 8: The time evolution of the populations p44(t) (solid line), 
pas{t) (dashed line) and paa{t) (dash-dotted line) for 9 — 0, 
ri2/A = 0.25 and jl 1. fi2. 

The above considerations are supported by the analysis of the 
time evolution of the population of the excited states of the 
system that is illustrated in Fig. [8] It is quite evident from the 
figure that at the time t ss 4/7 when the entanglement starts 
to build up, the antisymmetric state is the only excited state of 
the system being populated. Clearly, the effect of the delayed 
creation of entanglement is attributed to the slow decay rate 
of the antisymmetric state. The state decays on the time scale 
of (7 — 712 that is much shorter than the decay time of the 
symmetric and the upper states. 

Concluding this section, we note that sudden birth of entangle- 



ment is characteristic of initial states which include the two- 
photon state with both qubits inverted. The sudden death of 
entanglement is not found in the small sample Dicke model 
that ignores the evolution of the antisymmetric state. It is also 
not found in a system with initially only one qubit excited. 
In experimental practice, the initial conditions of both atoms 
inverted that can be done by using a standard technique of 
a short tt pulse excitation. One could also use a short 7r/2 
pulse excitation which leaves atoms separable and simultane- 
ously prepared in the superposition of their energy states that 
includes the state with both atoms inverted. 



VII. TWO ATOMS IN A SINGLE-MODE CAVITY 

The systems we have considered so far involved two atoms 
interacting in free space with a common environment that ap- 
pears as multi-mode reservoirs to the atoms. In this section, 
we consider a different scheme where the atoms are confined 
to a single-mode cavity, as illustrated in Fig. |9] It is reason- 
able to expect that inside the cavity the atomic dynamics will 
be strongly affected | HIM [Ml- 

We discuss the entanglement of two identical atoms which, in 
general, could be in nonequivalent positions inside the cavity 
field 1 90]. We model this situation by assuming that the atoms 
experience different amplitudes of the cavity field which, on 
the other side, leads to different coupling strengths of the 
atoms to the cavity field. We derive the master equation for 
the density operator of the atoms by using the technique of 
adiabatic elimination of the cavity mode. We work in the 
limit of a large detuning of the cavity field frequency from the 
atomic resonance and/or a large cavity damping. With such 
an approximation, we can eliminate the dynamics of the cav- 
ity mode leaving the master equation for the atoms alone. The 
analogy of this system with that of a single two-level atom 
driven by a detuned coherent field is exploited and discussed. 
Next, we discuss in detail entangled properties of the system 
and the dependence of the entanglement on the evolution time 
and the position of the atoms inside the cavity mode. 




FIG. 9: Schematic diagram of the system composed of two atoms 
coupled to a single-mode standing-wave cavity. Each atom is located 
in a separate half- wavelength of the cavity mode. 

We begin by a brief outline of the derivation of the master 
equation for the density operator of two atoms coupled to a 
cavity mode. As before, the atoms are modeled as two-level 
system (qubits) with upper levels \ei), [i = 1, 2), lower lev- 
els \gi), connected by an electric dipole transition and sepa- 
rated by energy hwQ. The atoms are coupled to a standing- 
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wave cavity mode with the coupling constants .g(r'i), hav- 
ing units of frequency and being proportional to the electric 
dipole matrix elements of the atoms. In addition, the atoms 
are damped at the rate 7 by spontaneous emission to modes 
other than the privileged cavity mode. Moreover, we assume 
that the atoms are separated by many optical wavelengths and 
therefore ignore the direct interaction between them, i.e. the 
collective parameters 712 — and fli2 — 0. This also allows 
a selective preparation of the atoms such that a given set of 
initial conditions for the atomic states is achieved. The cavity 
mode is damped with the rate k and its frequency uic is sig- 
nificantly detuned from the atomic transition frequency ljq, so 
there is no direct exchange of photons between the atoms and 
the cavity mode. The behavior of the total system, the atoms 
plus the cavity mode, is described by the density operator px, 
which in the interaction picture satisfies the master equation 



dpT 
dt 



'j[Ho,Pt\ - ]^nCcPT - ]pCaPT, 



where 



^0 = E Wj)aS'^ 



+ -iAt 



H.C.] 



(90) 



(91) 



is the Hamiltonian describing the interaction, in the rotating- 
wave (RWA) approximation, between the cavity field and the 
atoms, and 



CcPt — o) apT + pTd^a ~ lap^a 



l^<^PT = Y.^{ptS^. + {St S^PT^) ' (92) 

are the Liouvillian operators representing losses in the sys- 
tem by the cavity damping and by spontaneous emission, re- 
spectively. The operators S*^ k^^ ) ^'^'^ ^'"^ '^he usual rais- 
ing (lowering) and population difference operators of the j th 
atom, respectively, a and are the cavity-mode annihilation 
and creation operators, A = — t^^o is the detuning of the 
cavity-mode frequency from the atomic transition frequency, 
and Tj is the position coordinate of the jth atom within the 
cavity mode. 

To isolate the atomic dynamics, we introduce the photon num- 
ber representation for the density operator with respect to the 
cavity mode 



Pt 



E 



|m) (n| 



(93) 



where pmn are the density matrix elements in the basis of the 
photon number states of the cavity mode. To avoid population 
of higher levels we assume that the cavity mode is strongly 
detuned from the atomic transition frequency. In this case, 
we may assume that only the two lowest energy levels, the 
ground state (n = 0) and the one-photon state (n — 1) of the 
cavity mode can be populated. Thus, neglecting population in 
all but the two lowest levels, we find from the master equa- 
tion (|90] l the following equations of motion for the diagonal 



field-matrix elements, the population of the photon-number 
states 



Poo 



2 1 



J = l 
2 



Pii = -*E (djSj poi - gjPwSj) 



(94) 



and the off-diagonal density matrix elements, the coherences 

2 

Poi = 

1 



1 \ 

-K + i A j poi - « E 5j ('^j^^ii ~ PoaSj') 
^ .7=1 



-.ll^aPm 



Pio 



K - iA j pio + « E 3] {piiSj - Sj Poo) 



- -jl^^aPlO- 



(95) 



Here, gj = g{rj) and we have removed the fast oscillating 
terms by introducing a rotating frame through the relations 



Poi = Poi exp(iAi), pio = pio exp(-«At). 



(96) 



Equations ( |94l ) and ( |95] l are cumbersome, although involving 
only two equations, because of the coupling which exists be- 
tween the diagonal and off-diagonal matrix elements. More- 
over, the field-matrix elements are still operators with respect 
to the atomic variables. 

In order to simplify the procedure and also to minimize spon- 
taneous emission of photons to modes other than the privi- 
leged cavity mode, we assume that the cavity frequency is 
highly detuned from the atomic resonances and/or the cavity 
field is strongly damped, i.e. 



A,K > 7,gj. 



(97) 



Under such conditions, the most populated state of the cav- 
ity field is the ground state |0), and then we can make a fun- 
damental physical assumption: The coherence poi changes 
slowly in time, in the sense that poi ^ t^oPoi- so that we can 
apply the so-called adiabatic approximation, which states that 
for a large detuning and/or cavity damping 



poi w and pio « 0. 



(98) 



With such an approximation, the coherences obey the follow- 
ing relations 



^ "j^ XI {s^pii ~ POOS'; 
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Efj (P^^Sj ^ Sj Poo) ' 



(99) 
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where /C = k/2 — zA. 

Substitution of Eq. (l99T l into the equations of motion for the 
populations, poo and pn, leads to 



1 2 ^ ^ 

= ^ X! 3'9*j [sis J ,poo] - --yCaPoQ, (100) 



and 



P^^^ ^Yl 9i9j [Si Sj.pii] - --iLaPii- (101) 



Since 



Poo + Pii = Ttf{pt) = P 



(102) 



is the reduced density operator of the atoms, we find by adding 
Eqs. ( 1 1001 ) and ( IIOII ) and neglecting the population pu, as 
the cavity mode will never be populated, that the adiabatic 
elimination procedure leads to the following master equation 
for the density operator of the atoms 



dp 
'dt 



= iY.5, \StS- , p] + z E [StSj , p] 



i=l 
2 



where the coefficients of the equation are 



(103) 



A2 



(104) 



The master equation (1103b is the principal result of this sec- 
tion. The major importance is that it includes both, single- 
atom and two-atom collective terms. The coefficients and 
5i represent respectively the decay rate induced by the cav- 
ity and the frequency shift of the energy levels of the ith 
atom. The shift is an analog of a dynamic Stark shift jgOh . 
The multi-atom terms, and Vlij {i ^ j) represent, respec- 
tively, the collective damping and the collective shift in en- 
ergy separation of the levels of atom i due to its interaction 
with the atom j through the cavity mode. The term fiy is an 
analog of the familiar dipole-dipole interaction between the 
atoms 1^,111,12111^]. This shows that one can engineer the 
interaction between distant atoms by the adiabatic elimination 
of the cavity mode. It is interesting to note, that the engineer- 
ing depends on what king of the cavity is used. The Stark 
shift 5i and the collective shift depend on the detuning 
A and vanish when A = 0. The damping terms 7^ depend 
on the cavity damping n and vanish in the good cavity limit 
of gi ^ K, or in the limit of A k. Consequently, no ad- 
dition atomic damping and the collective damping result from 
the adiabatic elimination of the cavity mode. In other words. 



the interaction with a strongly detuned cavity mode does not 
create the collective incoherent interaction of the atoms with 
the environment. 

The new feature of non-equivalent positions of the atoms in- 
side the cavity mode is in the shift of the atomic levels 6i and 
in the additional damping rates that can have different val- 
ues for different atoms when gi 7^ g2- As a consequence 
of unequal coupling strengths, the system composed of two 
identical atoms located at different positions inside the cav- 
ity mode may behave as a system of two non-identical atoms 
of different transition frequencies and different damping rates. 
In what follows, we choose the reference frame such that 

5(^1) = 50 and g[r2) = 5o cos (fcri2) , (105) 

where ri2 ~ \r2 — | is the distance between the atoms. This 
choice of the reference frame corresponds to a situation where 
atom 1 is kept exactly at an antinode and the atom 2 can be 
moved through successive nodes and antinodes of the standing 
wave. This choice, of course, involves no loss of generality as 
the dynamics of the system depends on the relative rather than 
the individual positions of the atoms. 

The remaining term ^jCaP represents spontaneous emission 
of the atoms into modes other than the cavity mode. For atoms 
located inside a cavity, the damping rate 7 depends on the 
cavity mode function, also known as the Airy function of the 
cavity |[36l l9Tl - l93ll . The function depends on the solid an- 
gle subtended by the cavity mode, and is often modeled as a 
simple Lorentzian with a maximum value at the cavity res- 
onance frequency. The Lorentzian has a half width at half 
maximum of k. In other words, the cavity "tailors" the vac- 
uum modes surrounding the atoms to those orthogonal to the 
mirrors, which modifies the structure of modes available to 
the atoms for spontaneous emission, thus changing the spon- 
taneous emission rates. For an open cavity in which the cavity 
mode subtends a small solid angle, 7 is approximately equal 
to the free-space damping rate. For a small size cavity, where 
the cavity mode subtends a large solid angle, the rate 7 is very 
small and can be neglected. In this case, the atoms interchange 
photons entirely with the cavity mode and any losses in the 
system occurs only through the cavity muTors with the rate 
determined by the cavity damping k. 



A. Case of a good cavity 

Consider first the good cavity limit of A ^ ^ k, 7. Now 
because the cavity damping is negligible, the atoms are cou- 
pled only through the coherent dipole-dipole interaction term 
ilij and there is no the collective damping, 7^ — 0. The 
dynamics between the atoms will simplify to the coherent ex- 
change of an initial excitation. We shall explicitly demonstrate 
how the dynamics of the system can be expressed in terms of 
the Bloch vector components of a two-level system driven by 
a detuned coherent. Next, we will illustrate how the concur- 
rence C{t) and its dependence on the relative position ri2 of 
the atoms inside the cavity mode can be readily interpreted in 
terms of the variables of the two-level system. 
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As we have already learnt, the concurrence is determined in 
terms of matrix elements of the density operator of the system. 
Thus, we use the master equation ( 1102b and find that in the 
product-state representation (fT2] ). the equations of motion for 
the density matrix elements are in the form 

P44 = -27P44, pii = 7 {p22 + /O33) , 
P23 = - (7 - i^l2) P23 + i^l2 (P22 - P33) , 
P32 = - (7 + i5l2) P32 - i^l2 (P22 " P33) , 
P22 = -7P22 + i^l2 (P23 - P32) , 

P33 = -7P33 - i^l2 (P23 " P32) ■ (106) 

The parameter 612 = Si — 62 is a difference between the 
single-atom Stark shifts. This parameter is of central impor- 
tance here as it determines the relative variation of atomic 
transition frequencies with position of the atoms inside the 
cavity mode. Note that the parameter S12 is present only when 
<?i 92, and vanish for gi = 52- 

Equation ( 11061 ) can be converted into two separate sets of 
equations if we introduce the following linear combinations 

W = P23+P32, V ^ i{p23 ~ P32), 

W = P22-P33, S = p22 + P33, (107) 

In terms of the new variables, Eq. ( 11061 ) splits into a set of 
equations 

P44 = -27P44, pii = 7Pss, 

Pss = -IPss + "ijPM, (108) 

involving only the populations, and 

u = —-fu + S12V, 

V = —JV — S12U — 2fli2W, 

w = --fw + 2ni2V, (109) 

involving the coherences and the population difference be- 
tween the states |'I'2) and [^'3). 

It is interesting to note that the equations of motion ( |109t are 
of exactly the same form as the optical Bloch equations of 
a two-level system driven by a detuned coherent field [I94ll . 
Here, the dipole-dipole interaction ili2 plays the same role as 
the Rabi frequency of a driving coherent field, and S12 appears 
as a detuning of the field from the driven transition [94]. How- 
ever, the analogy is not absolute because, in contrast to the 
case of a "real" two-level atom driven by a detuned laser field, 
the parameters ili2 and 612 in Eq. ( |109t are associated with 
the same quantities gi and g2, i.e. both depend on the cou- 
pling constants of the atoms to the cavity field. Consequently, 
the parameters f7i2 and 612 can not be independently adjusted. 
Moreover, the damping rate 7 that appears in Eq. ( 11091 ) repre- 
sents spontaneous decay out of the two-level system, actually 
to the auxiliary state l^*!), not to the ground state of the two- 
level system. 

The dynamics of the effective two-level system can be easily 
visualized in the Bloch vector model, as shown in Fig.[TOl In 
this model, the system and the driving field are represented by 
vectors in a three-dimensional space, and the time evolution 



is simply visualized as a precession of the system-state vector 
about the driving field. In terms of the Bloch vector, Eq. ( 11091 ) 
can be written as 

^ = -^B + nBy<B, (110) 

at 

where fls = 2rii2, 0, ^12) is the pseudofield vector and 
B = {u,v,w) is the Bloch vector The quantities u and v 
are, respectively, the real and imaginary parts of the coherence 
between the levels [^'2) and |^'3), and w is the population 
inversion. The solutions of Eq. ( II lOl l describe the damped 
precession of the Bloch vector around (Ib- The Bloch vector 
makes a constant angle 6 — tan^^(— 2r2i2/(5i2) with (Ib, 
it rotates around it tracing out, in the absence of damping, a 
circle on the Bloch sphere. When the Bloch vector does not 
intersect the "north pole" B = (0,0, 1) or the "south pole" 
B = (0,0,-1), the inversion w{t) 7^ ±1 and then, as we 
shall see, an entanglement between the atoms occurs. 
Figure[TT]shows an example of the Rabi nutation due to a de- 
tuned field. The system was initially in the state |5'2). We see 
that the Bloch vector periodically intersects the "north pole", 
w{t) = 1, but never intersects the "south pole", w{t) — —1, 
when the atoms are not in the equivalent positions inside the 
cavity mode. This means that complete population inversion 
does not occur during the evolution of the system. Viewing 
this result from the perspective of the population transfer be- 
tween two atoms, no complete population transfer between 
the atoms can occur. 

Let us now consider how the concurrence varies in time and, 
in particular, with the relative position of the atoms inside the 
cavity mode. We assume that initially the system is prepared 
in an arbitrary one-photon state. In this case, it is enough 
to solve the set of equations ( 11091 ). and a simple calculation 
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FIG. 10: The Bloch vector model of a system composed of two prod- 
uct states. The vector B rotates with frequency on the sphere be- 
tween the north pole corresponding to the excited state and the south 
pole corresponding to the ground state of the system. Any other point 
corresponds to a superposition of these states. 
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shows that he solutions, vaHd for arbitrary initial conditions, 
are of the form 

,,, 2Q12A + 5i2(voa sin at + B cos at) 
u{t) = e ^\ 



the concurrence can be written as 



( vqu cos at + B sin at) ^ 
v[t) = e ' , 



a 



(111) 



,^ — 6 12 A + 2ni2 ivQa sin at + B COS at) 
w{t) = e ^\ 



where a — AVl\2 + (5^2 is the detuned Rabi frequency, and 

^ = 2O12U0 - (5i2Wo, B ^ 5i2Ua + 2^i2Wq, (112) 

are constants given in terms of the initial values wq = w{0), 
uq = u(0) and vq = v{0) of the population inversion and co- 
herence in the system. The density matrix elements oscillate 
in time with the detuned Rabi frequency a, and all are equally 
damped with the rate 7 due to the decay of the population to 
the ground state 

Since the u{t) and v{t) components of the coherence are re- 
lated to the interaction between the atoms, their properties 
should be reflected in the entanglement between the atoms. 
For an initial single-quantum excitation, the density matrix el- 
ements pi4{t) and p44{t) are always zero and then, according 
to Eq. (|29] |. the concurrence is solely determined by the crite- 
rion C2{t): 



C{t)=C2{t)=2\p23{t)\. 



(113) 



Non-zero values of C{t) are recognized as the signature of 
entanglement between the atoms. We see from Eq. ( I113l l 
that the basic dynamical mechanism for entanglement cre- 
ation in this system is the coherence P2?,{t)- That is, the origin 
of the entanglement in the system can be traced back to the 
time evolution of the coherence P23{t)- Utilizing the relation 
P23{t) — {u{t) — iv{t))/2, and through the relation 



u\t)+v^{t)+w\t) 



- o-27t 



(114) 




FIG. 11: Population inversion w{t) as a function of time for 7 = 
and different displacements Srd of the atom 2 from an antinode of 
the standing wave: Sra = O.IA (solid line), Srd ~ 0.18A (dashed 
line). The system was initially in the state |^'2). 



Cit) = y/l-w^{t) e" 



7* 



where the following notation has been introduced 



—5i2A + 2Qi2 {voa sin at + B cos at) 



(115) 



(116) 



Equation (II 15b gives us a transparent way to interpret the en- 
tangled dynamics of the system. It shows that the concurrence 
can be completely determined only by evaluating the contri- 
bution of the population inversion. When the population is in 
the lower or upper level, w{t) = ±1, and then C{t) = 0, 
whereas C{t) achieves its optimum value C{t) = 1 when 
u;(i) — 0. Therefore, we can interpret the entanglement as 
a consequence of a distribution of the population among the 
energy levels. 

It is easy to see from Eq. ( II 15b that the evolution of the con- 
currence depends crucially on whether gi = 32 or gi §2- If 
91 — 92, then S12 = and one finds from Eq. ( II 15b that the 



concurrence vanishes periodically at the particular times 



to — mr/a, 



0,1,2,.. 



(117) 



On the other hand, if gi 7^ 32 then 612 7^ 0, and we find that 
the concurrence behaves quite differently such that it vanishes 
only at times 



ts — 2tQ — 2mr/a, n = 0, 1, 2, 



(118) 



which shows that the entanglement exsists on the the time 
scale twice as long as for the case of equal coupling constants. 
This rather surprising result can be understood in terms of spa- 
tial localization of the energy emitted by an initially excited 
atom. For equal coupling strengths 612 = because the en- 
ergy levels of the atoms are equally shifted by the interaction 
with the cavity mode. In this case the energy emitted by the 
first atom oscillates with frequency 2rii2 such that at times 
tn = mr/a {n = 0, 1, 2, . . .) it is fully absorbed by the sec- 
ond atom. Since at these times the energy is well localized in 
space as being completely absorbed by the localized atoms, 
the entanglement, which results from an unlocalized energy, 
is zero. The situation changes when gi 7^ 52- In this case 
the energy levels of the atoms are unequally shifted by the 
interaction with the cavity mode. Consequently, due to the 
frequency mismatch, the energy emitted by the atom 1 is not 
fully absorbed by the atom 2, leading to a partial spatial delo- 
calization of the excitation at discrete times t — mr/a, where 
n = 1, 3, 5, . . .. Consequently, at these times a partial entan- 
glement is observed. The entanglement vanishes every time 
the excitation returns to its initial state, i.e. when it returns to 
atom 1. 

Above considerations are illustrated in Fig. [121 which shows 
the concurrence as a function of time and the position ri2 of 
the atom 2 between two successive nodes of the standing wave 
cavity field. We see oscillations in the transient evolution of 
the entanglement that follow the Rabi flopping of population 
back and forth between the atoms. In other words, this os- 
cillation reflects nutation of the atomic populations which, in 
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turn, can be associated with the precession of the Bloch vector 
B about the driving field vector Hb with frequency a. The 
most interesting feature of the transient entanglement seen 
in Fig.[T2]is that in the case of unequal coupling constants, the 
initially unentangled system evolves into an entangled state, 
and remains in this state longer than for the case of equal cou- 
pling constants. 



cosine shape, i.e. the shape of the cavity mode function. It 
is given by the product of two terms, one the absolute value 
of the cavity mode function | cos kr i2\ and the other the time- 
and position-dependent diffraction structure. In other words, 
the concurrence is in the form of position and time depen- 
dent diffraction pattern. 



B. Entanglement diffraction pattern 



Let us now briefly consider the sensitivity of the entanglement 
on the relative position of the atoms inside the cavity mode. 
We use the solution (11091 ) to evaluate the concurrence as a 
function of ri2 for a fixed time t. We will discuss the variation 
of the concurrence with ri2 for two initial conditions. Our 
examples of the initial conditions are firstly the separable state 
l^'a), and secondly, the maximally entangled symmetric state 

Our first example considers the initial separable state IvPs). In 
this case, we readily find from Eq. (II 15b that the variation of 
the concurrence with position of the atoms inside the standing 
wave is given by 

X |cosfcri2|e-^^. (119) 

Here, rf = (1 + cos^ fcri2)T/2, we have make use of a scaled 
time variable r = (2(7q/A) t and dimensionless quantity F = 
(A/2.go^) 7- 

The concurrence depends on time, the relative distance be- 
tween the atoms and decays exponentially in time with the 
rate F. It exhibits an interesting modulation of the amplitude 
of the harmonic oscillation. Somewhat surprisingly our re- 
sults show that, in general, the concurrence does not have a 




FIG. 12: Concurrence as a function of time t and the position ri2 
of the atom 2 between two successive nodes of the standing wave 
cavity field for 7 = 0.1. The system was initially in the separable 
state 1*2). 
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FIG. 13: Variation of the concurrence C(r 12) with the position of the 
atom 2 inside a half-wavelength of the cavity field for 7 = and 
different times r: (a) r = 7r/2, (b) r = In /2, (c) r = 3l7r/2. The 
system was initially in the state j^&a) (wq = 1,uo = = Q). 



Figure [T3] shows the concurrence as a function of ri2 for 
different times r corresponding to the evolution intervals at 
which the entanglement is maximal for the idealized case of 
gi = 92- At early times, the entanglement is seen to occur 
over the entire range of ri2 inside the half-wavelength of the 
cavity field. In addition, the concurrence is a bell-shaped func- 
tion of position without any oscillation. As time progresses, 
oscillations emerge and consequently the region of ri2 where 
the optimum entanglement occurs, becomes narrower The 
evolution of the concurrence tends to become increasingly os- 
cillatory with ri2 as time increases, and the optimum entan- 
glement occurs in a still more restricted range of ri2. As a 
result, the entanglement oscillates with position faster than 
the cosine function, and the oscillations are more dramatic for 
larger times. Only at very early times (r <C 1), the oscilla- 
tions are not modulated by the diffraction pattern and the con- 
currence reduces to | cos/cri2|, but for longer times, C{r 12) 
may vary slower or faster than the cosine functions. Another 
interesting observation is that within the diffraction structure 
itself, the magnitude of the concurrence exhibits a succession 
of nodes and of antinodes. As a consequence, the entangle- 
ment may be completely quenched even for locations of the 
atom close to the antinode of the cavity mode, and alterna- 
tively may achieve its optimum value even for locations of the 
atom close to a node of the cavity mode. 

In the second case, we assume that initially the system was 
prepared in the maximally entangled state \^s)- Now, the 
general solution ( II 151 ) leads to the concurrence that can be 
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written as 



C{ri2) 



sin id 



{kr 12) cos^ikr 12) 



-Ft 



(120) 

The concurrence exhibits only the second order diffraction 
pattern. This may result in a weaker dependence of the con- 
currence on the position of the atoms than that observed with 
the initial unentangled state. Figure[T4]shows how the concur- 
rence evolves with the position of the atom 2 inside the cav- 
ity mode and time t, in the absence of spontaneous emission 
(7 = 0). Notice that the initial entanglement is reduced for 
positions of the atom different from the antinode of the cavity 
mode. That is, at some points outside the antinode the concur- 
rence may become zero indicating the complete reduction of 
the initial entanglement. 




FIG. 14: Variation of the concurrence C(ri2) with the position of the 
atom 2 inside a half-wavelength of the cavity field for 7 = and 
different times r: (a) r = 7r/2, (b) r = 7tv/2, (c) t = 3l7r/2. The 
system was initially in the state |^'s) (uq — 1,vo = wq = 0). 

The strong modulation of the concurrence in time, seen in Fig- 
ures [13] and [14] can be understood as a consequence of the 
uncertainty relation between the evolution time and energy, 
AEAt > h. For the increasing time At, the uncertainty of 
the energy AE decreases which means that the energy (exci- 
tation) becomes more localized. The increase in the localiza- 
tion of the energy results in a degradation of the entanglement. 
In other words, with an increasing time, one can in principle 
obtain more information about the localization of the atoms 
inside the cavity mode. 

In concluding this section, let us determine the explicit ana- 
lytical form of single photon entangled states of the system. 
The calculations can be easy done using the density matrix of 
the system, which in the basis of the product states ( fT2l i has 
the form 



Pit) 



Pii[t) 

P23(<) 

P33(<) 





(121) 



where we have assumed that only one photon is present in the 
system, i.e. pn{t) = for all times. 

We see that the density matrix (11211) is not diagonal in the ba- 
sis of the product states ( fT2l i due to the presence of the coher- 
ences P2?,{t) and P32{t)- Thus, the product states (fT2] i are not 
the eigenstates of the system. To obtain the actual eigenstates 

We 



of the system, we have to diagonalize the matrix (1121b 
first find the eigenvalues of the matrix, from which we next 
determine eigenvectors, i.e. eigenstates of the system. Pro- 
ceeding in the usual way, we calculate the determinant of the 
matrix p{t), and find a characteristic (polynomial) equation 

X[X-pn{t)]{\^-[p22{t)+P33{t)]X 
+P22(t)p33W- |P23(t)P} =0. (122) 

This is a simple fourth order polynomial whose the roots, cor- 
responding to the eigenvalues of the matrix (11211) . are 



Xiit) = piiit), Xiit) = 0, 

A2,3W = l[p22it)+P33{t)]±^U, 

Where U = [(p22(i) - P33(t))' + 4|p23(<)P] 
following relations are satisfied for all times 

P22{t) + P33{t) = 1 - Pll(t), 
Mt)- P33it)f =W^t), 
4\p23{t)\^=V^t)+u\t), 



(123) 
Since the 



and 



we can write the roots (11231 ) in the form 

Ai(t) = pii(t), Xi{t)^0, 
A2,3W-i[l-pii(i)±e-^*] 



(124) 



(125) 



(126) 



It is seen from Eq. (1126b that in the presence of spontaneous 
emission (7 7^ 0), three of the eigenvalues are non-zero, indi- 
cating that the system is in a mixed state. In the Bloch picture, 
this corresponds to the fact that the Bloch vector B is not con- 



stant. According to Eq. ( 1125b . the magnitude of the Bloch vec- 
tor decays in time with the rate 27. However, in the absence 
of spontaneous emission, 7 = 0, and then the population of 
the ground state pii{t) = for all times. Consequently, all 
the eigenvalues are zero except X2{t), that is equal to one. In 
this case, the magnitude of the Bloch vector is conserved and 
the system remains in a pure state for all times. 
Following the standard procedure of diagonalisation of matri- 
ces, we can easily find the eigenstates of the matrix ( 1121b cor- 
responding to the eigenvalues (1123b . The normalised eigen- 
states are of the form 



\Mt)) = piiimi), 

|^2(t)> =;^[p23(t)|*2) 
l^3(t)> -;^[p23(t)|*2> 



(A2W-p22(t))|*3)], 
(A3(i)-P22(<))|*3)], 

(127) 
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where AAj,3 - {|p23(0l' + [A2,3(i) - /022(t)]}'^'. Thus, the 
system initially prepared in a one-photon state evolves into a 
mixed state determined by the states (11271) . 
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FIG. 15: Magnitude of the coherence P2sit) as a function of the 
normalised time at for 7 = and different ri2: ri2 = 0.3A (solid 
line), (b) ri2 — 0.4A (dashed line), ri2 = 0.5A (dashed-dotted line). 
The system was initially in the separable state |*I'2). 



It is evident from Eq. ( I127l i that the states |i?2(t)) and |t?3(t)) 
are in the form of entangled states that are produced by the 
coherence P23{t)- Thus, a non-zero value of the coherence 
P27i{t) leads to entangled states of the atoms. In other words, 
whenever p23 [t) 7^ 0, the atoms are entangled. 

In Fig. [15] we plot | p23 (i) | as a function of time for 7 = and 
different relative positions ri2 of the atoms inside the cavity 
mode. The figure demonstrates that the coherence develops 
in time and remains different from zero for all times except 
for some discrete times at which the population is completely 
in the states \^2) or |5'3). Interestingly, when the atom 2 
is not located precisely at an antinode of the cavity field, the 
coherence remains non-zero for longer times indicating that in 
this case the atoms can be entangled on the time scale much 
longer than that predicted for the atoms precisely located at 
the antinodes of the cavity field. 

As above, we can explained this result in terms of the degree 
of localization or non-locality of the energy. One can see from 
Eq. (1106b that for equal coupling strengths, 5i2 = 0, and then 
the initial energy oscillates with frequency 2fii2 between the 
atoms. When absorbed by one of the atoms, the energy is then 
well localized in space as being at localised objects. Conse- 
quently, the entanglement, which results from a non-locality 
of the energy, is zero. The situation changes when 612 7^ 0. 
According to Eq. ( IIO6I 1, in this case the energy levels of the 
atoms are unequally shifted. Due to the frequency mismatch, 
the energy radiated by the atom 1 is not fully absorbed by the 
atom 2, leading to a partial spatial delocalization of the exci- 
tation at discrete times corresponding to the time at which the 
atom 2 is excited. Consequently, at these times an entangle- 
ment is observed. 



C. Case of a bad cavity 

In the preceding section, we have discussed the dynamics of 
two atoms coupled to a single-mode cavity field. We have 
considered the good cavity regime, where we have ignored 
the effects of the cavity damping. In the present section, we 
continue the analysis of the dynamics of the atoms coupled to 
a single-mode cavity, but now we consider the opposite limit, 
the bad cavity regime where the dynamics are primary deter- 
mined by the cavity damping k. As we shall see, a fast damp- 
ing of the cavity mode may result in a dark state geometry of 
the system that there might exist an entangled non-decaying 
(trapped) state. We demonstrate the occurrence of a station- 
ary entanglement in the system. 

If we include the cavity damping, k 7^ 0, we then readily 
find from the master equation ( 1103b that the density matrix 
elements satisfy a set of four coupled differential equations 
that can be written in a matrix form as 



= AY, 

dt 



(128) 



where F is a column vector having components 
(P22, P33, u, v), and A is a 4 x 4 matrix 



V 



-72 


-712 

-2^12 





-71 

~7l2 

2ni2 



'2 712 
-2 712 

-S12 



-i^l2 
^12 
-5l2 

-h ) 



(129) 



with 5i2 ^ S2~Si,^ ~ 71 +72, and the remaining coefficient 



given in Eq. ( 1104b . 

The equations of motion ( 1128b are first-order differential equa- 
tions with the coefficients dependent on the relative position 
of the atoms inside the standing-wave cavity mode. The set 
of equations (1128b is easily solved by matrix inversion. The 
method, however, depends critically on whether det[A] 7^ 
or det[A] = 0. It is easily verified that the determinant of the 
matrix A is of the form 



det[A] = rii2(72 - 71) +7i2'5i2- 



(130) 



Further, using the explicit forms of the coefficients, Eq. ( 11041 ). 
we find the determinant is equal to zero. It is always zero 
independent of the values of the parameters involved. An im- 
mediate consequence of det[A] = is that the density matrix 
elements can have non-zero steady-state values that depend on 
the initial values of the density matrix elements |95-97]. This 
fact is connected with the existence of a linear combination of 
the density matrix elements 



Pec = 71P22 + 72P33 " 712M, 



(131) 



which is a constant of motion, i.e. pcc = 0. 

The existence of the constant of motion reflects the presence 

of a linear combination of the states of the system which does 



not decay. From Eq. ( 1128b . it is straightforward to show that 
the state corresponding to the constant of motion is of the form 



72 



71 + 72 



1*2 



71 



71 +72 



1^-3), (132) 
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It is easily verified that the state ( 11321 ) satisfies the trapping 
condition of pcc = 0. In other words, spontaneous emission 
from this state does not take place, so that any initial popula- 
tion in the state \^c) will stay there for all times. 
However, the channel |'I'4) — > Y^c) can be open for spon- 
taneous emission, i.e. transitions into the state \^c) from the 
two-photon state |4'4) could be allowed for spontaneous emis- 
sion. Hence, the state |5'c) could be populated by spontaneous 
emission from the upper state [^'4). It is easy to find from the 
master equation ( 1103b that the inclusion of the state |'I'4) into 
the dynamics of the system leads to the following equation of 
motion for the population of the state |^'c): 



( 72 - 71) ' 
71 +72 



-P44- 



(133) 



We see that the state is populated with the rate propor- 
tional to the difference between the damping rates induced 
by the cavity damping. According to Eq. ( 11041 ). the damping 
rates are different from each other only when gi 7^ 52, i-e. 
when the atoms are in nonequivalent positions inside the cav- 
ity mode. Otherwise, the transition rate from the state |'I'4) to 
I'I'c) is equal to zero indicating that the state becomes decou- 
pled from all the states so that cannot be accessible by spon- 
taneous emission. In this case, the state can be regarded 
as a dark state in the sense that the state is decoupled from the 
environment. 




FIG. 16: Concurrence as a function of time t and the position ri2 of 
the atom 2 between two successive nodes of the standing wave cavity 
field for 7 = and A/k = 30. The system was initially in the state 

Let us now to compute the concurrence with the initial state 
|/) corresponding to an excitation of the system into the sep- 
arable state |^'2). On making use of Eqs ( l29b and (I128I I. and 
remembering that with the initial state \^2), the coherence 
PiA{t) and the population pn{t) are always zero, we obtain 

2gi52 



g2g-^(7l+72)(l-i»?)t 



where = A/k. 

If the system is initially prepared in the state l^'s), the concur- 
rence takes the following form 



C|*3)(t) 



2.91.92 

g2g-i(7l+72)(l-i'7)t 
322-3(71 +72)(l+i»))i 



.g2_g2^-(7i-f72)* 



(135) 



The most obvious characteristic of these results is that the 
concurrence depends on the initial state even in the station- 
ary limit of < ^ 00. For the steady-state, we have 

2^/?52 



C|*2>(i ^ 00) - , — -2 , 

W +.92) 

H*3>(* ^ 00) -> — -2 . 

(gi + gi) 



(136) 



+ ff2e" 



-^(7l+72)(l+ir))t 



(134) 



Evidently, the system decays to a stationary state which de- 
pends on the initial conditions. Only in the limit of gi — g2, 
which corresponds to the atoms located at the equivalent po- 
sitions inside the cavity mode, the concurrence decays to 
the same stationary value, C|^2^(t 00) = C^^^^•j{t — 
00) — 1/2. Otherwise, when gi 7^ g2, the concurrence de- 
cays to different stationary values, which strongly depend on 
whether gi > 92 or 91 < 92. 

Figures [16] and [17] show the time evolution of the entangle- 
ment and its dependence on the position of the atom 2 inside a 
half-wavelength of the standing-wave cavity mode. The most 
interesting feature is that the entanglement decays to a non- 
zero stationary state which depends on both, the relative po- 
sition of the atoms inside the cavity mode and the initial state 
of the system. One can see from the figures that the stationary 
entanglement can be larger for positions of the atom different 
from the antinode of the cavity mode. Again, this effect can 
be explained in terms of the frequency shift of the atomic lev- 
els that leads to the frequency mismatch resulting to a spatial 
delocalization of the initial excitation. 

So far we have considered only the evolution of the con- 
currence from initial one-photon separable states. We have 
demonstrated that a system of two atoms located inside a 
standing wave cavity mode and initially prepared in a sepa- 
rable one-photon state, can decay to an entangled state. The 
amount of entanglement depends on the relative position of 
the atoms inside the standing wave mode. In fact, the amount 
of the stationary entanglement depends on the initial popula- 
tion of the trapping state |5'c). Since the initial population of 
the state \'^c) cannot day, it stays there for all times. 
A somehow different situation arises if the system is initially 
prepared in the separable state |^'4). In this case, the initial 
population of the state \'^c) is equal to zero. Nevertheless, the 
system could decay to a stationary entangled state determined 
by the population of the trapping state |5'c). This possibility 
arises from the fact that under some circumstances the transi- 
tion channel [^'4) — > \'^c) can be allowed for spontaneous 
emission. Using Eqs. ( i29l ) and ( 11031 ). it is straightforward 
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to show that for the initial state |/) = {'^4), the steady-state 
value for the concurrence is given by the expression 

C{\^i),t ^00)^ p^^^. (137) 

(71+72)'' 

It is readily seen that the atoms can be entangled in the steady- 
state. However, it happens only if 71 7^ 72, i.e., when the 
atoms are located in nonequivalent positions inside the cavity 
mode (gi 7^ 92)- Otherwise, when 71 =72, which corre- 
sponds to gi ~ §2, the system decays to a separable state. 

VIII. TWO ATOMS IN SEPARATE CAVITIES 

Our third model for creation of entanglement between two 
atoms consists of two separated subsystems, each composed 
of a single mode cavity containing a single two-level atom [ 491 
Us, 86, 98]. In general, the two subsystems may not be identi- 
cal in that the cavity frequencies and the coupling constants 
of the atoms to the cavity modes could be different. We 
study the entanglement evolution of two remote atoms in- 
teracting independently with a cavity field, as in the double 
Jaynes-Cummings model. We consider two cases: the first 
where there is only one excitation present in the entire sys- 
tem; the second where there is a superposition of two excita- 
tion levels, zero excitation and two excitations, one in each 
of the two remote systems. These two cases exhibit com- 
pletely different properties of the entanglement evolution. As 
we shall see, a complete and abrupt loss of atom-atom en- 
tanglement, the entanglement sudden death, is found to occur 
for the two-photon Bell states provided the initial entangle- 
ment is not maximum, though the full entanglement is peri- 
odically regained. We present complete entanglement solu- 
tions, for the case of asymmetric atom-cavity couplings and 
off-resonant interactions. We will show that asymmetric cavi- 
ties and cavity-atom detunings can prove advantageous in en- 




FIG. 17; Concurrence as a function of time t and the position ri2 
of the atom 2 between two successive nodes of the standing wave 
cavity field for 7 = and A/k = 30. The system was initially in 
the state 1*3). 



abling a control of entanglement. For example, entanglement 
sudden death can be eliminated with the use of suitable de- 
tunings of the atomic transition frequencies from the cavity 
mode frequencies. Thus the general effect of detuning is to 
stabilize entanglement. The situations of asymmetric atom- 
cavity couplings and off-resonant interactions correspond to 
realistic experimental situations, where it may be difficult to 
produce identical cavities. The unequal coupling constants for 
example may arise when atoms are not in equivalent positions 
inside the cavities. 

We begin by introducing the Hamiltonian of the system, which 
now is composed of two separated single mode cavities each 
containing a single two-level atom, as shown in Fig. [18] The 
total Hamiltonian H for the atoms and the cavity fields, in 
the electric-dipole and rotating-wave approximations, can be 
written in the Jaynes-Cummings (JC) form as 

H = HF + HA + H,nt, (138) 

where 

Hp = huJi (^a\ai + + fiwz ^a^a2 + (139) 

is the Hamiltonian of the cavity fields. 

Ha = nwoSl + hwoS^ (140) 
is the Hamiltonian of the atoms, and 

Hint = hgi{a\S^ + aiSl ^ +%2 (oaS'a" +a25^^ (141) 

is the interaction Hamiltonian between the atoms and the cav- 
ity modes. Here, gi {i — 1,2) is the coupling constant be- 
tween the ith atom and the field mode of the ithe cavity. As 
usual, we have denoted by 5"^^ and the raising and low- 
ering operators of the ith atom, and by aj (aj) the creation 
(annihilation) operators of the mode of the jth cavity. 




SUBSYSTEM 1 SUBSYSTEM 2 



FIG. 18: Schematic diagram of a system composed of two distant and 
non-interacting subsystems 1 and 2. Each subsystem is composed of 
a single-mode cavity containing a two-level atom. 

In contrast to what we have considered in the above sections, 
we now proceed to analyze the dynamics of the system that is 
completely isolated from the environment, so that the dissipa- 
tion effects are not present. On the one side, this restricts the 
applicability of the results to the case of a good cavity. In fact, 
this is not an overly restrictive limitation regarding the recent 
advances in the cavity QED. On the other hand, the absence 
of the dissipation means that the evolution of the system is 
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purely coherent and is governed by the Schrodinger equation. 
It also admits analytical treatment of the concurrence. 
We consider separately the cases of single and two-photon ex- 
citations of the system. In the first, we assume that only a 
single excitation is present and find the state vector of the sys- 
tem by solving the Schrodinger equation 



(142) 



where 



leW) =di(t)|a)+d2WI6)+d3WI6)+d4Wie4) (143) 

is the state vector of the system given in terms of the product 
states ( l40l l of the atoms and the cavity modes at time t. The 
coefficient di {t) determines the probability amplitude of the 
ith state at time t. The problem is to find the time evolution of 
the probability amplitudes provided that the initial state of the 
system is known. Of particular interest are the values of di (t) 
for an imperfect matching of the atoms to the cavity fields. It 



is easily verified that with the Hamiltonian (1138b . the coeffi- 
cients satisfy the following differential equations 



di 
d2 



-ig2di, 



da 
dn 



2iAid3 
2zA2d4 



igidi, 
ig2d2, 



(144) 



where 2Aj = (cjq — Wj) represents the detuning between the 
atomic transition frequency and the jth cavity frequency. 
Equations ( 11441 ) form decoupled pairs of simple differential 
equations that can be solved by using e.g. the Laplace trans- 
form technique. A solution of the equations, valid for an arbi- 
trary initial state is readily found to be 

di(t) = e'"^'' {di{0) cos{nit) 

-z[5idi(0)+/3id3(0)]sin(17ii)}, 



d2{t) 



d3{t) 



e 



iA2t 



{d2(Q) cos(f}2i) 



i[S2d2{0)+l32dimM^2t)}., 



{dsiO) cosiflit) 



i[(Sid3(0)-/3idi(0)]sin(l]iO}, 



d4{t) ^ e'"^'' {d4{0) cos{n2t) 



+ i [(52^4(0) - M2{0)] sm{n2t)} , (145) 



where fti = \/^| + A|^ {i = 1,2) is a detuned Rabi fre- 
quency, dj (0) (j = 1, 2, 3, 4) are the initial values of the prob- 
ability amplitudes, 6i — Ai/ili, and (3i — gi/fti represent, re- 
spectively, the scaled (dimensionless) frequency detuning and 
coupling constant between the ith atom and the ith cavity. 
The probability amplitudes oscillate sinusoidally with the 
Rabi frequency fti, and their dynamics is strongly affected 
by the modulation term, proportional to sm{ilit), whose am- 
plitude depends on the detuning Si and the coupling con- 
stant (3j between the atom and the corresponding cavity mode. 
The modulation terms vanish periodically at times t„ — 
mr/rii, n = 0, 1, . . .. 

In general, the time evolution of the probability amplitudes is 
quite complicated and not easy to interpret. However, one can 



see from Eqs. ( 1145b that the detuning enters the solutions in 



an antisymmetric way, whereas the coupling strength enters 
the solutions in a symmetric way. As we shall see, this dif- 
ference will be evident in the features of the time evolution of 
entanglement in the system. 

Our interest is to evaluate the concurrence to quantify entan- 
glement between the different parts of the two two-qubit sub- 
systems 1 and 2. By denoting the two atoms as A and B, and 
the corresponding cavity modes as a and b, we may distin- 
guish six pairs of sub-systems AB, ab, Aa, Ab, Ba and Bb. 
We shall quantify entanglement in each pair by the concur- 
rence. To calculate the entanglement between any two qubit 
pair ij, for example the atoms A and B, we take the trace over 
the other sub-sytems, to evaluate the reduced density matrix, 
from which the concurrence can be calculated. In general, the 
reduced states are two qubit mixtures, and have a density ma- 
trix of the general X-form 



I ait) 

b{t) z{t) 

z*{t) c{t) 

V d{t) 



(146) 



which, according to Eq. ( I29b gives for the concurrence 



C,y (t) = 2max{0,|z(i)| - ^a{t)d{t)]. 
These concurrences are found to be evaluated as 



(147) 



CAB{i) = 2\dimd2{t)l Cabit) = 2|d3Wp4(t)|, 
CAa{t) = 2|di(t)||d3(t)|, Casit) = 2 1 ^2 ) | W I , 
CAb{t) = 2\di{t)\\di{t)\, Csbit) = 2\d2{t)\\di{t)\. 

(148) 

Note that there is no the threshold term when only one-photon 
states are involved in the dynamics of the system ll99l] . 
As an example of the sensitivity of the system to imperfect 
matching of the atoms to the cavity modes, consider the evo- 
lution of the concurrences for the following initial state 



1 

72 
1 

72 



(16) ±16)) 

(|ei) I52) ± Iffi) lea)) ® |0)i 10)2 



(149) 



Note, that the initial entanglement is restricted to just the 
two atoms. Hence, for the case of equal coupling constants, 
gi — g2 = g but non-zero detunings, the bipartite concurrence 
measures are 



CAB{t) — cos^ i^gt) - 
Cabit) = -^sin^ngt) 



ain^ (ilgt) 



CAait) = CAbit) = CBait) - Csbit) = - sin 



cos^ i^gt) 



sin^ i^gt) 



(150) 
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where we have introduced a dimensionless Rabi frequency, 
the same for both subsystems, ft — {l + S^ Y^^, and the scaled 
detuning S = A/g. 
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FIG. 19: Time evolution of the pair concurrences for equal cou- 
pling constants, gi ~ 92 ~ g and exact resonance Ai — 
A2 — 0: CAB{t) (solid line), Cab{t) (dashed line), and 
CAa{t),CaB{t),CAbit),CbBit) (dashed-dotted line). The atoms 
were initially in the entangled state \'i/s). 

The expressions dlSOb for the concurrence measures are plot- 
ted in Fig. [191 We see periodic oscillations of the concur- 
rence, which reveal the periodic transfer of entanglement from 
atoms to cavity modes, and vice versa. It is interesting that the 
transfer process of entanglement from the atoms to the cav- 
ity modes does not involve just the pairs CAsit) and Cabit)- 
By inspection of the time evolution of the concurrences in 
Fig. [19] we find that the initial maximal entanglement be- 
tween the atoms is not only totally transfeiTed to the cavity 
modes, but at the same time an additional pairwise entangle- 
ment is created during the evolution. This can also be seen by 
summing the pair concurrence measures to find that at times 
t„ = mr/Afl, where n = 1, 3, 5..., the total pair concurrence 
is larger than one. As time progresses, a part of the entan- 
glement is transferred into the other pairs of the sub-systems. 
Then, after a further interval, the entanglement is completely 
transfeiTed into the cavity modes. The additional entangle- 
ment vanishes at times i„ = mr/2n, n ~ 1,2, 3, when the 
transfer process is completed. 

Another interesting prediction of ( llSOI l is that complete entan- 
glement transfer from the atoms to the cavity fields requires 
exact resonances. To examine the effect of detuning, we plot 
in Fig. |20]the atom-atom concurrence CAB{t) as a function 
of time and the detuning A. It is evident from the figure that 
in the case of non-zero detuning A, the localized initial en- 
tanglement is not completely transferred to another pair of 
qubits. Detuning increases the oscillation frequency and de- 
creases the minimum entanglement between the atoms. It is 
interesting, however, that the atomic entanglement returns to 
its initial maximum value periodically with the detuned Rabi 
frequency. This can be understood in terms of the localiza- 
tion of the initial energy. When A ^ 0, only a part of the 
initial energy, proportional to ljq is transferred to the cavity 
modes leaving the excess energy unlocalized. Thus, we may 



conclude that a non-resonant coupling of the atoms to the cav- 
ity modes can be used to stabilize entanglement between the 
atoms. 




FIG. 20: The atom-atom concurrence Cab (t) plotted as a function 
of the normalized time gt and detuning A/g for equal coupling con- 
stants, gi = g2 = g. 

We now turn to the problem of entanglement dynamics when 
two excitations (one in each subsystem) are present in the sys- 
tem. Using the Schrodinger equation, we set up equations of 
motion for the probability amplitudes and then solve them to 
study the dynamics of the entanglement. We define the state 
vector of the system in terms of the basis states ( |4TI ) as 

\x) = di{t)\xi)+d2{t)\x2)+d3{t)\x3)+d4{t)\x4), (151) 

and find from the Schrodinger equation that the probability 
amplitudes satisfy a closed set of four coupled equations of 
motion 

di = -i(Ai + A2)di - ^(32^2 + 51^3), 

d2 = -«(Ai - A2)d2 - i(.g2rfi + 51^4), 

da = i(Ai - A2) J3 - iigidi + .92^4), 

di = i{Ai + A2) J4 - j(.9irf2 + .924)- (152) 

To remove the fast oscillating terms in Eq. (I152l i. we have 
introduced a rotating frame through the relations 

J,; =d,c-*("i+"^)*, i = 1,2,3,4. (153) 

As before for the single excitation, we allow for the possibil- 
ity of imperfect couplings of the atoms to the cavity modes 
by introducing detunings Ai and A2, and unequal coupling 
constants gi and 92- 

Although the set of equations of motion ( I152l l involves four 
coupled equations, it is analytically solvable. We can put 
Eqs. ( 1152b into a matrix form and solve them by the matrix 
inversion. However, due to the complexity of the general so- 
lution, we study separately the time evolution of the proba- 
bility amplitudes for two cases: In the first case, we con- 
sider non-zero detunings but equal coupling constants, i.e. 
Ai = A2 = A and gi = 92 = .9- In the second, we con- 
sider zero detunings, Ai = A2 = but unequal coupling 
constants, gi 7^ 52- 
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We now investigate the evolution of the bipartite entanglement 
in the six pairs of the subsystems. We stress that the general 
state of the system, Eq. ( 11511 ). cannot give any entanglement 
in the atom-atom subsystem, regardless of the values of the 
coefficients involved. This is readily seen by evaluating the 
reduced density operator for the atom systems. By tracing the 
density operator of the total system over the cavity modes, we 
obtain the reduced density operator of the atoms alone 

PAB =Tr cavity ^ {n\{m\p\m)\n) , (154) 

m.n— 0, 1 

where \n) and |m) refer to cavity modes a and b respectively, 
and find that 

PAB = Ml(i)n4'4)(*4| + M2(i)n4'2)(*2| 

+ M3wn*3)(*3i+M4(i)n*i)(*ii- (155) 

Evidently, the density operator represents is a mixture of 
separable states and hence cannot be entangled. The atom- 
atom bipartite entanglement is always zero, Cab(0 = 0. 
The same conclusion applies to the bipartite entanglements 
Cab{t),CAbit) and Csait)- The only non-zero concurrence 
possible is for the qubit pairs CAa{t) and Csbit)- In other 
words, an entanglement can be created during the evolution 
between the atoms and the cavity modes to which they are 
coupled. The amount of entanglement that can be created in 
these pairs can be determined from the relations 

CAa{t) - 2\di{t)dl{t) + d2{t)dl{t)l (156) 

and 

CBb{t) - 2\d^{t)dl{t) + d:i{t)dl{t)\. (157) 

Equations show that the source of entanglement between the 
atoms and the corresponding cavity modes is in the interfer- 
ence between the probability amplitudes. 




FIG. 21: Time evolution of the concurrence CAait) for the initial 
separable state and different detunings: A = (solid line) and 
A = g (dashed line). 

Figure |2T] shows the entanglement creation and evolution be- 
tween the atom A and the cavity mode a, with an initial sep- 
arable state Ixi). Similar entanglement properties appear be- 
tween the atom B and the cavity mode h. We see that an en- 
tanglement between the atom and the cavity mode is created 



during the transfer process of the excitation between them. 
The maximum entanglement is observed at times when the 
excitation is equally shared between the atom and the cavity 
mode, and vanishes when the excitation is completely located 
at either the atom or the cavity mode. In terms of the Rabi 
oscillations, the system becomes disentangled at every half of 
the Rabi cycle of the oscillation, i.e. at qt = mr/Q, {n — 
0,1,2,...). 

Similar to the case of single excitation states, the system be- 
comes entangled for a longer time when the frequencies of 
the atom and the cavity field are detuned from each other. In 
other words, an imperfect matching between the atom and the 
cavity field leads to a more stable entanglement than in the 
case of the perfect matching. In terms of the Rabi oscillations, 
the system becomes disentangled only at every Rabi cycle of 
the oscillation, i.e. at gt — 2mi/Q,. This is easy to under- 
stand if one considers entanglement as resulting from a super- 
position corresponding to a delocalization of energy. When 
the frequencies of the atom and the cavity mode are different, 
the initial amount of energy localized in the atom is not com- 
pletely transferred to the cavity mode. A part of the energy, 
not absorbed by the cavity mode remains delocalized, which 
results in a nonzero entanglement. 

The above analysis has showed that no pairwise entanglement 
is possible between any parts of the two subsystems 1 and 2 
when we have with certainty a single excitation in each of the 
subsystems. Nevertheless, we can generate entanglement in 
the qubit pairs by including in the wave function an auxiliary 
state \xo), the ground state for which there is no excitation, 
and consider the time evolution of a state vector of the system 
containing the auxiliary state 

\x{t)) = dlWIXl>+rf2(i)|X2>+d3WIX3) 

+ d4WIX4) +do(i)lxo). (158) 

With the auxiliary state included, we find that the bipar- 
tite concurrence measures between the subsystems are of the 
forms 

CAB{t) = 2max{0, \di{t)\\do{t)\ - 1^2(^)11^3 W|}, 
Cab{t) = 2max{0,|d4(t)Po(<)| - \d2(t)\\d^{t)\] , 
CAb{t) = 2max{0,|d2(i)||do(t)| - Mi(t)||d4(t)|}, 
CaB{t) = 2max{0, \d^{t)\\d^{t)\ - \di{t)\\di{t)\] . 

(159) 

Evidently, an entanglement between the pairs of subsystems 
is possible only if the ground state is included. Since the state 
Ixo) is not an eigenstate of the Hamiltonian of the system, the 
probability amplitude d[){t) is a constant of motion, i.e. the 
amplitude does not evolve in time, do(t) = do(0). 
Equation ( I159l l shows several interesting properties of the 
concurrence and differences between the single and double- 
excitation cases. Firstly, the concurrences appear as differ- 
ences of products of the absolute values of the probability 
amplitudes, i.e. they involve the threshold term for entangle- 
ment. This gives a possibility for sudden death and revival 
in the behavior of entanglement that were absent in the case 
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of the one-excitation. Secondly, there is evident a competi- 
tion in the creation of atom-atom entanglement between two 
pairs of states (|xi), |X4)) and (IX2), Ixs))- For example, en- 
tanglement creation in concurrence pairs involving the states 
Ixi) and 1x4) is diminished by the presence of population in 
the states |x2) and jxa), and vice versa creation of entangle- 
ment involving the states |x2) and jxs) is diminished by the 
presence of population in the states |xi) and |x4). In terms of 
the population transfer, entanglement creation by a simulta- 
neous exchange of two photons is diminished by one-photon 
exchange processes, and vice versa. The competition between 
these one and two-photon processes is the source of the phe- 
nomenon of entanglement sudden death. 




FIG. 22: The atom-atom concurrence Cab as a function of the nor- 
malized time gt and detuning A/g for gi — g2 = g and the initial 
state Ixsd), given by Eq. J43t with a = 7r/12 and l3 — 0. The effect 
of the detuning is to remove the sudden death of entanglement. 

Figure |22] illustrates the atom-atom concurrence as a function 
of time and the detuning A, calculated from an initial con- 
dition of the atoms prepared in the non-maximally entangled 
state ( |43] ). We can see that in the case of A = 0, the ini- 
tial entanglement is lost in a discontinuous fashion and can 
revive periodically during the evolution. On the other hand, 
when A 7^ 0, the amplitude of the oscillation decreases and 
the entanglement sudden death dissapears. Thus, an imperfect 
matching of the cavities to the atomic transition can prevent 
the onset of the entanglement sudden death, to provide a sta- 
bilization of entanglement. The principle of the effect, which 
is similar to that showed in Fig. |20l may again be readily be 
understood by realizing that a nonzero detuning inhibits the 
fuU entanglement transfer from the atoms to the cavity modes. 



IX. TRIGGERED ENTANGLEMENT EVOLUTION 

We now proceed to analyze the problem of how one could 
switch on (trigger) an evolution of stable or "frozen" entan- 
glement. This is an important problem for controlled trans- 
mission on demand of an initial entanglement stored, for ex- 
ample, in a trapped state. As we have seen, trapped states have 
the property that an initial population stored in these states 
cannot decay due to zero transition dipole moments to and 



from the remaining states of a given system. Then an interest- 
ing practical question arises, how one could access the trap- 
ping state to trigger an evolution of the encoded information. 
In the following, we address this question and consider two 
schemes introduced in the previous sections. In the first, we 
consider two atoms located inside a single mode cavity and 
initially prepared in a maximally entangled state. In the sec- 
ond, the atoms are located in separate cavities and the system 
is initially prepared in an entangled state such that the initial 
entanglement is equally shared between all the pairs of the 
sub-systems. The concurrence is then monitored as a function 
of time and atom-cavity-field detunings. 
Let us begin with the scheme, introduced in Section 7, that 
involves two atoms located inside a single mode cavity. Sup- 
pose that the atoms are initially prepared in the collective sin- 
gle excitation symmetric state 



1 



\9i) 



(160) 



which is a maximally entangled state. Alternatively, we may 
say that an information has been initially encoded into the sys- 
tem and is stored in the maximally entangled symmetric state. 
For the initial state ( 1160b . the initial values of the Bloch vec- 
tor components are uo = 1, vq = wq ~ 0. It follows from 
Eqs. (II 1 lb and ( II 15b that initially at i = 0, the concurrence 
C(0) = 1, and then for t > the time evolution of the concur- 



rence is of the form 



C{t) = 



.^12 . 
I — sm at 

a 



.-It 



, (161) 



where 7 is the spontaneous emission rate and 5i2 results from 
a non-equivalent position of the atoms inside the cavity mode. 
We see from Eq. ( 11611 ) that in the absence of the spontaneous 
emission and with 612 = 0, the concurrence does not evolve in 
time. In other words, the concurrence is a constant of motion 
that an initial entanglement encoded into the system at i = 
will remain there for all times. The evolution of the concur- 
rence can be triggered by switching on a non-zero detuning 
6i2- In this case, the concurrence oscillates in time and the 
amplitude of the oscillation increases with 612- 




FIG. 23: Concurrence as a function of the normalized time at and 
the position ri2 of the second atoms inside a half-wavelength of the 
standing-wave cavity mode. The overall slow decay of the concur- 
rence is due to spontaneous emission with the rate 7/0 = 0.1. 
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Figure |23] displays the concurrence as a function of time and 
the position ri2 of the atom "2" inside a half-wavelength of 
the cavity mode. We see that the initial entanglement does 
not evolve in time when the atom is located at ri2/A = 0.5, 
the antinode of the cavity field. This situation corresponds to 
the case of both atoms located in equivalent positions, antin- 
odes of the standing wave of the cavity field. Consequently, 
9i = 32 and then di2 = 0. The absence of the oscillations in 
the concurrence versus time is linked to the fact that the Bloch 
vector B and the il b vector are initially aligned at t = and 
remain aligned for all times when ri2/X = 0.5. By break- 
ing the symmetry between the coupling constants, i.e. by a 
dislocation of the atom "2" from the antinode of the standing 
wave, gi ^ g2, and then 612 7^ 0. As a result, the concur- 
rence oscillates in time indicating a continuous oscillation of 
the excitation between the atoms. Here, the oscillation of the 
concurrence comes from the precession of the Bloch vector 
about the vector Ob with frequency a. Thus, by an imperfect 
coupling of the atoms to the cavity mode, one can trigger an 
evolution of the "frozen" entanglement. 
Consider now the second scheme, introduced in Section 8 that 
involves two atoms each located in separate single-mode cav- 
ity. We concentrate on the gi = 52 case and will demonstrate 
the effect of non-zero detunings Ai = A2 = A on the evo- 
lution of an initial entanglement. We use the same strategy as 
in the above example and assume that initially the system was 
prepared in a superposition state such that the concurrence re- 
mains a constant of motion under the condition that the modes 
of the two cavities are resonant with the atomic transition fre- 
quencies. 



It is easy to show from Eqs. (I145I I and ( I148l l that with the 
initial state (I162l i and zero detunings, the two-qubit entangle- 
ment between all of the pairs of the sub-systems remains con- 
stant in time. In other words, the initial entanglement is frozen 
in the initial state. It is interesting to note that the entangle- 
ment remains frozen in the initial state independent of whether 
91 = 92 orgi ^ 52. 

The dependence of the entanglement of an arbitrary qubit pair 
on time t and the detuning A is illustrated in Fig. |24] The 
evolution of the concurrence exhibits several interesting prop- 
erties. First of all, the concurrences vary in time only for 
nonzero detuning. Thus, a nonzero detuning A triggers an 
evolution of the entanglement. Secondly, the time evolution 
of the concurrences is not symmetric with respect to the sign 
of the detunings. A large and even maximal entanglement can 
be created when the detuning is positive, whereas the initial 
atomic entanglement is reduced and can even be suppressed 
when the detuning is negative. When, for example, the atom- 
atom entanglement is maximum, CAB{t) — 1, then the entan- 
glement between the other qubit pairs is zero and vice versa, 
when the atom-atom entanglement is zero then entanglement 
between one of the other qubit pairs is maximal. This result 
also implies a possibility of a controlled evolution of the sys- 
tem towards the maximum entanglement between one of the 
qubit pairs. Finally, the results show that the maximum en- 
tanglement of an arbitrary qubit pair can be created when the 
detuning is matched to the coupling constant, IS. — g. Other- 
wise, the entanglement is reduced. 

X. STEERED ENTANGLEMENT TRANSFER 
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FIG. 24: Concurrence as a function of the normalized time gt and 
the detuning A of the cavity modes from the atomic resonances for 
51 = 92- 

Assume that the system is initially prepared in a superposition 
state with a uniform population distribution over the available 
energy states 

l^o) = i(|ei)+c^'l6)±(l6)-e^1C4))), (162) 

where 9 and are arbitrary phase factors. In this case, 
|dj(0)| = 1/2 (i = 1, . . . ,4) and then, according to Eq. (fl48] l. 
the initial entanglement is equally shared between the six 
qubit pairs, so C„/j(0) = Cab(O) = Caa(O) = C„b(0) = 

Cpa=Cpb{Q) = 1/2. 



The possibility of triggering an asymmetric evolution of an 
initial entanglement by changing the sign of the detuning A 
implies that one can also engineer the direction of evolution 
of entanglement by controlling the detuning. By this we mean 
that an initial entanglement can be transferred to a desired "lo- 
calized" pair of the qubits by a suitable choice of the detun- 
ings. The localized atom-atom entanglement in this context 
means that the entanglement exists solely between the two 
atoms. As it is clear from Fig. |24l a suitably chosen posi- 
tive value of the detuning can channel entanglement entirely 
into the atoms. Therefore, the entanglement transfer can be 
controlled by varying the frequency of the cavity mode. 
Another method of transferring an entanglement between 
qubit pairs in a controlled way is to create an asymmetry be- 
tween the coupling constants gi and 92- We will describe how 
to achieve a steered evolution of an initial entanglement to a 
desired pair of qubits using the scheme involving two cavities 
each containing a single two-level atom. We assume that the 
system is prepared initially in a superposition state 

l^o) = ^(iei> + l6)), (163) 

where and 1^2) are the single excitation states of the ba- 
sis (|40] |. The initial state involves the superposition of the 
atoms only, so that initially at t — the atoms are maxi- 
mally entangled, Cab{^) = 1, and the other pairs of qubits 
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are disentangled. During the evolution, the atom-atom entan- 
glement can be completely transferred only to one of the pairs 
of qubits or can be redistributed between different pairs. A 
couple questions then arises; Can we steer the transfer pro- 
cess such that the entanglement could be transferred to a de- 
sired pair of qubits? Can we engineer the transfer process by 
changing only the system's parameters, such as the coupling 
constants of the atoms to the cavity modes? 
Let us proceed to answer these questions by considering the 
time evolution of the concurrence. Provided the system is pre- 
pared initially in the state |^o)> we find from Eqs. ( 11451 ) and 
(1148b that for the case of exact resonances Ai = A2 = 
and unequal coupling constants gi ^ g2, the entanglement 
measures take the forms 



CAB{t) = 


cos (fligt) 


|cos {fi2gt) 


Cab{t) = 


sin (fligt) 1 


sin (n2gt) 1 


CAa{t) = 


cos (fligt) 


|sin (Qigt) | 


CAb{t) = 


cos {^igt) 


|sin {Q,2gt) 1 


Csait) = 


cos {^2gt) 


|sin [VLigt) 1 


CBbit) = 


sin ($72.9^) 1 


cos {^l2gt) 



(164) 



Here we have introduced dimensionless Rabi frequencies 

0.1 = (I + u/g) and = (1 - u/g), with g = {gi + 52)/2 
and u= {gi- ,g2)/2. 

We see that the entanglement measures are in the form of the 
product of two harmonic functions oscillating with different 
frequencies. At t = 0, Ca_b(0) = 1 and it is our purpose to 
determine which of the qubit pairs can be maximally entan- 
gled at t > 0. Since sin (i7i(2).9^) and cos (i7i(2)3^) cannot 
be simultaneously equal to one, it rules out the possibility of 
the maximum concurrence equal to unity for the pairs Aa and 
Bh. To further distinguish which of the remaining four con- 
currence measures can be equal to one at t > 0, we note that if 
the ratio 172/51 is not an integer number or a fraction of an in- 
teger number, none of the concurrence measures can be equal 
to one at i > 0. It implies that no complete transfer of the 
initial maximal entanglement is possible to any of the qubit 
pairs. The complete transfer is possible only if the ratio is an 
integer number or a fraction of an integer number. 
However, the destination to where the initial entanglement can 
be completely transfered depends on whether the ratio is an 
even or an odd integer number. If the ratio 52 / ffi is an even 
number, the initial maximal entanglement between the atoms 
can be completely transferred only to the atom-field qubit pair 
Ba. On the other hand, if the ratio is an odd integer number, 
the initial entanglement between the atoms can be completely 
transferred only to the field-field qubit pair ah. From this it 
follows that 



Cas(O) = 1 



CBa{t) 
Cab{t) 



for 32/51 even, 
for 52/51 odd. 



(165) 



We illustrate this situation in Fig. |25] where we plot the con- 
currence of the different qubit pairs as a function of time 
for exact resonances but unequal coupling constants. When 
51 — 52, we obtain the complete transfer of the entangle- 
ment only between the AB and ab qubit pairs. For gi 7^ g2. 



the initial entanglement between the atoms can be completely 
transferred between AB and the qubit pairs ab and Ba. The 
transfer of the entanglement between AB and the other pairs 
also occurs, but is not complete. 




FIG. 25: Concurrence for qubit pairs plotted as a function of the 
normalized time gt = \{gi + g2)t for A = and different ratios 
of the coupling constants 32/31: (a) (72/51 ~ 1, (b) 32/51 ~ 2, 
(c) 52/51 ~ 3, (d) 52/51 ~ 4. In all figures the solid line is for the 
atom-atom concurrence CAB{t). The dashed line is for the concur- 
rence measure (a) Cab{t), (b) Csait), (c) Cab{t) and (d) Csait). 



The reason for this feature of the entanglement transfer can 
be understood intuitively by noting that, for example, for 
52/51 = 2 the Rabi frequency g2 of the population oscilla- 
tion in the cavity system 2 is twice that of the Rabi frequency 
gi for the population oscillation in the system 1. This means 
that over a complete Rabi cycle g2t = tt, the initial popula- 
tion in the system 2 returns to the atom, but at the same time 
the population makes a half Rabi cycle in the system 1, i.e. 
the excitation in system 1 will be in the cavity mode. Thus, 
CBa{t) = 1 at that time, with the concurrence in the other 
qubit pairs equal to zero. 

In summary of this section, we have seen that in the scheme 
involving two separate cavities, each containing a single two- 
level atom, we can achieve a controlled (steered) transfer of 
an initial entanglement to a desired pair of qubits. The results 
show that the initial entanglement can be completely trans- 
ferred not to all but only to some of the qubit pairs. The com- 
plete steered transfer of entanglement can be done by properly 
adjusting the cavity frequencies to the atomic transition fre- 
quencies or by a proper engineering of the coupling constants 
between the atomic dipole moments and the cavity modes. 



XI. BEYOND THE ROTATING-WAVE APPROXIMATION 

The most interesting feature of transient entanglement is the 
phenomenon of entanglement sudden death. We have already 
seen that an entanglement encoded in spin correlated states 
can undergo an abrupt discontinuous evolution, while an ini- 
tial entanglement encoded in a spin anti-correlared (single ex- 
citation) state decays asymptotically in time without any dis- 
continuity. The lack of discontinuity has been referred to zero 
population of the two-photon state involved in the threshold 
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term for the concurrence. If initially only a single excitation is 
present in the system, the population of the two-photon state 
is always zero. As a result, no threshold term is present in 
the concurrence measure and consequently an initial entan- 
glement evolves continuously without any discontinuity ll99ll . 
However, the analysis were restricte d to the evolution of the 
system under the RWA Hamiltonian flOO*! . In this section, we 
return to the problem of the transient evolution of entangle- 
ment in the system of two atoms located inside a single mode 
cavity. In contrast to what we have discussed in Section 7, 
we now examine the transient evolution of entanglement be- 
yond the RWA. We shall demonstrate that under the non-RWA 
interaction, an entanglement initially encoded into a single ex- 
citation state can un derg o a discontinuity during the evolution 
of the system 1 101-103]. We interpret this result as a conse- 
quence of the principle of complementarity between the evo- 
lution time and energy. 

Let us reconsider the system composed of two identical two- 
level atoms located inside a single-mode standing-wave cav- 
ity. Our purpose is to compare the dynamics of an entangle- 
ment encoded initially in a spin anti-correlated state of the 
system under the RWA and non-RWA Hamiltonians. We as- 
sume that the atoms are strongly coupled to the cavity mode 
and consider the situation where the strength of the coupling 
is controlled by varying the position of the atoms inside the 
standing wave. Moreover, we will work in the short time 
regime of ujot « 1, which is much shorter than a typical spon- 
taneous emission time. Hence, we will neglect the sponta- 
neous emission from the atoms. 

The dynamics of the system are determined by the master 
equation of the density operator ps of the total, qubits plus 
the cavity field system 

-dt - -p^-^^^^P^i 

- -K.[a)aps+ psa)a-2apsa)) , (166) 

where k is the damping rate of the cavity mode, and 



stant varies with the position of the atoms inside the cavity 
mode as 



1 

^^nRWA = -zf^O ^ 5*1 + hWcah 



go sm 



(168) 



where n = L/X determines the length of the cavity in units of 
the cavity wavelength A, and we have assumed that the atoms 
are placed symmetrically about the center of the cavity such 
that ri + r2 = L and r2 — ri = d is the distance between the 
atoms, as illustrated in Fig. |26l When the atoms are close to 
antinodes of the cavity mode, d w (n — i)A, and then both 
atoms experience the peak coupling strength go to the cavity 
mode. We shall call this situation, a strong coupling regime. 
The coupling constants decrease with the displacement of the 
atoms from the antinodes and become very small, g{rj) « 
for distances d w (n — 1)A. In this case, the atoms are weakly 
coupled to the cavity mode, and we shall call this situation, a 
weak coupling regime. Thus, the strength of the coupling of 
the atoms to the cavity mode can be controlled by varying the 
positions of the atoms inside the standing wave. 
The interaction part of the Hamiltonian (1167b is in the non- 
RWA form. It contains both, the energy conserving terms S^^a 

and Sj'a^ as well as energy non-conserving terms S~a and 
Sj'a^ also known as the counter-rotating terms. The counter- 
rotating terms describe processes in which a photon is anni- 
hilated as the atom makes a downward transition, determined 
by the Sj'a term, or a photon is created as the atom makes an 
upward transition, determined by the S^a^ term. 





FIG. 26: Schematic diagram of the system composed of two distant 
atoms symmetrically located about the antinode of the standing-wave 
cavity mode. 



+ [9irj)SJa^ + g*ir,)aS^] (167) 

i=i 

is the non-RWA Hamiltonian of the system. The first term 
on the right-hand side of Eq. (11671 ) is the Hamiltonian of the 
atoms, the second term is the Hamiltonian of the cavity field, 
and the third term the interaction Hamiltonian between the 
atoms and the cavity field in the electric-dipole approxima- 
tion. The operators a and are the usual annihilation and 
creation operators of the cavity field, 5'J and 5*1 are the atomic 
spin operators for the jth atom, and lUc is the cavity mode fre- 
quency which, in general, can be different from the atomic 
transition frequency ojq. 

The parameter g{rj), which appears in Eq. (1167b is the cou- 
pling constant between the cavity mode and the jth atom lo- 
cated at a position rj along the cavity axis. The coupling con- 



The counter-rotating terms, although at a first look seem to be 
unphysical since could lead to a violation of energy conserva- 
tion, continue to create a debate on their significance and im- 
porta nce in the interaction of atoms with the electromagnetic 
field ill 041 - 1 lodl . since they have proven to produce important 
dynamical effects. Examples include the Bloch-Siegert shift 
of atomic frequencies |107], quantum chaos |108], bifurca- 
tions in the phase space 1109], a fine structure in the opti- 
cal Stern-Ger lach e ffect fl ICill . generation of photons from the 
vacuum field flTT'] and entanglement betw een at omic ensem- 
bles with no any initial excitation present ll 12ll . We should 
point out that the effects are notable for strong couplings, 
go « iOQ, and appear over short times, t < 1/ujq. 
Consider the time evolution of the concurrence when the dy- 
namics of the system is determined by the non-RWA Hamilto- 
nian. The required atomic density matrix elements are found 
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by solving numerically the master equation ( 11661 ) in the com- 
bined n photon basis of the cavity field and the product state 
basis of the atoms, Eq. (fT2l i. We compute the atomic density 
matrix elements by tracing the density operator ps over the 
cavity field assumed to be in the vacuum state |0). Provided 
that the atoms are prepared initially in the single excitation 
collective symmetric state |/) ~ l^'s), the concurrence is then 
determined by the Ci (t) criterion 

C{t) = max{0, Ci{t)} 

= 2max{o,|p23W|- VPii(*)P44(t)} , (169) 

It is apparent that the discontinuity or threshold behavior of 
the concurrence requires a nonzero population of both, the 
ground and the two-photon states of the system. For the ini- 
tial state l^'s), /3ii(0) ~ /O44(0) = 0, and it is predicted 
that (t) ~ for all times as there are no external sources 
present. In this case, the initial entanglement is expected to 
evolve asymptotically without any discontinuity. 




FIG. 27: Concurrence as a function of the dimensionless time ujt and 
the distance d/X between the atoms for go = iu, {u — uo) — O.Oluj, 
K = O.lij and n — 1/2. The atoms were prepared initially in the 
single excitation collective symmetric state l^/s). 

Let us examine if the expectations of the asymptotic evolution 
of the initial entanglement are manifested when the dynam- 
ics of the system undergoes the non-RWA Hamiltonian. Fig- 
ure |27] displays the time evolution of the concurrence for the 
initial single excitation collective symmetric state j^^) and 
for the peak coupling strength of the atoms to the cavity mode 
go ~ oj. We see that for small distances between the atoms, 
corresponding to a strong coupling of the atoms to the cav- 
ity mode, g{ri,2) ~ go, the initial entanglement undergoes 
the discontinuity, the sudden death behavior. The entangle- 
ment disappears quite rapidly over a very short time, ujt ^ 1. 
For large distances, where the coupling constants g{rj) are 
very small, the discontinuity disappears and the entanglement 
decays asymptotically in time, the behavior predicted before 
under the RWA. Clearly, the counter-rotating terms which be- 
come important in the strong coupling regime, are responsible 
for the entanglement sudden death. 



The discontinuity in the time evolution of entanglement, the 
sudden death effect seen in Fig. |27] seems puzzling at first, 
because it appears to contradict the predictions based upon a 
simple argument that with a single excitation present in the 
system, it is impossible to achieve the discontinuity as the 
population j044(t) = for all times. To resolve this prob- 
lem, we plot in Fig. |28]the time evolution of the population 
of the two-photon state |5'4) for the same parameters as in 
Fig. |27] with c? = 0. It is seen that initially unpopulated two- 
photon state |^'4) becomes populated in a comparatively short 
time. Thus, it is clear now that the entanglement sudden death 
is associated with the non-zero value of the population P44 (t) 
which is developed during the evolution governed by the non- 
RWA Hamiltonian. 
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FIG. 28: Time evolution of the population p44(t) calculated for both 
atoms located at the antinode of the standing wave {d = 0), with 
go = id, {uj — (jjo) — O.Olcj, K = O.loj and n — 1/2. The atoms 
were prepared initially in the single excitation collective symmetric 
state 



To explain more clearly the origin of the population of the 
two-photon state, we offer two complementary views of the 
underlying physics. The first is obtained by the derivation of 
the reduced density operator for the atoms using state vectors 
of the combined atoms plus the cavity field system. The sec- 
ond is provided by the principle of complementarity between 
the evolution time At and uncertainty AE in the energy of the 
system. 

Assume that there is only a single excitation in the system that 
evolves under the non-RWA Hamiltonian, that includes the 
counter-rotating terms aSj' and S^a^ . In this case, we must 
include the processes that do not strictly conserve excitation 
number, when one of the atoms goes to the excited (ground) 
state by emitting (absorbing) a photon. If these processes are 
included, then the Hilbert space of the system, the atoms plus 
the cavity field, should be spanned by six-state vector [^'2) <8) 
|0), |*3)'»|0), l^-i)®!!), |«'2)®|2), |«'3)®|2) and 1^-4)011). 
If we evaluate the reduced density operator p{t) for the atoms 



38 



by tracing over the cavity mode, we find 

2 

p{t) = Tr cavity [ps{t)] ^y^^{i\ps{t)\i) 

i=0 

= Pll(t)|*l)(^'l|+P22(t)|*2>(«'2| 
+ P33(i)|«'3)(*3| +/044WI*4)(*4| 
+ P32(t)|«'3)(*2| +P23WI*2)(«'3| 
+ pi4W|«'l)(*4|+P4l(i)|*4)(*l|, (170) 

where \i) refers to the cavity mode and pij{t) are popula- 
tions {i = j) of the atomic states and coherences {i ^ j) 
between them. It is evident from Eq. ( 1170b that the evolution 
of the atoms under the non-RWA Hamiltonian involves the 
two-photon state, that population of the state |\E'4) becomes 
possible. 

An alternative explanation of the physical origin of the popu- 
lation of the two-photon state is provided by the principle of 
complementarity between the evolution time At and uncer- 
tainty AE in the energy 

AtAE>h. (171) 

Clearly, for the evolution time of the order 1/lu, that has been 
considered here, the energy of the initial excitation is spread 
over on energy interval of order hu, the order required to 
achieve a non-zero population of the upper state \'i'4) of the 
two-atom system. 

In summary of this section, we have demonstrated a somewhat 
surprising result that an entanglement initially encoded into a 
single excitation state can undergo the phenomenon of sudden 
death. This non-intuitive behavior of the entanglement occurs 
in the very short time of the evolution and arises solely from 
the counter-rotating terms in the interaction Hamiltonian of 
the system. When the dynamics of the system are governed 
by the RWA Hamiltonian, the entanglement decays asymptot- 
ically without any discontinuity. A simple explanation of this 
unexpected feature is provided in terms of the quantum prop- 
erty of complementarity, which is manifested as a tradeoff be- 
tween the knowledge of energy versus the evolution time of 
the system. 

XII. ENTANGLING ATOMIC ENSEMBLES 

The phenomena described in the previous sections arise from 
entangled and disentangled dynamics of simple systems com- 
posed of only two atoms. Particular attention has been paid 
to the anomalous "strange" behavior of entanglement. Exact 
analytical results have been obtained which allowed for the 
complete understanding of the underlying physics. 
It has become apparent that a large scale quantum information 
processing and quantum computation are, in general, much 
more complex and require to create special entangled states 
involving a large number of strongly coupled atoms. Natu- 
rally, when more than two atoms are involved, the situation is 
much more complex, and requires completely different pro- 
cedures than those already considered. In this connection. 



we now turn to multi-atom systems and address in this and 
the following sections the question of a simple procedure for 
creation of entangled states in multi-atom systems, such as 
macroscopic ensembles of cold atoms. The reason for the in- 
terest in atomic ensembles is twofold. On the one side, atomic 
ensembles are macroscopic systems that are easily created in 
the laboratory. On the other hand, the collective behavior 
of the atoms enables to achieve almost a perfect coupling of 
the ensembles to external fields without the need to achieve a 
strong field-single atom coupling. 

Our main effort will be devoted to the study of the cre- 
ation of entanglement between M atomic ensembles. We 
shall make use a simple but powerful procedure, initiated 
by Gao-xiang Li ill 13ll . which constructs multi-mode entan- 
gled states from the vacuum by internal dynamical processes 
rather than being injected into the ensembles from external 
sources. The procedure is associated with the collective dy- 
namics of the atomic ensembles subjected to driving lasers 
of a suitably adjusted amplitudes and phases. It employs the 
Holstein-Primakoff representation of angular momentum op- 
erators, which expresses atomic systems in terms of orthog- 
onal collective bosonic modes, and then prepares separately 
each of the modes in a desired entangled state. 
The physical system we first consider for creation of multi- 
mode entangled states consists of M ~ 2 one-dimensional 
atomic ensembles trapped along the axis of a high-Q ring 
cavity, as illustrated in Fig. |29] The cavity is composed of 
three mirrors that create two mutually counter-propagating 
modes, called clockwise and anti-clockwise modes, to which 
the atoms are equally coupled. The cavity modes are degen- 
erate in frequency, i.e. a;+ = a;_ = uic- The cavity is damped 
with the rate k that is assumed small to achieve a high finesse 
at a relatively large size of the cavity. 




FIG. 29: Configuration of a ring cavity with two atomic ensembles, 1 
and 2, for the preparation of multi-mode entangled states. The cavity 
modes are damped with the same rate n. The driving laser fields are 
injected through the cavity mirrors and co-propagate with one of the 
cavity modes. 

The atomic ensembles contain a large number of identical 
four-level atoms, each composed of two stable ground states, 
|0j„), |ljri)' ^nd two excited states \ujn), \sjn), where the 
subscript jn labels jth atom of the nth ensemble. The ground 
state |0j„) of energy i?o = is coupled to the excited state 
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\sjn) by a laser field of the Rabi frequency r2s„ and frequency 
ujls that is detuned from the atomic transition frequency by 
As — {uJs — ^Ls), where lOs = Es/h and Eg is the energy of 
the state \sjn)- Similarly, the ground state \ljn) of energy 



El = huJi is coupled to the excited state 



of energy 



Eu = fu^u by an another laser field with the Rabi frequency 
and the angular frequency lolu that is detuned from the 
atomic transition \ljn) — \ujn) by Au — (w„ — wi — wlu)- 
The frequencies ujlu and ujls of the laser fields are matched 
close to the cavity frequency, so that the wave numbers of the 
laser fields ku and kg are approximated hy ku ~ kg ~ k. A 
schematic diagram of the atomic levels together with the con- 
figuration of the driving lasers and cavity couplings is shown 
in Fig. [30] 

To secure some simplification in notation, we shall assume 
that both ensembles contain the same number N of the atoms. 
The atoms interact with the cavity modes of degenerate fre- 
quencies ojc that simultaneously couple to the \ujn) \^jn) 
and \ sjn) |ljri) transitions, with the coupling strengths gun 
and gsn, respectively. We assume that the atom-field coupling 
strengths are uniform through the atomic ensembles. This 
is consistent with current experiments involving ring cavities 
and and large samples of trapped atoms [56]. Practical sizes 
of the atomic samples are ~ lO'^ nm that is smaller than the 
cavity mode radius wo = 130^m, and also are much smaller 
than the practical length of a single arm of the cavity. 
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FIG. 30: Atomic level scheme. Two laser fields of the Rabi fre- 
quencies Qu and f2s drive the atomic transitions \ljn) ^ l^jn) and 
\Ojn) — >■ I Sj„), respectively. The atomic transitions \ljn) \sjn) 
and |0j„) \ujn) are coupled to the single-mode cavity with the 
coupling strengths gun and gsn, respectively. 



In what follows, we will assume that the laser fields and the 
cavity frequencies are far detuned from the atomic transition 
frequencies. This will allow us to adiabatically eliminate the 
upper states of the atoms and obtain an effective dispersive 
type interaction of the driving lasers and the cavity modes with 
the atomic transitions. 



A. Effective Hamiltonian of the system 

The Hamiltonian of the system, in the rotating-wave approx- 
imation, and transformed into the interaction picture has the 
following form 



H — Ho + Hal + Hac, 



(172) 



where 



2 N 



Ho = fiiJc(a+a+ + aLa_) + ^ X! X! {^u\ujn){ujn\ 

n=l j=l 

+ UJs\Sjn){Sjn\ + Wi 1 1 -,„) (lj„ | } (173) 

is the free Hamiltonian of the atomic ensembles and the cavity 
modes, 

2 N 



Hal = {r!„„(x,„)e-^("-+'^"")Vjn)(lj„| 

ri=l j=l 

+ r!,„(a;j„)e-'("^-+*=")*|s,„)(Oj„|+H.c.} (174) 



is the interaction Hamiltonian between the atoms and the driv- 
ing fields, and 



2 N 



Hac = XI {9Lixjn)\ujn){0jn\a+ 

n=l j = l 

+9uni^jn)\Ujn) {Ojn\a- + fjtni^ 3n)\^ jn) {^3 



+9sniXjn)\Sjn){ljn\a^ +H.C.} 



(175) 



is the interaction Hamiltonian between the atoms and the two 
cavity modes. Here, a± and Oj. are the annihilation and cre- 
ation operators associated with the two counter-propagating 
modes of the cavity; clockwise (+) and anti-clockwise (— ) 
propagating modes, and k is the wave number of the cav- 
ity modes. The Hamiltonian retains only the energy con- 
seving terms, which play the dominant role in the interac- 
tion. We have denoted the energies of the atomic levels by 
huji {i = l,u,s) and have set the energy of the ground state 
\Qjn) equal to zero. The parameters ^lun{xjn) and ^lsn{xjn) 
are spatially varying Rabi frequencies of the driving laser 
fields, and (j>un, <t>sn are their phases. The coupling constants 
of the atomic transitions to the cavity fields, guni^jn) and 
gfj^{xjn), are also dependent on the atomic position. In what 
follows, we will use the plane traveling wave representation 
for the laser fields and the cavity modes, in which 



' 9sni-'^jn) = 5s 



(176) 



where we have assumed that the coupling constants g.an 
and gsn are the same for the two cavity modes. This is accept- 
able if the modes have the same polarization and geometry, 
which is feasible with the practical situations ll54l - l56ll . 
We now make an unitary transformation and few standard ap- 
proximations on the Hamiltonian ( 11721 ) to eliminate the ex- 
plicit time dependence and the excited states to neglect atomic 
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spontaneous emission and to obtain an effective two-level 
Raman-coupled Hamiltonian. Then, we will make the trans- 
formation of the atomic operators into the field (bosonic) rep- 
resentation. These operations will be done in the following 
steps. In the first step, we make the unitary transformation 
U = cxp{iHl)t/h), with 

2 N 

+ h'^^KUu + UJl)\Ujn){Ujn\ 



71=1 j=l 

^^Ls\Sjn){Sjn\ + ^1 1 Ijn) (Ijn I } i 



(177) 



and assume that the laser frequencies satisfy the resonance 
condition ujls — ^Lu = 2wi. This specifically chosen trans- 
formation allows us to eliminate the explicit time dependence 
of the Hamiltonian. 

In the second step, we introduce detunings of the laser fields 
from the atomic transition frequencies 



(wL-u + Wi), As 



(178) 



In order to eliminate spontaneous scattering of photons to 
modes other than the privileged cavity mode, we assume as- 
sume that the detunings are much larger than the Rabi frequen- 
cies, cavity coupling constants and the atomic spontaneous 
emission rates, i.e. 



(179) 



where 7^ and 7„ are the total spontaneous emission rates from 
the states \sjn) and |wj„), respectively. 
The assumption about large detunings allows us to perform 
the standard adiabatic elimination of the atomic excited states, 
\sjn), \ujn), and obtain an effective two-level Hamiltonian 
involving only the ground states of the atoms 
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+ /3.„e-^^=" (4^)t„^. + ji'^\y_)+H.c.}, (180) 

where 6c = uJc — [ujls — 1^1) is the detuning of the cavity 
frequency from the Raman coupling resonance. 
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are the spatially dependent macroscopic inversion and polar- 
ization operators, and 



2A 



2A 
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are the coupling strengths of the effective two-level system to 
the cavity modes. 

The first line of Eq. (II8OI 1 represents the free energy of the 
atomic ensembles. The second line represents an intensity de- 
pendent (Stark) shift of the atomic energy levels, and the re- 
maining two lines represent the interaction between the cavity 
field and the atomic ensembles. The three collective atomic 
operators, and which appear in Eq. ( 11801) . 

arise naturally for the position dependent atomic transition 
operators and appear in a cavity with two mutually counter- 
propagating modes. In the case of a single-mode cavity, the 
Hamiltonian involves only the operator ll42l l43i IstIi . 
The operators crj"' , o-j""*^ and cr^"'' are the standard Pauli 
spin— 1/2 operators, which satisfy the well known commuta- 
tion relations [crj"-*, CT^™''^] — 2a'"J^^ SjiSmn- However, the col- 
lective operators J^^^ , J^^l^ and ji"' do not in general satisfy 
the angular momentum commutation relations. The reason is 
in the presence of the phase factors exp{imkxjn), so that the 
commutation relations are satisfied only in the small sample 
limit of kxjn ^ 1, at which cxp{imkxjn) ~ 1. 



The essential feature of the effective Hamiltonian ( 1180b is 
that the atoms interact dispersively with the cavity mode. 
This means that the cavity mode will remain unpopulated 
during the evolution. Note the presence of nonlinear terms, 
a^jj, aJ~, which are analogous to the counter-rotating 
terms. These terms will play the crucial role in creation of 
a squeezed state between the atomic ensembles. 
The Hamiltonian (11801 ) is general in the parameter values. To 
avoid unessential complexity due to the presence of the free 
energy and the Stark shift terms, we will work with a simpli- 
fied version of the Hamiltonian by choosing the frequencies 
of the driving lasers and the cavity field such that 
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With this simplification, we find that the Hamiltonian ( 11801 ) 
reduces to 



^ n— 1 

+ /3,„e-^^="(j(;!)^ai. + jiltat_)+H.c.}. (184) 

This equation is in the form of a non-RWA Hamiltonian of 
two extended atomic ensembles independently coupled to two 
counter-propagating cavity modes. The parameters of the 
Hamiltonian are a function of the detunings and Rabi fre- 
quencies of the two highly detuned laser fields co-propagating 
with the clockwise cavity mode and thus could be controlled 
through the laser frequencies and intensities. In the small 
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sample case of kxjn ^ 1 and under the single-mode ap- 
proximation, the Hamiltonian simplifies to the standard Dicke 
model llislQl. 



In the final step we shall reformulate the Hamiltonian ( 1184b in 
terms of bosonic variables by adopting the Holstein-Primakoff 
representation of angular momentum operators |70]. In this 
representation, the collective atomic operators, J^^IK J^nl 

(n) 

and Jz are expressed in terms of annihilation and creation 



operators C^^l and C^jl' of a single bosonic mode. Pro- 
vided the atoms in each ensemble are initially prepared in 
their ground states { |0j„) }, and taking into account that due to 
large detunings of the driving fields, the excitation probabil- 
ity of each atom is low during the laser-atom-cavity coupling. 
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— 1/2, the collective atomic operators can be 
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are collective bosonic operators with the operators 
and foj"^^ obeying the standard bosonic commutation relation 

ri,(") j,(™)ti _ A A 
[Oj ,0^ \ — OjiOnm- 

Note that the collective bosonic operators do not in general 
commute, i.e. 
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Since, in the Hmit of iV > 1: 
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we obtain 
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trapped inside a ring cavity can be expressed now as the in- 
teraction between the cavity modes and three orthogonal field 
modes; a collective mode c'^ solely coupled to the cavity 
mode a+, which co-propagates with the driving lasers, and 

two modes c'^^^k solely coupled to the cavity counter- 

propagating mode a_ . 



The Hamiltonian ( 1190b holds for the laser fields co- 
propagating with the clockwise mode only. Following the 
same procedure as above, we can easily show that in the case 
of the driving fields co-propagating with the anti-clockwise 
mode a_, the effective Hamiltonian takes the form 
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We see the complete symmetry between the two cases that re- 
versing the direction of the propagation of the laser fields from 
clockwise to anti-clockwise is equivalent to the exchange of 

a+ ^ a_ and k 



-kin the Hamiltonian ( 1190b . 
Our objective is to prepare the atomic ensembles in a desired 
entangled state. To achieve it, we consider the evolution of 



the system under the effective Hamiltonian ( 1190b including 
also a possible loss of photons due to the damping of the cav- 
ity mode. This is the only damping which we will consider as 
we have already eliminated spontaneous emission by choos- 
ing large detunings of the driving lasers. 
With the cavity damping included, the properties of the sys- 
tem are determined by the density operator p whose the time 
evolution is governed by the master equation 



where 



CcP = X! (^^iPal - alttip ■ 



pal a, 



(192) 



(193) 



is the Liuivilian operator representing the damping of the cav- 
ity field modes with the rate k. 

In what follows, demonstrate how to generate on demand 
multi-mode entangled states in ensembles of cold atoms lo- 
cated inside a two-mode ring cavity. 



which shows that in the limit of small separations between the 
atoms, the collective bosonic operators are orthogonal to each 
other. 

In terms of the collective bosonic operators cj^^} and C^^k^ 
the effective Hamiltonian (1184b can be written as 



+ 



((iune-''^-C^2k + A„e-'^^=" Ci^t) al +H.c.} 
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The important property of the bosonic representation is the 
fact that effective Hamiltonian of ensembles of cold atoms 



B. One and two-mode entangled states 

First we illustrate a procedure which constructs entangled 
states in a single (M = 1) ensemble of cold atoms located 
inside a two-mode ring cavity. We focus on the creation of 
single and two-mode entangled states and attempt to charac- 
terize the entanglement in terms of unitary operators called 
single and two-mod e squ eezed operators. The squeezed oper- 
ators are defined as ill 14ll 
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exp 



S±k{ii) = exp(^CiC2kC^2k~^iClcU,) , (194) 
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where and are complex one and two-mode squeezing 

parameters, respectively. 

We adopt the procedure of Li ill 13ll . which constructs 
squeezed states from the vacuum by a unitary transformation 
associated with the realistic dynamical process determined by 
the master equation (|190t . In the procedure the atomic en- 
sembles, initially in the ground state, are subjected to different 
sequences of laser pulses of specifically chosen Rabi frequen- 
cies and fisn and phases and (f>sn- The laser pulses 
prepare the ensembles in a desired state that then decays to a 
steady-state with the rate k. We will assume in all our consid- 
erations that the initial ground state of the atomic ensembles 
corresponds to all the atoms being in their ground states 1 j„ ) . 
The squeezed operators can be easily associated with evolu- 
tion operators for the effective Hamiltonians ( 1190b and (1191b . 
which are already expressed in terms of the collective bosonic 
operators Cg^^ and C^2k- squeezing parameters are then 
given in terms of the coupling strengths /?„„ and /3s„ and the 
phases 0„„ , ipsn, so they can be adjusted and controlled by the 
driving laser fields. 

The construction procedure is done in two steps. In the first 
step, we adjust the driving lasers to propagate in the clockwise 
direction, along the cavity mode a+. In this case, the dynam- 
ics of the system are determined by the Hamiltonian (1190b . 
We then send series of laser pulses of phases (pui = 4>si = 
and arbitrary Rabi frequencies ilui and ilsi, but such that 
f^ui > Psi- With this choice of the parameters of the driv- 
ing lasers and under the unitary squeezing transformation 



with 



So{-S.o)S±k{~S.i)pSQ{S,o)S±k{^i) = p, 



the master equation ( |192b becomes 



dt 
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where 



state of the system, we confine our attention only to the eigen- 
values of Eq. (1197b . which are of the form 



77+ = ± 
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Evidently, both eigenvalues have negative real parts which 
means that the system subjected to a series of laser pulses 
up to a short time t will then evolve (decay) to a stationary 
state that is a vacuum state. Thus, as a result of the interaction 
given by the Hamiltonian ( 11981 ). and after a sufficiently long 
evolution time, the modes a±, Cok and C2k will be found in 
the vacuum state, whereas the mode C-2k will remain in an 
undetermined state. The state of the mode C-2k will be de- 
termined in the next, second step of the preparation process. 
In order to estimate the time scale for the system to reach 
the steady-state, we see from Eq. ( 11991 ) that as long as 
V l^ui ~ f^si > th6 time scale for the system to reach 
the steady state is of order of ^ 2/k. Thus, as a result of the 
cavity damping the system, after a sufficient long time, will 
definitely be found in the stationary state. 
In summary of the first step of the preparation, we find that 
in the steady-state, the density matrix representing the state of 
the transformed system is in the factorized form 



p(r ~ 2/k) = /5„ (?) pc_ 
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where 



Pv = |0a+ , Oa_ , Oco. , Oc,J (Oa+ , Oa_ , Oco. , Oc.. I (201) 

is the density matrix of the four modes prepared in their vac- 
uum states, and pc^^k 1^ density matrix of the mode C-2k 
whose the state has not been determined in the first step of the 
procedure. The ket |0a^ , Oq , Qcak i ^C2k) represents the state 
with zero photons in each of the modes. 
Thus, we are left with the problem of the preparation of the 
remaining collective mode C_2fc in a desired squeezed vac- 
uum state. This is done in what we call the second step of the 
preparation, in which we first adjust the driving lasers to prop- 
agate along the anti-clockwise mode a_. We then send series 
of pulses of frequencies, phases and amplitudes the same as in 
the above first stage. As a result of the coupling to the cavity 
mode a_, the interaction is now governed by the Hamiltonian 
( 11931 ), and therefore after the unitary squeezing transformation 
the Hamiltonian of the system takes the form 



^2i-/3?i(«+^ofc + a^-C2fc+H.c.) . (198) = S^{-^^)S±k{-S.i)HeS^{io)S±k{ii) 



It is seen that under the squeezing transformation, the Hamil- 
tonian represents a simple system of two independent linear 
mixers, where the collective bosonic modes Cofe and C2k lin- 
early couple to the cavity modes aj^ and a_, respectively. The 
mode C-2k is decoupled from the cavity modes and therefore 
does not evolve. In other words, the state of the mode C_2fe 
cannot be determined by the evolution operator for the Hamil- 
tonian ( |198b . The important property of the transformed sys- 
tem is that the master equation ( 1 197b is fully soluble, i.e. all 
eigenvectors and eigenvalues can be obtained exactly. Hence, 
we can monitor the evolution of the bosonic modes towards 
their steady-state values. Since we are interested in the steady- 



= ^ Pli-Pli{alCok+alC-2k + ^-c) . (202) 

As above in the case of the coupling to the cavity mode a+, the 
Hamiltonian ( 1202b describes a system of two independent lin- 
ear mixers. Hence, the state of the system will evolve during 
the interaction towards its stationary value, and after s suitably 
long time, ^ 2/k, the transformed system will be found in the 
vacuum state. 

Thus, after the second step of the preparation, the transformed 
system is found in the pure vacuum state determined by the 
density matrix of the form 



(203) 
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where 

- |0a+,0a_,0co.,0c,„0c_,J (204) 

represents the vacuum state of the transformed system and the 
ket 1^1') represents the final stationary state of the system. 
If we now perform the inverse transformation from |^) to j^"), 
we find that the system is in the muhi-mode pure squeezed 
state 

I*) = 5o(eo)5±fc(ei)|0co.,0c.„0c_,,)®|0,+ ,0a_).(205) 
The density operator representing the multi-mode squeezed 



state (1205b is 

P = 5o(eo)5±fc(ei)|{0})({0}|5o(-Co)^±fc(-a), (206) 

where |{0}) = |0c„.,0c..,0c_.J® |0a+,0a_>. 
It shows that in the steady-state the cavity modes are left in 
the vacuum state and the atomic ensemble is prepared in sin- 
gle and two-mode squeezed states. In other words, we have 
found that the collective mode Cq/c is prepared in the one- 
mode squeezed vacuum state 5'o(Co)|Ocofc)' whereas the col- 
lective modes C±2k are in the two-mode squeezed vacuum 
state <S'±fe(^i)|0c2fc, Oc_2fc) associated with the superposition 
of two, position dependent counter-propagating modes. It is 
interesting to note that only a single set of the laser parame- 
ters, the Rabi frequencies and phases, was required to create 
the pure squeezed states. 

The most interesting of the results of this section is that the 



state ( 1205b cannot be factorized into the product of states of 
the individual modes. Thus, we may conclude that the atomic 
ensemble, after the interaction with the sequences of the laser 
pulses is prepared in single and two-mode entangled states. 



C. Four-mode entangled state 

Thus far our discussion has been limited to creation of en- 
tangled states inside a single atomic ensemble. It is partic- 
ularly interesting to generate entangled states between two 
(M = 2) atomic ensembles. Since the cavity is composed of 
two counter-propagating modes, the interaction between the 
collective modes and the pair of cavity modes can produce 
squeezing of the multi-mode variety, with a possibility to cre- 
ate a four- mode squeezed state 1,1 isl ll 16ll . 
We now proceed to illustrate how it is possible to create a 
four-mode squeezed (entangled) state between two atomic en- 
sembles by a separate addressing of the collective modes with 
sequences of laser pulses of suitably chosen Rabi frequencies 
and phases. We again adopt the procedure of Li and demon- 
strate how to prepare the four collective bosonic modes C^^k 
in the following pure four-mode squeezed state 

I*) =^±fc(O|0^(i),0„(i, ,0^(2,,0„(2, ), (207) 



where the unitary operator has the form 

(c\ - o^T^/ f /r^(i)tr^(i)t , ^(i)t^(2: 
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and ^ is the squeezing parameter which we take to be a real 
number. As in the above cases of one and two-mode squeez- 
ing, the squeezing parameter can be adjusted and controlled 
by the driving laser fields. 

Before introducing the procedure, we would like to point out 
that a transformation of the density matrix of the system with 
the unitary operator (1208b does not lead to a density ma- 
trix describing a simple system of independent linear mix- 
ers. Therefore, we first transform the state j^*) to a new state 
1$) — T\'^), with a unitary operator T chosen such that it 

in) 

transforms the field operators C)^2k operators 
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which are linear combinations involving the operators of dif- 
ferent ensembles only, i.e. 



_ '-'2k +^'-'2fc 
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where A = (1 + %/5)/2. 

With this carefully chosen transformation, we find that the 
state 1$) takes a form 



1 



.(2)tr^(l)t 
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where S^l{X^) and S^^l{—^/X) are squeezing operators in- 
volving field modes of the ensemble 1 and 2, respectively. 
As we shall see, the advantage of working with the trans- 
formed state 1$) rather than the state \'^) is that the den- 
sity matrix transformed with the operators S^l{X£,) and 

(2) 

S]_l{~-£_/ X) describes a simple system of a linear mixer of 
a cavity mode with one of the field modes. 
We now proceed to perform the construction of the four-mode 
squeezed state (1207b . which is done in four steps, each sepa- 
rately addressing one of the collective modes. Since the col- 
lective modes are orthogonal to each other, an arbitrary trans- 
formation performed on one of the modes will not affect the 
remaining modes. 

In the first step, we adjust the driving lasers to propagate in the 
direction of the clockwise mode a+. In this case, the dynamics 
of the system are described by the Hamiltonian ( 11901 ), and 
after a unitary transformation of the density matrix 

4'i(e/A)5W(-Ae)TpTt4i)(A05i'i(-e/A) = Pi(212) 
with the squeezing parameter 



^ 2X ""^A,, 
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and with the laser Rabi frequencies and phases such that 



4'u„ = 0s„ = 0, 71 = 1, 2, 



with 

H2 



(214) 
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we find that the master equation of the transformed density 
matrix is of the form 



H.c. 
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dt 
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where 
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We see that the choice of the laser parameters ( 1214b results in 
the collective mode c'^^ , out of the four modes involved in 
the squeezing operator ( 1208b . being effectively coupled to the 
cavity modes, the other modes are decoupled. It may appear 
surprising that only one of the four collective modes is effec- 
tively coupled to the cavity mode. However, this is merely a 
consequence of the form of the effective interaction Hamilto- 
nian ( |190b . It has an obvious advantage that we can separately 
address each of the collective modes. One may also notice 
from Eq. ( 1216b that apart form the collective mode Cj^^ the 



We see that the choice of the laser parameters, Eq. (1217b . has 
resulted in the collective mode C_2k effectively coupled 
to the cavity modes, with the other modes decoupled. The 
dynamics of the mode are determined by the master equa- 
tion ( |219b . which is of the similar form as Eq. ( |197b . Thus, 
we can use the same arguments as in the previous step, and 
conclude that after a sufficiently long time, ^ 2/k, also the 

mode C_2j: will be found in a stationary vacuum state. 
In the next step, we change the direction of propagation of the 
driving lasers back to the direction of the clockwise mode a+. 
Similarly as in the above two steps, we first perform the uni- 



tary transformation of the density matrix, Eq. (1212b . with the 
laser parameters chosen such that 



and 



in) 

modes Cqj. are also coupled to one of the cavity modes. Since 

in) 

the modes Cgj. are not involved in the multi-mode squeezed 
state (1207b whose construction we are interested in, we do not 
consider their evolution. 

Fortunately, the master equation ( 12151 ) is of a similar form 



= 0S1 = 0, = 



A / fiu2 + fis^ 
^ 2 ""^/J^, 
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as the master equation ( 1197b . Therefore, we can foUow the 

same arguments as before to conclude that the mode C^J? will collective mode C'^^^ will be found in a stationaiy vacuum 



With this choice of the laser parameters (1221b . we find that 
now the collective mode C^jJ is the only mode coupled to the 
cavity modes and the dynamics of the mode follows the same 
pattern as the modes considered in the above two steps. Thus, 
we may conclude that after a sufficiently long time, also the 

(2) 



■'2k 

evolve in time and definitely after a sufficiently long time, 
^ 2/k, the mode will be found in a stationary vacuum state. 
The other three modes will remain in undetermined states. 
We now turn off the lasers propagating in the direction of 
the clockwise mode, and perform the second step, in which 
we first adjust the driving lasers to propagate in the direction 
of the anti-clockwise mode a_. In this case, we adjust the 
squeezing parameters such that 



state. 

Finally, in the fourth step. 



we prepare the remaining 



mode ci^2fc the vacuum state. To do this, we again change 
the direction of the driving lasers to propagate in the direction 
of the anti-clockwise mode and choose the laser parameters 
such that 



TT, 



4>u„ = </'s„ = 0, 71 = 1, 2, 
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and 



to obtain the same value for the squeezing parameter now 
given by 
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For this choice of the parameters, and bearing in mind that 
the dynamics of the system are now governed by the Hamil- 
tonian ( |191b . we find that after a unitary transformation, 
Eq. ( 1212b . the master equation of the transformed density ma- 
trix takes the form 



Following the same procedure as in the previous three steps, 
one can easily show that the choice of the laser parameters, 
Eq. (1223b . results in a master equation for the transformed 
density operator determined by the interaction Hamiltonian 
involving the collective mode cf^jfe- Thus, as a result of the 
damping of the cavity modes, the collective mode will evolve 
towards a vacuum state. 

In this way, the final state of the transformed system is a four- 
mode vacuum state determined by the density operator 



dt 



P2 = -i[H2,P2] + CcP2, 



(219) 



p(r^8/«;) = |$)($|, 



(225) 
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where 
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represents the vacuum state of the transformed system and the 
ket I represents the final stationary state of the system. 
Note that the density matrix ( 12251 ) cannot be factorized into 
the product of the density matrices of the individual modes. 
Thus, we may conclude that the atomic ensemble, after the 
interaction with the sequences of the laser pulses is prepared 
in four-mode entangled states. 

Finally it may be mentioned that the procedures which con- 
struct entangled states by using only simple unitary transfor- 
mations offer further interesting prospects for the study of 
multi-mode entanglement and creation of continuous variable 
entangled cluster states. Examples of cluster states and a pro- 
cedure of their creation are discussed in the next section. 



Xm. CLUSTER STATES OF ATOMIC ENSEMBLES 

As the last issue of this review, we wish to make a further the- 
oretical application of the bosonic continuous variable model, 
introduced in the preceding section, to discuss a practical 
scheme for creation of cluster states of separate atomic en- 
sembles. Cluster states have been recognized as suitable for 
one-way quantum computation and possess a high persistency 
and robustne ss of their entanglement with regard to decoher- 
ence effects ill ITlll isll . There have been many proposals for 
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FIG. 3 1 : Configuration of a ring cavity and four atomic ensembles 
for the preparation of entangled continuous variable cluster states. 
External pulse lasers of the Rabi frequencies f^un and Q,sn couple to 
only one (propagating) mode of the cavity and interact dispersively 
with the atoms. 



icant interest is now in the possibility of creating continuous 
variable cluster states of distinct atomic ensembles. 
We present here the results of the study of a practical scheme 
for creation of continuous variable entangled cluster states of 
four separate atomic ensembles located inside a hig h-Q r ing 
cavity. Specifically, we employ the method of Li ill 13ll to 
demonstrate how to prepare the so-called continuous variable 
linear, square, and T-type cluster states of the atomic ensem- 
bles using the sequence of laser pulses and the cavity dissipa- 
tion. 

We shall consider a system consisting of four atomic ensem- 
bles located inside a high-finesse ring cavity. The cavity is 
composed of four mirrors that create two modes, called prop- 
agating and counter-propagating modes, to which the atomic 
ensembles are equally coupled. External pulse lasers that are 
used to drive the atomic ensembles couple to only a single 
propagating mode, as it is illustrated in Fig. |3T| We assume 
that the atoms are homogeneously distributed inside the en- 
sembles. In this case only the forward scattering occurs that 
allows us to neglect the coupling of the atoms to the counter- 
propagating mode and work in the single mode approxima- 
tion. In other words, the cavity mode propagates with the laser 
fields. The effective Hamiltonian of the system, Eq. ( |190l l, can 
then be simplified to the form 



He. = 



H.c. 
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where C„ = C, 
The condition of the coupling of the atomic ensembles to only 
a single mode of the cavity can be easily fulfilled in trapped 
room-temperature atomic ensembles where fast atomic os- 
cillations over the interaction time lead to a collectively en- 
hanced coupling of the atoms to a singl e mode that is practi- 
cally collinear with the laser fields [123]. 
We now proceed to discuss the detailed procedure of the 
preparation of continuous variable four-mode cluster states of 
atomic ensembles located inside a single-mode ring cavity. In 
particular, we shall show how to deterministically prepare a 
linear cluster state, a square cluster state, and a T-shape clus- 
ter state, as shown in Fig.[32l 

We begin with a brief explanation of how one could dis- 
tinguish that a given state belongs to the class of cluster 
states. Simply, a given state is quantified as a cluster state 
if the quadrature correlations are such that in the limit of in- 
finite squeezing, the state becomes zero eigenstate of a set of 
quadrature combinations [1241 
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creation of cluster states that involve squeezed light sources, 
a network of beam splitters and a measurement by the homo- 
dyne detection [119-122]. Essentially all these proposals are 
based on linear optics schemes and measurements on the sys- 
tem. Apart from the schemes based on linear optics, a signif- 



where x and p are the position and momentum operators 
(quadratures) of a mode a, and the modes h are the nearest 
neighbors Na of the mode a. In what follows, we quantify a 
given state as a cluster state by evaluating the variances of lin- 
ear combinations of the momentum and position operators of 
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the involved field modes. If the variances vanish in the limit of 
the infinite squeezing then, according to the above definition, 
a given state is an example of cluster state. 



A. Preparation of a linear cluster state 

Let us first discuss in greater detail the problem of the prepara- 
tion of a four-mode linear cluster state involving four separate 
atomic ensembles, as shown in Fig. [32l a). The essential idea 
is as follows. We begin by making an unitary transformation 

dL„=TiCnTl n= 1,2,3,4 (229) 

that transfers the C„ operators into linear combinations 
1 



V2' 
1 



dl = y={zCi+C2) 

ds = - — {Cs + tCi) 

- 



{iCi -C2- 2iC3 ~ 2C4), 



1 

1 



(2Ci + 2iC2 + C3- iCi 



(230) 



Since the operators (i„ commutate with each other, the com- 
bined modes are orthogonal to each other In this case, the 
system consists of completely decoupled modes, which will 
allow us to prepare each mode separately in a desired state. 
In other words, an arbitrary transformation performed on the 
operators of a given mode will not affect the remaining modes. 
Next, we prepare the field modes in a desired state by using 
laser pulses of equal length and suitably chosen magnitudes of 
the Rabi frequencies and phases. More concretely, for a clus- 
ter state, all the modes dn should be prepared in a squeezed 
vacuum state. We employ the fact that the modes can be sep- 
arately prepared in the squeezed vacuum state. Afterwards, if 



C — C — © — €) 

(a) 



f) 



^ — €) 

(b) 



FIG. 32: Examples of continuous variable four-mode cluster states: 
(a) linear cluster state, (b) square cluster state, and (c) T-shape cluster 
state. 



the variances of specific combinations of the quadrature com- 
ponents of the field operators vanish in the limit of the infi- 
nite squeezing, the desired state is a linear continuous variable 
cluster state. 

In the first step of the preparation, we send a set of laser pulses 
driving the the atomic ensembles 1 and 2 only, the lasers driv- 
ing the atomic ensembles 3 and 4 are turned off. We choose 
the Rabi frequencies and phases of the turned on lasers as 

^un = — =72^!, n = l,2, 
r 

3 1 

0nl = 't'sl = 2^^' '^"2 = '/>s2 = TT. (231) 

After the first step of the preparation, the mode di is left in a 
state that is described by the density operator 



which obeys a master equation 



dtP' 



^di +ra^d\] +H.c.,pi 



(232) 



LaPi, (233) 



where /3 = 2(3^(1,2) = 2^s(i,2), and r e (0, 1). 

If we now perform the single-mode squeezing transformation 

for the mode di as 



Pi^S\{£,)p^S,{0. 
with the single-mode squeezing operator lll25ll 



S'„(^)=exp 



{dl - 4^) 



(234) 



(235) 



where = 
becomes 



tanh (r), we find that the master equation ( 1233 



dt 



Pi 



ad\ + a^di, 



Pi 



LaPi. (236) 



It can be shown from Eq. (1236b that in the steady state, the 
cavity mode will be in the vacuum state and the mode di will 
be in a squeezed vacuum state. This is easy to understand, 
the master equation ( 1236b represents two coupled modes with 
the cavity mode linearly damped with the rate k. Since the 
remaining modes ^2,(^3, and d^ are decoupled from the mode 
dl, they remain in an undetermined state determined by the 
density operator Pdadadil^)- 

In order to estimate the required time to reach the steady state, 
we calculate eigenvalues of Eq. ( 12361 ) 



± 



(237) 



from which we observe that as long as — r"^ > n/2, 
the time for the system to reach its steady state is of order of 
^ 2/k. Therefore, the system will definitely evolve into the 
steady state, provided the interaction time is sufficient long, 
for example, t — A/k. The time t = A/k determines the time 
scale in our protocol for the preparation of the cluster states. 
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By taking the inverse unitary transformation, it follows that in 
the steady state, the total system is in a state determined by 
the density operator 

Pl(T)=5i(O|0a,0<i,^)(0a,0rf,J5l(O®Pd,,d,3d,,(T).(238) 

Briefly summarize what we have obtained after the first 
step of the preparation. The application of suitably chosen 
laser pulses leaves the mode di prepared in the single-mode 
squeezed vacuum state, with the cavity field found in the vac- 
uum state, and the remaining combined modes d2, d^, and d^ 
left in the states related to their initial states. 
In the second step, we turn off the first series of driving 
lasers and sequentially send another series of laser pulses with 
different parameters to preparation of the combined bosonic 
mode d2 in the single-mode squeezed vacuum state so that a 
similar linearly mixing interaction between the cavity mode 
and another combined bosonic modes arises. As before, for 
the first series of pulses, a single-mode squeezed vacuum state 
for this combined bosonic mode can be prepared due to the 
cavity dissipation. The Rabi frequencies of the second series 
of the laser pulses, which are turned on during the time of 
t e [r, 2t) for the preparation of the combined bosonic mode 
d2 in the single-mode squeezed vacuum state 5*2(^)10^2), are 



n, TO = 3, 4, (239) 



and the phases of the driving lasers 



t>u2 



^3= ^ 

i>s2 = (j^ui = (f>s4 = 0. 



(240) 



The specific choice of the Rabi frequencies and the phases 
ensures the mode d2 to be prepared in the state described by 



the density operator pi satisfying the same master as Eq. ( 1236b 
but with di replaced by d2- 

In the third step, which is performed during the time of t e 
[2t, 3t), the combined bosonic mode ds is being prepared in 
the single-mode squeezed vacuum state S'3(^)|0ci3). We turn 
off lasers driving ensembles 1 and 2, and the laser driving the 
ensembles 3 and 4 have Rabi frequencies and phases 



y«3 



= \/2rj, m = 3,4, 
1 



's3 



0u4 = (Psi = TT. 



(241) 



The fourth series, laser pulses are turned on during the time 
of t Q [3t, 4t). In this final series, the combined bosonic 
mode d4 is prepared in the single-mode squeezed vacuum 
state 54(^)|0d^j). The laser pulses required to achieve this are 
of the Rabi frequencies 



TO = 3, 4, (242) 



and phases 



y«2 



-TT, <ps2 = -TT, 0„4 = -TT, (p^i = -^T^ , 

c^sm^O, m = l,3. (243) 



Therefore after enough long time 4r, the system evolves into 
the state determined by the density operator 



pi(4T) = |$i)(<l>i|0|O,)(O,|, 



(244) 



where 



®'53(0|Od3)'^^4(0|Od,). (245) 

This means that the transformed state is in the form of a 
quadripartite squeezed vacuum state. Once all the combined 
modes d„ are prepared in the corresponding single-mode 
squeezed vacuum state, the explicit form of the state created 
can be obtained by reversing the unitary transformation Ti. 
This leads to a pure linear continuous variable quadripartite 
state of the form 



exp 



10 



- A{Ci + iC2){C3 - iC4) 

- 8i(CiC2 + C3C4)]-H.c.}|{0}), (246) 



where 



|{0}) = |0e,, 0,2, 0^3,0, 



(247) 



The question remains as to whether the state is an exam- 
ple of the linear continuous variable cluster state. To examine 
this, we introduce the quadrature amplitude and phase com- 
ponents of the four modes 

qn = (Cn + Cl)/V2, = -i{C^ ~ Ci)/V2, (248) 

from which we easily find that the variances of the linear com- 
binations of the components, evaluated according to the defi- 
nition ( 12271) . are 



-2? 



V{p2 - 91 - qz) 



V{P3 - 92 - 94) - 2"^"'*' ^^P^ " '^3) - e-2«, (249) 

where V{X) — {X"^) — (X)^. Clearly, afl the four variances 
tend to zero in the limit of infinite squeezing, ^ — > 00. We 
therefore conclude that the state is a four-mode linear 
continuous variable cluster state. 

We may summarize that by an appropriate driving the four 
atomic ensembles, the four combined bosonic modes will 
be ultimately prepared in four single-mode squeezed vacuum 
states. It can be done in four steps by appropriately choos- 
ing the laser parameters such as phases and intensities. Then, 
by applying an inverse unitary transformation, the pure linear 
continuous variable entangled state for the four atomic 
ensembles can be obtained. 
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B. Preparation of a square cluster state 



the state 



We now turn to the calculation of the realization of an entan- 
gled continuous variable square cluster state \'^s)- We follow 
exactly the same route as we employed in the preparation of 
the continuous variable linear cluster state. 

A square continuous variable cluster state may be obtained by 
performing a unitary transformation /„ = T2CnTl to obtain 
superposition modes determined by the following operators 



(iCi +2C3 + 2C4), 



/l = 
h = 
h = 



(2Ci+2C2+iC3 + iC4), 



(250) 



Similar to the preparation of the linear continuous variable 
cluster state, we send series of laser pulses with different ap- 
propriate parameters, the phases and amplitudes. Each se- 
ries serves to prepare the combined bosonic modes /„ in a 
squeezed vacuum state. 

Table 1 Rabi frequencies of the laser pulses for the preparation of a square 
cluster state 
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The required Rabi frequencies and phases of the laser fields 
to create the square cluster state are listed in Tables 1 and 2. 
Different sets of laser pulses are chosen for different periods 
of time. Throughout this sequence of laser pulses each indi- 
vidual field mode is incoherently damped by the cavity losses 
with the rate k. 

Table 2 Phases of the laser pulses for the preparation of a square cluster 
state 
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After the above sequence of the laser pulses, the system is left 
in a state described by the density operator 



Pl{AT)=T2\^s){-^s\Tl®\Qa){Qa\ 



(251) 



exp 



1 
10 



[Cl 



cl 



cl 



- 4z(Ci + C2)(C3 + C4) 

- 8(CiC2 + C3C4)]-H.c.}|{0}). 



(252) 



Finally, we calculate the variances in the sum and difference 
operators following the general procedure for quantifying the 
cluster states and find that the variances for the square type 
state illustrated in Fig.[3Tl'b) are given by 



y{j>\ - 93 - 94) 



= y{v2 - 93 - 94) 

y(P3 - gi - 92) = y{.VA - 91 - 92) 
3 



-2? 



(253) 



In the limit ^ — > cx), it is easy to see that the variances tend to 
zero. Thus, according to the criterion (12281 ), the state \^s) is 
an example of a four- mode square cluster state. We conclude 
that the sequential application of the laser pulses, the continu- 
ous variable entangled quadripartite square cluster state \^s) 
is unconditionally produced. 



C. Preparation of a T-shape cluster state 

The final example illustrates the procedure to prepare a T- 
shape entangled cluster state of the form 



exp 



i{cl - cl - cl - cl 



2(CiC-i + C2C4 + C3C4) 
2iCi(C2 + C3 + C4)]-H.c.}|{0}). 



(254) 



To arrive at the cluster state ( 1254b . we first perform a unitary 
transformation = T'iCnT\ that transfers the bosonic oper- 
ators Cn into superposition modes 



hi 



V3 
2 

3 

^(C3-C4), 

I {iCi + C2 + C3 + Ci) . 



iC^^^{C2 + C^ + Ci) 



C2 ~ 2 + C4) 



(255) 



Inverting the transformation T2, we find that the system is in 



We will show that the transformation T3 leads to the field 
modes that can be prepared in a T-shape cluster state. 
Following the similar four-step procedure as in the above two 
examples, we first show that by a suitable choice of the Rabi 
frequencies and phases of the laser fields, one can prepare 
each individual mode ft,„ in a squeezed vacuum state. 
To achieve this, we choose a sequence of the laser pulses of 
the Rabi frequencies and phases as listed in the following Ta- 
bles 3 and 4. Of course, throughout this sequence of laser 
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pulses each field mode is damped to a steady-state by the cav- 
ity losses with the rate k. 

Table 3 Rabi frequencies of the laser pulses for the creation of a T-shape 
cluster state 

t f^iil ^sl ^u2 ^s2 f^uS f^sS ^u4 

[o,r) V3n vsn 

IT,ZT) U U —3— —3— -3- —3— -3- — 3— 

[2t,3t) v^Q V2rn V2n V^rU 

[3t,4t) n rfi n r-n n T-n n r-n 



able cluster states of effective bosonic modes realized in 
four physically separated atomic ensembles interacting col- 
lectively with a single-mode optical cavity and driving laser 
fields. We have demonstrated robustness of the scheme on 
three examples of the so-called continuous variable linear, 
square and T-type cluster states. These examples show that it 
should be straightforward in principle to create different types 
of cluster states with only few laser pulses. 



Table 4 Phases of the laser pulses for the creation of a T-shape cluster state 
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After the above sequence of the laser pulses, the system is left 
in a state described by the density operator 



pi(4r) = r3|*T)(*T|rJ ® \0a}{0a\. (256) 

Inverting the transformation T3, we find that the system is in 
the pure state ( 12541 ). 

The only thing left is to determine as to whether the state 
belongs to the class of cluster states. It is done by calculating 
the variances of the sum and difference operators from which 
we find that for the state |^'t), the variances are given by 

V(p2 - 91) V{p3 - qi) 

- l/(p4-gi) = e-2«, 

V{pi-q2-q3-q4) - 2e-2«. (257) 

Clearly, in the limit ^ — > 00, the variances tend to zero. Hence, 
we conclude that the state |5't) is an example of a four-mode 
continuous variable T-shape cluster state. 
In summary of this section, we have discussed a practi- 
cal scheme for the preparation of entangled continuous vari- 



XIV. SUMMARY 

We have shown that the current research on entanglement and 
disentanglement of multi-atom systems covers a wide class of 
situations ranging from isolated atoms located inside high-Q 
cavities to atoms coupled to a common multi-mode electro- 
magnetic field. The results show unexpected behaviors of en- 
tanglement, such as entanglement sudden death, spontaneous 
revival and sudden birth of entanglement, triggered entan- 
glement evolution by imperfection and steered entanglement 
transfer. We have explored the role of the irreversible process 
of spontaneous emission in creation of entanglement and in 
disentanglement of two atoms. We have shown that sponta- 
neous emission does not necessary lead to disentanglement of 
an initial entangled atoms. Under some circumstances, this 
irreversible process can entangle already disentangled atoms. 
Finally, we have presented models for creation of pure con- 
tinuous variable entangled and cluster states between macro- 
scopic atomic ensembles. 
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